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EFFECT OF CARBON ON THE CRITICAL COOLING RATE 
OF HIGH-PURITY IRON-CARBON ALLOYS AND PLAIN 
CARBON STEELS 


By Thomas G. Digges 


ABSTRACT 


This report describes the preparation of a series of high-purity iron-carbon alloys 
and two series of plain carbon steels, the members of each series differing only 
inearbon content. Determinations of the austenitic grain size and critical cooling 
' rates show how carbon affects the hardening properties and permit comparison 

of the hardenability of these alloys and steels. 
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I. INTRODUCTION 


' The composition and grain size of austenite are now recognized as 
factors that influence the critical cooling rate or hardenability of steel. 
In previous investigations on the effect of carbon on the critical cooling 

' tate, carried out at the National Bureau of Standards! with plain 

carbon steels, and by other investigators? with high-purity iron-carbon 

} alloys, due regard was not given to the simultaneous influence of the 

austenitic grain size. At the time the study was made at the National 

Bureau of Standards, it was recognized that wide variation existed in 

the hardening characteristics of different heats of plain carbon steels 

of apparently the same composition, but the marked influence of 
austenitic grain size on the hardenability of steel was not then recog- 
nzed. The experiments described in the present report were made 
to determine the influence of carbon on the critical cooling rate of 
both high-purity iron-carbon alloys and plain carbon steels. The 
method of preparing the alloys and steels used in this study eliminated 

_ the quality factor and variables other than carbon in the composition. 
Determinations were made of the critical cooling rates for both 

far iat eit. Penne hae » Quenching yt Lt — in relation to some quenching media 


. Esser, W. Bilender, and E.§ . 
4 mode erbat Bilender, and E. Spenle, Quenching diagrams of the tron-carbon alloys, Arch. Bisenhiittenw. 
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series of alloys and steels, quenched directly from the temperatupg 
required to establish a constant austenitic grain size with all carbo, & 
in solution. In these experiments, the only variable that influence; 
the critical cooling rate was the carbon content of the austenite, sing 
the time-temperature relationship used was such as to produce com. 
plete solution and uniform distribution of carbon in austenite of 
constant grain size. However, the study was extended to includ 
determinations of the critical cooling rates of the carbon steels quenched 
from the usually recommended hardening temperature range. In the 
latter experiments, the carbon content and grain size of the austenitp 
and total carbon varied with the different steels. 


II. PREPARATION OF MATERIALS 
1. IRON-CARBON ALLOYS 


A heat of high-purity iron was used as the basic material for the 
iron-carbon alloys. This heat was prepared by melting and solidifying F 
electrolytic iron in vacuo in a crucible of chemically pure magnesia 
bonded with a solution of 2 to 3 percent of magnesium chloride in 
water. The ingot, about 1% inches in diameter at the top, 1% inches 
in diameter at the bottom, and 3% inches long, had the surface imper. f 
fections removed before being split longitudinally into two approx. § 
mately equal sections preparatory to mechanical working. All the F 
specimens of the iron-carbon alloys were prepared from this ingot. 

Hot-working was carried out at an initial temperature of 2,100 to 
2,200° F, and the split halves of the original ingot were forged until § 
reduced in thickness to about %. inch and then hot-rolled to bars about F 
¥% inch thick. The scale formed on the surface of the bars during F 
hot-working was removed by grinding, before cold-rolling to the 4 
form of plates, 0.040 inch thick. Prior to carburizing, these plate | 
were cut into strips and annealed at 1,100° F for % hour either in vacw 7 
or in a reducing atmosphere, to remove the effects of cold-working. / 

The advantages of preparing high-purity iron-carbon alloys by 7 
carburization in gases have been pointed out by Mehl and Wells’ § 
They prepared a series of high-purity alloys from hydrogen-purified F 
carbonyl iron by carburization in mixtures of dipentene-hydrogen and F 
dipentene-benzene-hydrogen. These alloys were homogeneous in the F 
carburized condition. ; 

In the present experiments with propane as the carburizer, a study 
was made of the effects of temperature, rate of flow, concentration, F 


direction of flow, and agitation of the gas on the uniformity of carb F ‘ 


rization of the iron. The preliminary carburizing aspen with 
ropane as the carburizer gave unsatisfactory results. However, 
Sensase proved to be entirely satisfactory for uniformly carburuing 
the iron to the desired carbon content. 
Specimens of the high-purity iron were carburized in a mixture 
hydrogen and benzene (less than 0.003 percent of sulfur) within the 
temperature range of 1,670 to 1,700° F, with a gas flow of 90 to 100 
milliliters per minute. The assembly of the apparatus used is shown 
in figure 1. The hydrogen was purified by passing it through Ascanlé 
(B) and silica gel (D) contained in a vertical tube heated electrically to 
1,300 to 1,475° F, and then through magnesium perchlorate (£) 
phosphorus pentoxide (F). The dry hydrogen next passed through 


3 R. F. Mehl and Cyril Wells, Consitution of high-purity iron-carbon alloys, Metals Tech. T. P. 798,4 Gus 
1937). 
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the benzene (G) maintained at a temperature of 60 to 65° F by a water 
| bath (7) in a Dewar flask (H), then through the orifice (K) and 
around the specimen (iM) suspended vertically in a quartz tube by a 
wire of the same material as the specimen, and finally into the air at 
the exhaust end (NV) of the furnace. ‘ 

The desired gas flow was maintained during carburizing by regulat- 
ing the pressure of the system, as indicated by the oil manometer (J), 
with reducing and needle valves attached to the hydrogen tank. 

The furnace, constructed —, to obtain uniform heating of the 
specimens, consisted essentially of an Alundum tube wound with 


| Nichrome wire, with a compensating coil placed at each end of the 
' furnace over the main heating coil. The temperature was adjusted 
and controlled by independent adjustment of the current in each of 
the three branches of the heating circuit. By preliminary tests made 
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Figure 1.—Diagram of apparatus used for carburizing with hydrogen-benzene 
mixture. 


ion, | Prior to the carburizing treatment, the conditions necessary for 


| establishing a uniform temperature were determined, and the current 
Input necessary for maintaining equilibrium at the desired tempera- 
_ ture was ascertained. The procedure for carburizing was as follows: 

After obtaining a uniform temperature in the carburizing zone of the 
_ furnace with the gas flowing at the required rate, the high-purity iron 
| Specimen (0.040 to 0.041 by 0.275 to 0.300 by 8 inches long) cleaned 
with 1G emery paper and washed in xylene, was inserted in the cold 
zone of the quartz tube at the top of the furnace. After remaining 
in this zone 1 to 2 minutes, the specimen was lowered without touching 
the walls of the tube into the carburizing zone, held for a predeter- 
mined time, and returned to the cold zone to cool in a reducing atmos- 
phere to an estimated temperature of 200 to 400° F. It was finally 
allowed to cool to room temperature in air. The time of carburizing 


a the desired carbon content of the alloy, as is shown in 
Te 2. 
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The carburized bars were ground to a width of 0.250 inch, by remoy. 
ing approximately equal amounts from each edge, before being heated 
in vacuo at 1,670 to 1,700° F for 4 to 5 hours to produce uniform 
distribution of the carbon in the cross section. This heat treatment 
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FiacurE 2.—Relation of carbon content to time for carburizing specimens of high ¥ 
purity iron, 0.040 inch thick, at 1,670 to 1,700° F with hydrogen-benzene mizture. © 


Determinations of carbon were made on the bars after the treatment for homogeneity. 


was carried out in a clear quartz tube, with the bar in a horizontal 
position and with each end resting on a support of the same material 
as the bar. The evacuated quartz tube containing the specimen was 
allowed to cool in air from the high temperature to room temperature. 
This cooling was sufficiently rapid to form sorbite (or fine pearlite) in the 
alloys, as is shown for the 0.23- and 1.21-percent carbon alloys in figure 
3. Specimens subsequently used for the determination of the austenitic 
grain size and critical cooling rate had this initial structure of sorbite. 


2. CARBON STEELS 


A forging, 2 inches in diameter, of acid open-hearth steel containing 
0.20 percent of carbon was thestarting material for the preparation of the 
plain carbon steels. The composition of the forging is shown in table|. 


TABLE 1.—Chemical composition of the seal used as the basic material for the carbon 
stee 





Element Element 
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Fiat }.— Structure of tron-carbon alloys after treatment for homogeneity. 


1, 0.23 perce carbon; B, 1.21 percent of carbon. The sorbite, or fine pearlite, is representative of the 
initial stru e of the alloys and steels for the determination of austenitic grain size and critical cooling 
rate, Etcl with 1-pereent nitric acid in alcohol. 500. 
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Figure 4.—Structure of annealed iron-carbon alloys. 


A, 0.23 percent of carbon; B, 0.50 percent of carbon; C, 0.85 percent of carbon; D, 1.21 percent of carbon. 
Specimens annealed at 1,700° F. for 1 hour in vacuo after treatment for homogeneity. The lengtho! 
the micrographs represents practically the entire thickness of the bars. The left border of the micrographs 
shows the structure at the extreme edge of the specimens. Etched with equal parts of nitric acid (1 
percent) and picric acid (5 percent) in aleohol. X100. 
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To minimize the influence of segregation, a %-inch hole was drilled 
along the horizontal axis of this forging before preparing the specimens 
for carburizing. A disk from the cross section of the forging etched 
in hot acid showed no indications of segregation. However, subse- 
quent microscopic examinations of the structures of the quenched 
specimens used in the determination of the critical cooling rates 
Lowed, in some of these specimens, indications of segregation of 
manganese sulfide inclusions or banding of ferrite and carbide. This 
factor evidently was responsible for the wide variation later observed 
in the structures of a few of the specimens prepared and quenched 
under identical conditions. 

Two series of carbon steels were prepared from specimens cut from 
this forging. The first series was prepared by pack-carburizing 
specimens, 0.050 by 0.260 by 8 inches long, in a mixture of approxi- 
mately 85 percent of wood charcoal and 15 percent of barium car- 
bonate at 1,700° F for various periods of time. The carburized 
specimens were then pack-annealed at 1,700° F for 4 to 5 hours before 
normalizing from the usually recommended temperature range. The 
normalizing treatment produced a structure of sorbite in all the 
steels. The bars with this initial structure of sorbite were then 
ground to a thickness of 0.040 inch and a width of 0.250 inch, before 
preparing specimens for the determinations of the austenitic grain 
size and critical cooling rate. 

With the procedure used in preparing the steels by pack-carburizing, 
pack-annealing, and normalizing difficulty was encountered in obtain- 
ing steels of the desired carbon content with the carbon uniformly 


distributed. Furthermore, pack-carburized © orc often show an 


appreciable increase in the oxygen content. Because of these factors, 
the second series of carbon steels was prepared by carburizing speci- 
mens of the 0.20-percent carbon steel in hydrogen and benzene, and 
homogenizing by heating in vacuo. These treatments were carried 
out under the same procedure as that already described for the high- 
purity iron-carbon alloys. Cooling the 0.040-inch-thick bars from 
the temperature used in the homogenizing treatment was sufficiently 
rapid to form sorbite in all the steels. The specimens subsequently 
used for determination of the austenitic grain size and critical cooling 
rate had this initial structure. 


III. MICROSCOPIC EXAMINATION 


Pieces, % to % inch long, were discarded from each end of all the 
bars after the treatments for homogeneity. Microscopic examination 
of the cross section of each end of the bars (after the end pieces were 
removed) showed the carbon to be uniformly distributed. This is 
illustrated for the gas-carburized alloys and steels in figure 4, which 
shows the microstructure of four of the iron-carbon alloys annealed 
at 1,700° F for 1 hour in vacuo after the treatment for homogeneity. 

hese micrographs show the structure of the bar from one surface 
almost to the opposite surface. (The actual cross sections were 0.040 
to 0.041 inch thick.) This annealing treatment produced “abnormal”’ 
structures in all of the hypereutectoid iron-carbon alloys, as shown 
in the 1.21-percent carbon alloy by the coalesced masses of cementite 
partially or entirely surrounded by ferrite. 
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IV. CHEMICAL COMPOSITION 


Spectrochemical and chemical analyses were made on specimens of 
the iron-carbon alloys and the determinations of gas content were made 
by the vacuum-fusion method. The arc spectrum of the alloy wa; 
examined for the sensitive lines of Ag, Al, As, Au, B, Ba, Be, Bi, Q, 
Cb, Cd, Ce, Co, Cr, Cu, Fe, Ga, Ge, Hf, Hg, In, Ir, K, Li, Mg, Mp’ 
Mo, Na, Ni, Os, Pb, Pd, Pt, Rh, Ru, Sb, Se, Si, Sn, Sr, Ta, Th, Ty 
Tl, U, V, W, Y, Zn, and Zr. The elements found present other than 
iron were Ca, Co, Cu, Mg, Mn, Ni, Pb, and Si (table 2). 


TABLE 2.—I mpurities determined in the iron-carbon alloys 














Analysis ! (percent by weight) 
Element —. ff 
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1 Determinations were made on specimens from the bars after carburizing and heat treatment for homo > 
geneity except as follows: Values of nickel, cobalt, and copper by chemical analysis were obtained fromthe > 
electrolytic iron prior to melting, and the values of manganese and phosphorus by chemical analysis wer — 
obtained from the bars as carburized. Spectrochemical analyses were made by B. F. Scribner, chemical — 
analyses by W. H. Jukkola and J. L. Hague, and vacuum-fusion analyses by V. C. F. Holm, all members fF 
of the staff of the National Bureau of Standards. 

2 Spectrum lines of cobalt and manganese were also found. 


Chemical analyses of the electrolytic iron before and after vacuum 
fusion and of the carburized alloys showed no increase in silicon, sul- f 
fur, or manganese in the various operations incident to the preparation | 
of the alloys. Likewise, analyses by the vacuum-fusion method 
showed no difference in the oxygen, nitrogen, or hydrogen contents asa 
result of working or carburizing of the high-purity iron. 

The results of the analyses, given in table 2, show that sulfw, 
nickel, cobalt, and oxygen were the major impurities contained in the 
alloys. These elements amounted to.about 0.021 percent, whereas the 
total percentage (by weight) for all impurities determined was about 
0.030. 

Carbon was also determined in specimens cut from opposite ends 0/ 
the 7%-inch bars. Some of the alloys had the same pes content at 
both ends. Maximum variation in the carbon value determined 1 [ 
any of the alloys was 0.04 percent. This, however, was an extreme [ 
case, and it is believed that the values for carbon for the entire bar | 
are within 0.01 or 0.02 percent of the values shown in the vamous 
figures given in this report. In the alloys prepared later and in some 
of the carbon steels, the carbon was determined in only one section, 
This was adjacent to the specimens used for the determination af 
austenitic grain size and critical cooling rate. 

Determinations of carbon in specimens from opposite ends of eight 
different carbon steels, prepared by carburizing in hydrogen and ben- | 
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Figure 5.—Grain size of iron-carbon alloys at 1,700° F. 


A, 0.23 . . ‘ - 
‘Austenitic grain ave pol compe of carbon; C, 1.01 percent of carbon; D, 1.21 percent of carbon 
4 ot a “ S$ per s r ig ; STM egrai ; i 
with 1-percent nitric acid in alcohol. a inch at 100 diameters (ASTM grain number 4). Etched 
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FIGURE 7.—Grain size of carbon steels at different temperatures, 


A, 0.46 percent of carbon; B, 0.76 percent of carbon; C, 1.10 percent of carbon; D, 1.15 percent of carbon. 
Average grain size of 6 to 12 grains per square inch at 100 diameters at 1,625° F for the 0.46 percent of 
carbon, 1,660° F for the 0.76 percent of carbon, 1,680° F for the 1.10 percent of carbon, and 1,700° F for 
the 1.15-percent carbon steel. Etched with 1-percent nitric acid in aleohol. 100. 
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zene, showed a maximum difference of 0.03 percent in any one steel. 
Variations in carbon of this same magnitude were also shown in the 
annealed bars of the selected steels, which were prepared by pack- 
carburizing. Further minor changes occurred in the carbon content 
of some of the latter steels during the normalizing treatment. The 
values determined after this treatment were used in preparing the 
yarious diagrams for this series of steels. 

Complete analyses were not made on the carbon steels after the 
carburizing and homogenizing treatments. The proportions of the 
elements other than carbon, given in table 1, should not be materially 
affected by the various steps incident to the preparation of these 
steels. Possibly, the oxygen content varied for steels prepared by 
pack-carburizing. Variations in the oxygen content, however, had 
no detectable influence on the critical cooling rate of the steels. 


V. AUSTENITIC GRAIN SIZE 


The specimens for the determination of the austenitic grain size 
had the same initial structure and were of the same shape and size as 
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FiaurEe 6.—Relation of temperature to carbon content of the steels for establishing 
austenitic grain size of 6 to 12 grains per square inch at 100 diameters. 


those used for the determination of the critical cooling rate. They 
were heated in the same manner, held for the same time at the same 
temperature, and were then cooled in a manner suitable for outlining 
the grains. 

The austenitic grain size at 1,700° F was determined for the highest- 
carbon alloy (1.21 percent of carbon). This temperature was suffi- 
ciently high to dissolve all the carbon. Subsequent tests on the lower- 
carbon alloys showed that variations in carbon did not affect the grain 
size at this selected temperature of 1,700° F. 

Each of the iron-carbon alloys had an average grain size of 6 to 12 
grains per square inch at 100 diameters (ASTM grain number 4) at 
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1,700°F. Figure 5 shows the’grain'size forgthe,0.23-, 0.62-,"1.01-, and 
1.21-percent carbon alloys at 1,700° F. The grain size at 1,700° F 
was also determined for the 1.18-percent carbon steel. As shown jn 
figure 6, lower Retnparetaines were necessary to establish the same 
average grain size of 6 to 12 grains per square inch at 100 diameter 
in the lower carbon steels. Figure 7 shows that this grain size wag 
obtained at 1,625° F for the 0.46-percent carbon steel, at 1,660° F for 
the 0.76-percent carbon steel, at 1,680° F for the 1.10-percent carbon 
steel, and at 1,700° F for the 1.15-percent carbon steel. 

On heating within the usually recommended hardening temperature 
range, specimens of the hypoeutectoid carbon steels had approxi- 
mately the same average grain size of 12 to 24 grains per square inch 
at 100 diarheters (ASTM grain number 5). 

No differences were observed between the temperature—grain-size- 
critical-cooling-rate relations of the carbon steels prepared by pack- 
carburizing and those prepared by carburizing in hydrocarbon. 


VI. CRITICAL COOLING RATE 


The work by Davenport and Bain ‘* showed that the hardenability 
of a steel depends upon the stability of its austenite on cooling 
through the temperature range of 
ZZ gage chromel-alumel pa cel F. agetle gr kya 
‘ mperature range in which the aus- 
\ er tenite is most oly to decompose, 
and in which its decomposition prod- 
ucts are relatively soft. If, by 
tie A—t----—-t—A a continuous quench to room tem- 
perature, the austenite is cooled 
| through this temperature range with- 
a i ae r a out decomposing, it transforms only 
040° z in the range below about 300° f 
Ficgure 8.—Form and dimensions of and the product of dec omposition 
the specimens for the determination 18 martensite, and in this condition 
of critical cooling rate. the steel is fully hardened. At some 
cooling rate, called the critical cooling 
rate, the austenite just begins decomposing in the range 1,110 to 930° 
F. Obviously, the critical cooling rate of a steel is an index of the 
stability of the austenite in this temperature range. The slower the 
critical cooling rate (or the deeper the hardening) the more stable is 
its austenite in the temperature range 1,110 to 930° F. In the present 
experiments, the critical cooling rate was taken as the average cooling 
rate between 1,110 and 930° F, which produced in the quenched speci- 
men a structure of martensite with nodular troostite (fine pearlite) in 
amounts estimated to be between 1 and 3 percent. 

For the determination of the critical cooling rate, specimens 0.25 
inch square by 0.040 inch thick were quenched directly from the 
temperatures establishing the austenitic grain size. One wire of 
22-gage Chromel-Alumel thermocouple was spot-welded to the center 
of the flat face of the specimen and the other wire was spot-welded to 
the center of the opposite face (fig. 8). The specimens were heated 
in an atmosphere of nitrogen at a rapid rate, approximately 1 minute 


4 E. S. Davenport and E. C. Bain, Transformation of austenite at constant subcritical temperatures, Trans. 
Am, Inst. Mining Met. Engrs., Iron and Steel Div. 90, 117-154 (1930). 























Digges] Iron-Carbon Alloys and Carbon Steel 579 


being required to reach the desired temperature, held at that temper- 
ature for 15 minutes, and quenched in a bath which sealed the bottom 
end of the vertical tube of the heating furnace. 

To minimize oxidation, the end of the furnace tube extending into 
the quenching bath was either closed by means of a slide, or the bath 
was covered with a layer of quenching oil during the entire time of 
heating the specimens of all the alloys and steels prepared by carburiz- 
ing with hydrogen and benzene. This procedure was modified for 
the steels prepared by pack-carburizing. Specimens of the latter 
steels, quenched from the usually recommended hardening tempera- 
ture range, were heated in an atmosphere of nitrogen, and those 
quenched from temperatures establishing a constant austenitic grain 
size were heated in an atmosphere of nitrogen and illuminating gas. 

The desired cooling rates were obtained by quenching either in 
oil or in a sodium-silicate-water bath of varying concentration. 


IRON-CARBON ALLOYS 
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FicurE 9.—Effect of cooling time on the structure of iron-carbon alloys. 


Specimens 0.25 inch square by 0.040 inch thick, quenched from 1,700° F, the temperature establishing 
austenitic grain size of 6 to 12 grains per square inch at 100 diameters with all carbon in solution. The 
— a the cooling time (resulting in 1 to 3 percent of troostite) and carbon content is represented 

y the straight line. 


The bath was at approximately room temperature at the time of 
quenching. A photographic time-temperature cooling curve was 
obtained during the quench by means of a string-galvanometer 
apparatus described in a previous report.® 

or the correlation of critical cooling rate with structure, micro- 
scopic examinations were made of the quenched specimens on the 
cross section of the 0.040-inch-thick sheet, a short distance from the 
point of contact of the thermocouple wires (fig. 8). 

The critical cooling rate was determined for iron-carbon alloys, 
ranging in carbon from 0.23 to 1.21 percent, and for carbon steels, 
ranging in carbon from 0.43 to 1.18 percent. Each specimen was 
quenched directly from the temperatures establishing a constant 
austenitic grain size with all the carbon in solution. 

The relation between the time to cool from 1,110 to 930° F and the 
carbon content is shown for the iron-carbon alloys in figure 9, for the 


‘H.J. French and 0. Z. ——— Quenchi Sr carbon steels nation § _ 
for heat treatment, Trans. Am. Soc. Steel Treating 251-304 (1924). wera Sis r 
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Figure 10.—Effect of cooling time on the structure of carbon steels. 


Specimens 0.25 inch square by 0.040 inch thick, quenched from temperatures establishing austenitic grain 
size of 6 to 12 grains per square inch at 100 diameters with all carbon in solution. Steels were — 
by carburizing specimens of a 0.20-percent carbon steel in hydrogen-benzene mixture. The 
between critical cooling time (resulting in 1 to 3 percent of troostite) and carbon content is represented 
by the straight line. 
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Ficure 11.—Effect of cooling time on the structure of carbon steels. 
Specimens 0.25 inch square by 0.040 inch thick, quenched from temperatures establishing austenitic grait 
size of 6 to 12 grains per square inch at 100 diameters. Steels were prepared 


by pack 3 
mens of a 0.20-percent carbon steel. The relation between critical cooling time (resulting in 1 to3 percent 
of troostite) and carbon is represented by the straight line. 
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steels pnpenes by carburizing in hydrogen and benzene in figure 10, 
and for the steels prepared by pack-carburizing in figure 11. The 
estimated amounts of nodular troostite formed (1 to 3 percent for 
critical cooling time) with the various quenching rates are also shown 
in the figures, together with the curves representing the relation of 
critical cooling time to the carbon content. 

The critical cooling rate was also determined for the carbon steels, 
varying in carbon from 0.48 to 1.38 percent, quenched from the 
usually recommended hardening temperatures. For these experi- 
ments, the austenitic grain size of the hypoeutectoid steels was 
constant (ASTM grain number 5) and all carbon was in solution. 
For the hypereutectoid steels the carbon content of the austenite 
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Figure 12.—Effect of cooling time on the structure of carbon steels. 


Specimens 0.25 inch square by 0.040 inch thick, quenched from the usually recommended hardening tem- 
perature range. Steels were prepared by either pack-carburizing, or carburizing in hydrogen-benzene 
mixture, specimens of a 0.20-percent carbon steel. The relation between critical cooling time (resulting 
in 1 to 3 percent of troostite) and carbon content is represented by the curve. 


was constant, but not all carbon was in solution and the average 
austenitic grain size was smaller than that of the hypoeutectoid 
steels. The critical cooling rates in these experiments, therefore, 
were influenced by three variables, namely, austenitic grain size, 
carbon content of austenite, and total carbon. 

For the steels quenched from the usually recommended hardening 
temperatures, the relation between the time to cool from 1,110 to 
930° F and the carbon content is shown in figure 12. 

For convenient comparison, the relations of the critical cooling rate 
to carbon content of the alloys and steels are summarized in figure 13. 
For the selected conditions of constant austenitic grain size and solu- 
tion of all carbon in the austenite, the critical cooling rate decreased 
continuously with increase in carbon, for both the high-purity iron- 
carbon alloys (curve A) and the plain carbon steel (curve C). That 
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is, the hypereutectoid alloys had slower critical cooling rates (wen 
deeper hardening) than the alloys of either eutectoid or hypoeutectoiq 
compositions, and the hypereutectoid steels were deeper hardening 
than the steels of eutectoid and hypoeutectoid compositions. Hoy. 
ever, the effect of a small change in carbon content of the austenite of 
the lower-carbon alloys and steels on the critical cooling rate was mor 
marked than a corresponding change in the higher-carbon alloys and 
steels. For example, an increase in carbon from 0.23 to 0.33 percent 
of the alloys decreased the critical cooling rate by about 450° F pe 
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Figure 13.—Relation of critical cooling rate to carbon content of the iron-carbon 
alloys and carbon steels. 


second, whereas an increase in carbon from 1.10 to 1.20 percent de 
creased the critical cooling rate only about 50° F per second. 
_ An important feature brought out in these experiments, therefore, 
is the influence of the carbon content of the austenite on hardenability. 
Since this was the only variable in the riments on critical cooling 
rates of the alloys, the differences observed in their hardenability must 
be attributed to this factor. Likewise, the differences observed in the 
hardenability of the plain carbon steels (curve C) are also due solely to 
differences in carbon content of the austenite. | 
Comparisons of the critical cooling rates of the alloys and steds 
under similar conditions with respect to austenitic grain size 
content, and all carbon in solution (curves A and O, respectively) show 
that the alloys had considerably higher rates, that is, they were com- 
paratively shallow hardening. For any selected carbon content 
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within the range investigated, and with austenitic grain size of 6 to 12 
| grains per square inch at 100 diameters, the ratio of the critical cooling 

rate “4 the alloys to that of the steels was approximately constant, 
and was of the order of 3.5 to 1. : 

The high-purity iron-carbon alloys differ from the plain carbon 
steels in the number of nonmetallic inclusions and amounts of manga- 
nese. ‘The high-purity alloys were relatively free from nonmetallic 
inclusions and manganese. The carbon steels contained a relatively 
large number of nonmetallic inclusions and about 0.46 percent of 
manganese. These factors produce opposite effects on the critical 
cooling rate or hardenability. Nonmetallic inclusions act as nuclei 
for transformation of the austenite on cooling and thus tend to pro- 
mote shallow hardening, whereas elements such as manganese and 
chromium stabilize the austenite or retard its rate of decomposition 
and thus contribute to deep hardening. The manganese content, 
obviously, was the dominant of these two factors in controlling the 
reaction rates and has imparted to the steels comparatively deep- 
hardening properties. 

For the carbon steels quenched from the usually recommended hard- 
| ening temperature range (curve B), the critical cooling rate decreased 
with increase in carbon up to about 0.75 or 0.80 percent; thereafter the 
critical cooling rate increased rapidly with increase in carbon up to 
about 1.05 percent. Further increase in carbon from about 1.05 to 
| 1.40 percent had no appreciable influence on the critical cooling rate. 

Thus, the steel of about eutectoid composition had the minimum 
critical cooling rate, or was deeper hardening, than either the steels of 
hypoeutectoid or hypereutectoid compositions. As already stated, 
these critical cooling rates (curve B) were influenced by variations in 
austenitic grain size and carbon content of austenite, as well as by 
total carbon. 

In these experiments, the average austenitic grain size was approxi- 
mately constant and all carbon was in solution at the time of quench- 
ing the steels ranging in carbon from about 0.40 to 0.75 percent. 
The difference in their hardenability was caused by a difference in 
the carbon content of the austenite. With carbon in excess of about 
0.80 percent, the carbon content of the austenite remained constant 
(quenching temperature of 1,425° F), but differences existed in_the 
grain size, and in the undissolved or total carbon. The undissolved 
carbon or free carbides not only retard grain growth but also (like 
nonmetallic inclusions) act as nuclei for transformation during cool- 
ing. Furthermore, the reaction rate of the austenite is influenced by 
the number, size, and distribution of the free carbides. Assuming 
unformity in distribution and size of the carbide particles, the number 
will vary directly with increase in carbon above the saturation value 
for carbon in austenite at the quenching temperature. With steels 
containing carbon slightly above the eutectoid value, the free carbide 
particles, although relatively few in number, were sufficient to in- 
crease materially the reaction rates. The reaction rates continued to 
increase uniformly and rapidly as the number of free carbide particles 
increased until a value finally was reached (corresponding to a total 
carbon value of about 1.05 percent) beyond which further increase in 
the number of carbide particles had no marked effect on the critical 
cooling rate or hardenability, That is, the rapid increase in critical 
cooling rate with increase in carbon from about 0.80 to 1.05 percent 
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was due principally to variations in the number of free carbide 
particles. 

The marked differences in critical cooling rates of the hypereutectoiq 
steels (curves B and C of fig. 13) were caused by differences in ays. 
tenitic grain size and undissolved carbides at the time of quenching. 
On heating to the usually recommended hardening temperature 
the austenitic grain size was relatively small and some carbop 
remained undissolved. Both these factors contribute to shallow 
hardening. Obviously, to utilize full hardening power of the stegl 
it is necessary to heat the steel to a temperature sufficiently high to 
insure complete solution of the carbon. 

This is further illustrated by comparison of the critical cooling rates 
of the hypereutectoid alloys and steels (curves A and B, respectively) 
As already shown, under similar conditions, the alloys were com- 
paratively shallow hardening. However, the hypereutectoid alloys 
with the larger grain size and with all carbon in solution had slower 
critical cooling rates (were deeper hardening) than the carbon steels 
of similar total carbon content quenched from the recommended 
hardening temperature range. 

The important role of austenitic grain size in controlling harden- 
ability is shown by comparison of the critical cooling rates of the 
hypoeutectoid steels (curves B and C for carbon ranging from 0,40 
to 0.75 percent). For a given carbon content within the range of 
0.40 to 0.75 percent, the steel differed only in the size of the austenite 
grains at the time of quenching. The critical cooling rates represented 
by curve B were obtained with specimens having an average of 12 to 
24 grains per square inch at 100 diameters (ASTM grain number 5) 
and those by curve C, with 6 to 12 grains per square inch at 100 di- 
ameters, (ASTM grain number 4). The smaller grain size steels 
obviously required the higher critical cooling rates. For example, 
increasing the average austenitic grain size from ASTM grain num- 
ber 5 to 4 decreased the critical cooling rate of the 0.65-percent carbon 
steel from about 400 to 300° F per second, or by about 100° F 
second. This change in critical cooling rate is equivalent to that 
resulting from increasing the carbon in solution from 0.45 to 0.65 
percent (constant austenitic grain size of 6 to 12 grains per square 
inch at 100 diameters). 

Attempts were also made to determine the critical cooling rate of 
the basic steel (table 1) used in preparing the carbon steels. For an 
austenitic grain size of 6 to 12 grains per square inch at 100 diameters 
the critical cooling rate of this 0.20-percent carbon steel would be 
expected to be in the neighborhood of 800 to 900° F per second. In 
fact, some of the experimental results showed it to be within this range. 
However, similar tests in which the specimens were cooled at con- 
siderably higher rates, observed values up to 12,000° F per second, 
produced in the quenched specimens some areas containing free ferrite. 
In the latter tests, the ferritic areas in the quenched specimens resulted 
from nonuniform distribution of the carbon in the austenite at the time 
of quenching. Regardless of cooling rates, complete hardening of 
these particular specimens should not be obtained until the ferrite or 
at least the very low carbon areas in the austenite were eliminated. 

A structure containing nodular troostite was not observed in any of 
the quenched specimens of this particular 0.20-percent carbon stedl. 
If transformation of the austenite on cooling occurred at temperatures 
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above the range where martensite was formed, then the transformation 
product was free ferrite, not nodular troostite. With the 0.23-percent 
carbon alloy, the transformation product was nodular troostite. 

Results of investigations of the effect of carbon on critical cooling 
rate have been reported by Esser, Eilender, and Spenle® for high- 
purity iron-carbon alloys, and by French and Kiopsch’ for plain 
carbon steels. In the experiments by the former, the cooling rate 
between 1,470 and 1,290° F was determined by quenching thin wires 
of iron-carbon alloys in gas from about 1,800° F. These alloys were 
prepared from carbonyl iron and carbon. French and Klopsch deter- 
mined the critical cooling rate at 1,330° F for the centers of cylinders 
0.5 inch diameter by 2 inches long quenched from about 135° F above 
Ac; for hypoeutectoid steels and for about 135° F above Ac, for hyper- 
eutectoid steels. = 

They found that the critical cooling rates of both the high-purity 
iron-carbon alloys and plain carbon steels decreased with increase in 
carbon content up to about eutectoid compositions, but with further 
increase in carbon the critical cooling rates also increased. That is, 
when both the iron-carbon alloy and plain carbon steel were of about 
eutectoid compositions, they had the minimum critical cooling rates 
(were deeper hardening) than the hypoeutectoid or hypereutectoid 
alloys and steels. These results were influenced, however, by variables 
other than the carbon content of the materials. These two investiga- 
tions were made without due regard for the influence of austenitic 
grain size and, in addition, the quenching temperatures used by 
French and Klopsch were not sufficiently high to insure complete 
solution of all the carbon in the highest-carbon steel. 


VII. SUMMARY 


High-purity iron-carbon alloys were te from hot- and cold- 
worked specimens of a heat of vacuum-fused electrolytic iron by car- 
burizing in a hydrogen-benzene mixture. The carburized specimens 
were heated and cooled in vacuo to obtain uniform distribution of the 
carbon. Cooling during this homogenizing treatment was sufficiently 
rapid to form sorbite (or fine pearlite) in all of the alloys. Microscopic 
examination of the cross sections of the bars of the different alloys 
showed the carbon to be uniformly distributed. 

The alloys studied contained 0.23, 0.40, 0.50, 0.62, 0.73, 0.80, 0.85, 
1.01, 1.14, and 1.21 percent of carbon, respectively. These values are 
oo to be accurate within the limits of plus or minus 0.01 percent 
of carbon. 

A series of plain carbon steels was also prepared by pack-carburizing 
specimens of a 0.20-percent carbon steel fo owed by pack-annealing 
and normalizing. A second series of plain carbon steels was prepared 
by carburizing specimens of the same 0.20-percent carbon steel in 
hydrogen and benzene and homogenizing by heating and cooling in 
vacuo. 

The plain carbon steels varied in carbon from that of the basic 
material, 0.20 to 1.88 percent. These steels also contained about 0.46 
percent of manganese and 0.15 percent of silicon. 


Py H, Beer, W - Eilender, and E. Spenle, Quenching diagrams of the iron-carbon alloys, Arch. Eisenhiittenw. 
"HJ. French and 0. Z. ~~ Quenching eg for carbon steels in relation to some quenching media 
for heat treatment, Trans, Am, oc. Steel Treating 6, 261-04 (1924). eg 
57763 —38———2 
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Spectrochemical and chemical analyses, as well as determinations of 
gas content by a vacuum-fusion method, showed that sulfur, ni 
cobalt, and oxygen were the major impurities in the iron-carbop 
alloys. The total of all the impurities determined was about 0.039 
percent. 

Abnormal structures were produced in the hypereutectoid alloys 
annealed in vacuo at 1,700° F. 

All specimens used for determination of the austenitic grain size 
and critical cooling rate had an initial structure of An he (fine 
pearlite). 

Austenitic grain size was determined for the highest-carbon alloy 
(1.21 percent of carbon) at a temperature sufficiently high to insure 
complete solution of all the carbon. Subsequent determinations on 
the lower-carbon alloys showed that variations in carbon did not 
affect the grain size at this selected temperature of 1,700° F. Each 
of the alloys at 1,700° F had an average grain size of 6 to 12 grains per 
square inch at 100 diameters (ASTM grain number 4). 

The carbon steel containing 1.18 percent of carbon had the same 
grain size as the alloys at 1,700° F. Lower temperatures were neces- 
sary to establish the same austenitic grain size in the lower-carbon 
steels. On heating within the recommended hardening temperature 
range, specimens of the hypoeutectoid carbon steels had a constant 
average grain size of 12 to 24 grains per square inch at 100 diameters 
(ASTM grain number 5). No differences were observed in temper- 
ature—grain-size-critical-cooling relations of the carbon steels pre- 
pared by pack-carburizing and those carburized in hydrocarbon. 

The critical cooling rate was taken as the average cooling rate, be- 
tween 1,110 and 930° F, which produced in the quenched specimen 
a structure of martensite with nodular troostite (fine pearlite) in 
amounts estimated to be between 1 and 3 percent. 

The critical cooling rate was determined for each of the iron- 
carbon alloys quenched directly from 1,700° F, and for the carbon 
steels varying in carbon from 0.43 to 1.18 percent, quenched also 
from the temperatures establishing a constant austenitic grain size 
with all carbon in solution. 

For the selected conditions of constant austenitic grain size and 
solution of all carbon in the austenite, the critical cooling rate de- 
creased continuously with increase in carbon, for both the high-punity 
iron-carbon alloys and plain carbon steels. That is, the hypereutec- 
toid alloys had slower critical cooling rates (were deeper hardening) 
than the alloys of eutectoid and hypoeutectoid compositions, and the 
hypereutectoid steels were deeper hardening than the steels of eutec- 
toid and hypoeutectoid compositions. However, a small change in 
the carbon content of the austenite of the lower-carbon alloys and 
steels had a more marked effect on the critical cooling rate than did 4 
corresponding change in carbon of the high-carbon alloys and steels. 

The carbon content of the austenite, therefore, had a pronounced 
influence on the hardenability of the alloys and steels. Under similar 
conditions with respect to austenitic grain size, carbon content, and 
with complete solution of the carbon in the austenite, the high-purity 
iron-carbon alloys hardened less deeply than did the plain carbon 
steels. i fall. in hardenability was ine PN to differ- 


ence in manganese content of the alloys and s 
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For the carbon steels quenched from the usually recommended 
hardening temperature, the critical cooling rate decreased with in- 
crease in carbon up to about 0.75 to 0.80 percent, but increased rapidly 
with further increase in carbon up to about 1.05 percent. Changes in 
carbon from about 1.05 to 1.40 percent had no appreciable influence 
on critical cooling rate. These critical cooling rates, however, were 
influenced by variations in austenitic grain size, and carbon content 
of austenite as well as the total carbon content. 

Full depth of hardening was obtained only when all the carbon was 
in solution in the austenite. 


J. A. Bennett, a member of the staff of this Bureau, assisted in 
assembling the apparatus shown in figure 1 and in carrying out the 
carburizing experiments with gases. 


WasHINGTON, February 11, 1938. 
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HEAT OF COMBUSTION OF ISOPRENE 
By Ralph S. Jessup 







ABSTRACT 







The heat of combustion of isoprene has been measured by means of a bomb 
calorimeter which was calibrated with benzoic acid. The mass of isoprene burned 
in each experiment was determined from the mass of CO, formed in combustion. 
As the result of seven measurements the value 3156.9 +1.6 international kilo- 
joules per mole was obtained for the heat of combustion of liquid isoprene at 
25° C and under a constant pressure of 1 atmosphere, the products of combustion 
being liquid water and gaseous CQ:. Using the value 25.9 international kilo- 
joules per mole calculated by Bekkedahl, Wood, and Wojciechowski for the heat 
of vaporization of isoprene, the heat of combustion of gaseous isoprene is calcu- 
lated to be 3182.8 +1.7 international kilojoules per mole at 25° C and under a 
constant pressure of 1 atmosphere, the products of combustion being liquid water 
and gaseous CO.. Combining the data on isoprene with data reported previously 
on the heat of combustion of rubber hydrocarbon, there are obtained for the heats 
evolved in the reactions 


xCsHg(liquid isoprene) = (CsHs).(rubber hydrocarbon) 



















and 






xC;Hs(gaseous isoprene) = (C;Hs) .(rubber hydrocarbon) 






at 25° C and a constant pressure of 1 atmosphere, the values (74.8 +6.2)z and 
(100.7 +6.3)z international kilojoules respectively. 
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I. INTRODUCTION 


_ The work described in this paper was done as part of an investiga- 
tion of thermodynamic properties of rubber hydrocarbon and related 
substances. In previous papers data have been reported on the heat 
of combustion [8]! and on the entropy and free energy of formation 
of rubber a peer [1]. Reventon data on the entropy of iso- 
prene have een reported by Bekkedahl and Wood [2], who discuss 
the chemical similarity of isoprene and rubber hydrocarbon, and 
, Point out the desirability of a knowledge of the free energy of forma- 


‘Numbers in brackets relate to the literature references given at the end of this paper. 
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tion of isoprene for use in studying the thermodynamics of reactions 
involving the two substances. In order to use the entropy data jp 
calculating the free energy of formation, an accurate value for the 
heat of combustion of isoprene is required. 


II. PREVIOUS WORK 


The only data on the heat of combustion of isoprene which haye 
been found in the literature are those of Lebedev, Chochlovkin, ang 
Kalacheva [10]. They report 11 measurements of the heat of com. 
bustion of liquid isoprene, the results of which range from 11,142 to 
11,215 calories per gram, the mean result being 11,165 calories per 
gram. These authors do not give sufficient data to permit an estimate 
of the accuracy of their results, and they do not define their en 
unit, so that accurate conversion of their data to the energy unit 
used for the entropy data is not possible. For these reasons a new 
determination of the heat of combustion of isoprene was desirable. 
If it is assumed that the calorie used by Lebedev, Chochlovkin, and f 
Kalacheva is equal to 4.1833 international joules, their value for the 
heat of combustion of isoprene becomes 46.7 international kilojoules 
per gram, a value which is nearly 1 percent higher than that obtained 
in the present work. 


III. APPARATUS AND METHOD 


The bomb calorimeter and accessory apparatus, experimental pro- 
cedure, and method of calculating results have been deocullae in 
detail [4, 7, 9]. 

Samples of approximately 1 g of isoprene were inclosed in thir- 
walled glass bulbs which were fattened on opposite sides and filled 
completely with the liquid in the manner described by Richards and 
Barry [11]. The samples were ignited by means of an electric fuse 
of iron wire, as described previously [7]. The volume of the bomb 
used was determined to be 377 cm*®. One cm* of water was placed 
in the bomb before each experiment. 

In most experiments the mass of isoprene burned was determined 
in two ways: (1) By weighing the sample of isoprene, and (2) by 
absorbing the CO, formed in combustion in Ascarite (a sodium hy- 
droxide-asbestos mixture) and weighing it, and passing the remaining 
gas through a copper tube heated to redness to oxidize any CO whieh 
might be present, absorbing the resulting CO, in Ascarite and weighing 
[12,7]. The mass CO, obtained by oxidizing the CO was very small, 
usually less than 0.01 percent of the total CO, formed in the bomb. 


IV. CALIBRATION OF CALORIMETER 


The calorimeter was calibrated with benzoic acid (NBS Standard 
Sample 39e), using for the heat of combustion of this material the 
value 26,414 international joules per gram mass under standard 
conditions [9]. These standard conditions are: 

Temperature to which reaction is referred 
Initial absolute pressure of oxygen at 30° C 


Mass of benzoic acid per liter of bomb volume 3 g. 
Mass of water placed in bomb per liter of bomb volume-_- 3 g. 
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The results of the calibration experiments are given in table 1. 
The uncertainty of the mean of the nine experiments (defined as 


9 Vi a iy about three times the “probable error’) is 0.026 
ercent. é ; 

In the calibration experiments the final temperature of the calo- 
rimeter was 30° C, and as the value used for the heat of combustion 
of benzoic acid is that for 30° C, the values obtained for the energy 
equivalent of the calorimeter are for the initial system at the mean 
temperature of the experiments. In other words, the observed values 
| of energy equivalent represent the energy per degree required to raise 
the temperature of the calorimeter from its initial to its final tem- 
' perature (30° C) when the bomb contains 1 g of water, the charge of 
benzoic acid, and oxygen under a pressure of 30 atmospheres at 30° C. 
The values given in table 1 for the energy equivalent of the calo- 
rimeter are the observed values minus the heat capacity of the charge 
of benzoic acid (specific heat=1.26 j/g° C). In each experiment on 
isoprene the final temperature of the calorimeter was very nearly 
30°C, and the value used for the energy equivalent was the mean 
value given in table 1, plus the heat capacity of the charge of isoprene 
(specific heat=2.25 ig C). Since the energy equivalent of the 
calorimeter is that of the initial system the observed values for the 
heat of combustion of the isoprene are for the final temperature of the 
calorimeter, or 30° C. 


TaBLE 1.—Calibration of bomb calorimeter 





Observed- Observed- 
eau Yat | Deviation from mean|| Experiment Pi a Deviation from mean 
at 28.5° C at 28.5° C 


Experiment 





int. j/°C |int.j/° C| Percent int. j/° C |int.j/° C| Percent 
13498. 6 —1.0 —0. 007 13501. 2 +1.6 . 

—9.2 —. 068 13510. 2 +10.6 
+0.5 +. 004 13499. 3 —0.3 
—4.1 —. 030 
—0.1 —. 001 13499. 6 +3. 2 
+1.8 +. 013 
































V. MATERIAL 


The isoprene used was prepared by Bekkedahl, Wood, and Wojcie- 
chowski [3] by cracking dipentene vapor and was purified by frac- 
tional distillation. Their ebulliometric measurements on the freshly 
distilled isoprene indicated that the material was of a very high degree 
of purity. However, the measurements of CO, formed in combustion 
in the present investigation showed that the samples absorbed signifi- 
cant amounts of impurity upon exposure to the air during the filling 
of the glass bulbs used to contain the combustion samples. The iso- 
prene was kept in glass bottles, closed in some cases by ground-glass 
stoppers, in others by cork stoppers, and it is possible that there was 
also some contamination of the isoprene as a result of leakage of air 
past the stoppers. Correlation of the results of the CO, determinations 
vith the heat of combustion measurements indicated that there was 
little or no heat of mixing or combination of the impurity with the 
isoprene, the effect of the impurity being merely to reduce the heat 
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of combustion and the amount of CO, formed in combustion in the 
same proportion. For this reason values of heat of combustion of the 
isoprene obtained by using the mass of CO, formed in combustion ag q 
measure of the mass of isoprene burned, should be fairly reliable jn 
spite of the relatively large amounts of impurity in the samples, 

The oxygen used was ordinary commercial oxygen from which com. 
bustible impurities were removed by passing the gas over copper 
oxide maintained at a temperature of 700 to 800°C. The oxygen was 
then passed through Ascarite before being admitted to the bomb ip 
order to remove any CO, which might be present. 


VI. RESULTS AND CONCLUSIONS 


The results of heat of combustion measurements on several samples 
of isoprene are shown in table 2. Columns 1, 2, and 3 of this table 
give the sample number, experiment number, and mass (weight in 
vacuo) of isoprene burned in the experiment. Column 4 gives the 


TABLE 2.—Results of measurements on isoprene 





Ap- 
Ob- 

proxi- 
served mate 


COs tousd__|heatof| Qa 
COs eale "| com- r — 


tion, Qz “aon 


Remarks 








int. j/g 
46, 293 | Freshly distilled, 


Same sample as A; but had 
been kept in a bottle at 
room temperature 17 days. 


Had been kept in a bottle for 
unknown length of time. 


} Freshly distilled. 


Same lot of material as C; but 
had been kept in a bottle in 
icebox for 7 days. 


Freshly distilled inte anti- 
oxidant. 


Habel distilled, bulbs filled 
immediately. 


Same lot of material as Ei, 
_ filled about 3 hours 
ter. 


Same lot of material as Ei, 
13 bulbs filled about 26 hours 
later. 
14 |0. ‘ 253 46, 242 | Same lot of material as i, 
bulb filled 4 days later. 


15 1. 

16 /1. Freshly distilled, bulbs filled 
17 |1. 45, 996 immediately. 

18 j1. i 46,019 





























* Samples Es and E, were kept outdoors ina bottle at temperatures ranging from —3 to +11° C for 1 and 
4 days, respectively. i 
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ratio, r, of the mass of carbon dioxide formed in the complete com- 
bustion of the sample to the amount calculated from the mass of 
sample, on uy Pape the sample was pure isoprene. The atomic 
weights used in this calculation are C=12.01, H=1.0078, O= 16.000. 
Column 5 of table 2 gives the observed heat of combustion per gram of 
isoprene (Q,) under the conditions of the bomb process, when the 
amount of isoprene burned is taken as the mass determined by weighing 
the sample of isoprene. Column 6 oe values of Q,/r, that is, the 
heat of combustion of isoprene under the conditions of the bomb 
process when the mass of isoprene burned in each experiment is cal- 
culated from the observed mass of carbon dioxide formed in combus- 
tion of the isoprene. The values of Q,/r are therefore the values of Q, 
corrected on the assumption that the departures from unity of the 
corresponding values of r are due to inert impurity in the isoprene. 
Column 7 gives Washburn’s approximate correction [13], and column 
8 the value of —AU,, that is, the energy decrease per gram of pure 
isoprene for the reaction 


C;H; (liq isoprene)+7 O, (gas)=5 CO, (gas)+4H,0 (liq) 


at 30° C and under a total pressure of 1 atmosphere, assuming again 
that the departures from unity of the ratio r are due to inert impurity. 

Table 2 shows that the values of r are considerably lower than unity. 
There is also evidence that r decreased with time when the isoprene 
was kept in stoppered bottles, as previously mentioned. For example, 
for sample A the value of r dropped from about 0.998 to 0.983 when 
the sample was kept at room temperature for 17 days, and for sample 
C the value of r dropped from about 0.994 to 0.992 when the sample 
was kept in an icebox for 7 days. On the other hand, the value of r 
for sample Z did not show any appreciable decrease when the sample 
was exposed to temperatures ranging from —3 to +11° C for 4 days. 

It was thought that the departures from unity of the values of r 
might be caused by oxidation of the isoprene. To test this hypothesis 
a sample was distilled into a vessel containing an antioxidant, and one 
determination was made of the heat of combustion of the sample 
(experiment 8). Unfortunately, the attempt to measure the mass 
of CO, formed in combustion in this experiment was unsuccessful due 
to loss of part of the CO,. However, the value obtained for the 
heat of combustion per gram of sample (Q,) was so low as to indicate 
that this sample of isoprene contained more rather than less impurity 
than most of the other samples. It was concluded, therefore, that 
the impurity was not oxygen or an oxidation product but was probably 
water absorbed from the atmosphere. This conclusion is supported 
by unpublished data obtained by Bekkedahl and Wood in combustion 
analyses of sample H. The isoprene they used was contained in glass 
bulbs which were filled at the same time as some of those used in the 
calorimetric combustion experiments. The value for the ratio r 
obtained in their first experiment was 0.9926, which is in good agree- 
ment with the values of r for sample EZ in table 2. Their value for the 
hydrogen content of the sample was 11.83 percent as compared with 
the theoretical value 11.84 —e for pure isoprene. In a second 
combustion analysis part of the carbon dioxide was lost, but the 
observed hydrogen content was 11.79 percent. The fact that the 
carbon content of the isoprene was low and the hydrogen content 
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was nearly that for pure isoprene indicates that the impurity in th 
sample was water, since the hydrogen content of water is nearly thy 
same as that of isoprene. 

It will be seen from table 2 that the values —AU, for most of th F 
experiments are in satisfactory agreement, but that the values fo 
experiments 5, 7, and 12, although consistent among themselves, ay, 
higher st about 0.63 percent than those for the remaining exper. 
ments. This is shown more clearly in table 3, where the results of 
experiments 5, 7, and 12 are listed separately from those of the other 
experiments. 


TABLE 3.—Comparison of the results of experiments 5, 7, and 12 with those of othe 
experiments 





a 


Experiment —AUR Deviation from mean 








int. j/g 
46, 























The difference between the mean of experiments 5, 7, and 12 and 
the mean of the other group of experiments in table 3 was found to 
be about 85 percent of the energy of evaporation of isoprene, and this 
suggested that the glass bulbs used to contain the liquid isoprene had 
broken in experiments 5, 7, and 12 and allowed liquid to evaporate to 
fill the bomb with saturated isoprene vapor before the beginning of 
the calorimetric experiments. Unfortunately, an accurate calcula- 
tion of the effect of evaporation of isoprene in the bomb is not possi- 
ble, because no data are available on either (a) the effect on the 
partial pressure of the isoprene of the increase in concentration of 
water in the liquid phase resulting from evaporation of part of the 
isoprene, and possibly also from absorption of water from the water- 
saturated atmosphere in the bomb; or on (5) the effect of the pressure 
of the oxygen in the bomb on the density of the saturated isoprene 
vapor. However, it is possible to make fairly reasonable assump- 
tions regarding the possible effects of factors (a) and (6), and to cal- 
culate the thermal effect of evaporation of isoprene in the bomb on 
the basis of these assumptions. In table 4 are given the results of 
recalculating the data of experiments 5, 7, and 12 on the assumption 
that the bomb was filled with isoprene vapor. ‘The amount of 1s0- 
prene in the vapor phase was calculated in the following three ways: 
(1) The effect of the factors (a) and (6) mentioned above was neg 
lected. The vapor pressure of pure isoprene at the initial temper 
ture of the calorimeter in each experiment was calculated from 4 
vapor-pressure equation of the form log p=a—b/T, where the cot- 
stants were determined from values of boiling point and change of 
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boiling point with _— reported by Bekkedahl, Wood, and 
Wojciechowski [3]. The actual pressure of the isoprene vapor in the 
bomb was assumed to be lower than the pod pressure of pure 
isoprene by an amount calculated by means of Henry’s law from the 
amount of impurity (assumed to be water) indicated by the observed 
values of r (table 2). The specific volume of the isoprene vapor was 
calculated by means of the relation PV/7'=constant, where the con- 
stant was determined from the values of V and 7 at P=1 atmos- 
phere [3]. (2) The liquid isoprene was assumed to have absorbed an 
amount of water sufficient to reduce the density of isoprene in the 


vapor phase to 80 percent of the density of vapor in equilibrium with 


pe liquid isoprene. The effect of the presence of oxygen in the 


mb on the density of the isoprene vapor was neglected. (3) The 
presence of the oxygen in the bomb was assumed to have caused com- 


| plete evaporation of the isoprene. 


Taste 4.—Results of three methods of recalculating the daia of experimenis 6, 7, 


and 12 





Recaleculated values of ~AUz at 30° O 





Experiment 
Method 1 Method 2 Method 3 





int. j/g int. Ve int. j/g 
46, 278 rH 332 46, 227 


46, 354 391 46, 214 
46, 208 46, 268 46, 189 


46, 280 46, 330 46, 210 

















It will be seen from table 4 that the mean of the results recalculated 
by method 1 are in agreement within 0.03 percent with the mean of 
the first group of values given in table 3. The mean results obtained 
by methods 2 and 3 differ from the mean of the first group of values in 
table 3 by +0.14 and —0.12 percent, respectively. As the effects of 
the two factors neglected in method 1 are relatively small and of 
opposite sign they will partially cancel, and it may be concluded that 
the result obtained by method 1 is not seriously in error. The agree- 
ment of this result with the mean of the first group of values in table 3 
may therefore be considered as confirmation of the hypothesis that 
the high values of —AU, obtained in experiments 5, 7, and 12 were 
due to breaking of the glass bulbs containing the isoprene and the 
consequent evaporation of part of the isoprene. 

In view of the above discussion, it appears justifiable to discard the 
results of experiments 5, 7, and 12, and to adopt as the value for the 
decrease in energy for the reaction 


CsHs (liq)+7 O, (gas)=5 CO, (gas)+4 H,0 (liq) (1) 


at 30° C and under a pressure of 1 atmosphere, the value of —AU, 
obtained from the first group of experiments of table 3, namely, 
46,265 international joules per gram (weight in vacuo), or 3151.2 
international kilojoules per mole of C;Hs. The corresponding value 
of —AH for the above reaction is obtained by adding the value of 
~—A(pv) for the reaction, 5.1 international kilojoules per mole, to the 
value of —AUp, giving for —AH the value 3156.3 international 
ojoules per mole. 
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The uncertainty of the mean of the results given in the first group 


sa—CO« eo 
in table 3, defined as 24/ Zep is 0.035 percent. Combining this 


with the uncertainty of the energy equivalent of the calorimeter, 0,09 
percent, and the estimated uncertainty, 0.023 percent, in the valye 
used for the heat of combustion of benzoic acid, there is obtained for 
the uncertainty of the value given for the heat of combustion of liquid 
isoprene 





+ (0.023)?+ (0.026)?+ (0.035)?=0.049 percent. 


Using for the specific heats at constant pressure of liquid isoprene 
gaseous oxygen, gaseous CO, and liquid water the values 153.3 {2} 
29.18 [5], 36.62 [5], and 75.18 [6] international joules per mole, respec. 
tively, there is obtained the value —0.126 kilojoule per mole degre 
centigrade for the temperature coefficient (=—d/dT (—AH)) of the 
heat of combustion of isoprene at constant pressure. Hence the value 
derived from this investigation for the heat of combustion of liquid 
isoprene at 25° C and a constant pressure of 1 atmosphere to form 
liquid water and gaseous CO, is 3156.9+ 1.6 international kilojoules 
per mole. Adding to this the value for the heat of vaporization of 
isoprene calculated by Bekkedahl, Wood, and Wojciechowski {3}, 
25.9 +0.3 kilojoules per mole, there is obtained the value 3182.8+1,7 
international kilojoules per mole for the heat of combustion of gaseous 
isoprene at 25° . and at a constant pressure of 1 atmosphere, the 


products of combustion being liquid water and gaseous CQ. : 
The data on the heat of combustion of isoprene can be combined [ 


with data reported previously [8] for the heat of combustion of rubber 
to calculate the heats evolved in the reactions. 


xC;H, (liq isoprene)=(C;H;).(rubber hydrocarbon) (2) 
zC,;H, (gaseous isoprene)=(C,;H;),(rubber hydrocarbon) (8) 


Several values were reported in the reference cited for the heats of 
combustion of different samples of rubber. The value reported for 
the heat of combustion of “steam purified” rubber, which was the 
purest of the samples on which measurements were made, was 45239 
+90 international joules per gram (weight in vacuo), or (30813 
+6.0)z jalacemtinial kilojoules per mole at 30° C. Using the value 
128.7z joules per mole degree centigrade reported by Bekkedahl and 
Matheson [1] for the specific heat of rubber hydrocarbon at 27.5° (, 
and the values given previously in this paper for the specific heats of 
gaseous oxygen, gaseous carbon dioxide, and liquid water, the tem- 

erature coefficient of the heat of combustion of rubber hydrocarbon 
is calculated to be —0.1512 international kilojoule per mole. The 
heat of combustion of rubber hydrocarbon at 25° C is therefore cal- 
culated to be (3082.1 +6.0)z international kilojoules per mole. Com- 
bining this value with the values given previously for the heats of 
combustion of liquid and gaseous isoprene there are obtained the 
values (74.8 +6.2)z and (100.7 +6.3)z for the heats evolved in the 
reactions (2) and (3), respectively, at 25° C and a constant pressure 


of 1 atmosphere. 





Heat of Combustion of Isoprene 


VII. REFERENCES 


N. Bekkedahl and H. Matheson, J. Research NBS 15, 503 (1935) RP844. 

N. Bekkedahl and L. A. Wood, J. Research NBS 19, 551 (1937) RP1044. 

N. Bekkedahl, L. A. Wood, and M. Wojciechowski, J. Research NBS 17, 
883 (1936) RP951. 

H. C. Dickinson, Bul. BS 11, 189 (1914) 8230. 

Int. Critical Tables 5, 80 (1929). 

Int. Critical Tables 5, 113 (1929). 

R. 8. Jessup, J. Research NBS 18, 115 (1937) RP966. 

R. 8. Jessup and A. D. Cummings, J. Research NBS 13, 357 (1934) RP713. 

R. 8. Jessup and C. B. Green, J. Research NBS 13, 469 (1934) RP721. 

8. B. Lebedev, M. A. Chochlovkin, and A. V. Kalacheva, J. Applied Chem. 
(USSR) 9, 1605 (1936). 

(11] T. W. Richards and F. Barry, J. Am. Chem. Soc. 37, 993 (1915). 

{12} F. D. Rossini, BS J. Research 6, 37 (1931) RP260. 

(13) E. W. Washburn, BS J. Research 10, 525 (1933) RP546. 


COO o > i 





= 
S 





WasHINGTON, February 4, 1938. 























U. §. DepARTMENT OF COMMERCE NATIONAL Burgau oF STANDARDS 


RESEARCH PAPER RP1094 


Part of Journal of Research of the National Bureau of Standards, Volume 20, 
May 1938 





ELECTROMOTIVE FORCE OF SATURATED WESTON 
STANDARD CELLS CONTAINING DEUTERIUM OXIDE 


By Langhorne H. Brickwedde and George W. Vinal 





ABSTRACT 


Saturated Weston standard cells in which deuterium oxide is substituted for 
some or all of the normal water may be expected to have electromotive forces 
slightly different from those made with normal water. A study has been made of 
the characteristics of 12 cells, prepared with identical materials, except for the 
isotopic composition of the water. The initial samples of water contained 0.018 
(normal water), 1.036, 3.121, and 6.848 mole percent of D,0. The ratios of moles 
of D,O to H;O in the finished cells were less than in the original samples, because 
of dilution by the water of crystallization of the cadmium sulphate. The con- 
stancy of the electromotive force of each of the 12 cells is shown by the small 
average deviation of single observations, extending over several months, from the 
respective mean values. In no case did it exceed 0.6 microvolt. The decrease 
in electromotive force as a result of using the D,O amounted to 2.9 microvolts per 
mole percent. No significant difference was found in the temperature coefficient 
for the various D,O concentrations. Hysteresis measurements indicate a slight 
superiority for the cells containing added D,O. Measurements of the internal 
resistance of the cells at 20 and 28° C indicate a ratio of Ro/R.:3=1.24. The cells 
containing D,O are fully as reproducible and constant as those made with normal 
water, and in some respects they may be superior. 
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I. INTRODUCTION 


Saturated Weston standard cells in which deuterium oxide, D,O, is 
substituted for some or all of the H,O, may be expected to have an 
electromotive force slightly different from those made with normal 
water. The free-energy change, which determines the electromotive 
force of the cell, is the change accompanying the reaction: 


Cd(10% amalgam) + Hg,SO,(s)+-8/3H,O(in satd. CdSO, soln.) 
+52Hg(1)+CdSO,.8/3H,0(s) 


Since water enters into a part of the cell reaction involving a free- 
energy change, the electromotive force developed should depend on 
the isotopic composition of the water present. 

The three samples of deuterium oxide used in this work were pre- 
pared and standardized by the physical chemistry section of this 
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Bureau. These samples _were of sufficient amount to construc 
several cells of the same size as those which are used customarily fo; 
primary standards of electromotive force. 


II. EXPERIMENTAL PROCEDURE 


A study has been made of the characteristics of a group of 19 
saturated cells constructed on July 27, 1937. These cells were pre- 

ared in the same way and with identical materials except for the 
isotopic content of the water used in making the solutions of cadmium 
sulphate. The materials were purified according to methods de. 
scribed in an earlier publication,' and the cells were made in the usual 
manner. Four samples of water containing, respectively, 0.018 
(normal water), 1.036, 3.121, and 6.848 mole percent of D,O ? wer 
saturated with small CdSQ,.8/3H,O crystals. These crystals had 
been moistened with a small amount of H,SO,, so that the saturated 


CdSO, solution of normal water for the control cells was 0.0226 V : 
in acid. Since the crystals were well mixed and since the same ratio | 


of crystals to water was used in all cases, each solution should be of 
the same acidity. In accordance with this Bureau’s practice of 
referring to the acidity of a cell as that of the acid before saturati 

with neutral CdSO,.8/3H,O, these cells are designated as 0.03 N acid, 


The ratios of moles of D,O to moles of water in the finished cells F 


were less than in the samples initially. The dilution by the wate 


of crystallization of the cadmium sulphate used in saturating the [ 


solution and added as excess solid to the cells is not definitely known 
because of lack of knowledge of the rate and extent of exchange 


between the H,O of the solid CdSO,.8/3H,O crystals and the D,O of 
the liquid. 7 


TABLE 1.—Dilution of deuterium oxide by water of crystallization of the cadmium 
sulphate 





Assuming complete 


AE differ- Assuming no exchange exchange 


Initial D3:0 
concentra- | ®2¢e —_ 
tion consro 
D:0 con- |AE permole} D20 con- 
centration centration 








Mole % Mole % Mole % 
0. 018 0. 018 0. 018 
1. 036 . 85 
3. 121 


























* Control cells, containing normal water. 


If we assume no exchange between the liquid and crystals and 
attribute the dilution of the D,O to the HO of the dissolved crystals 
alone, the D,O mole percentages of the water in the solutions would 

1G, W. Vinal, D. N. Craig, and L. H. Brickwedde. Trans, Electrochem. Soc. 58, 146 (1935). 


2 These compositions were calculated from the densities of the samples, using the equation Nop 


=9.2385AS-+0.2597 (AS)*. In this equation, which was obtained from the data of L. G. Longsw — 


Chem. Soc. 59, 1483 (1937), by least-square computation, No.0 is the mole fraction of deuterium 


a 
AS is @—1~ 5.90705. 1, 
in which 0,99705s is the density, ao of D:0-free water. 





—> mom -f- 4 ©. © 4 


< 


Brickwedde) Emf of Cells Containing Deuterium Oxide 601 


be 0.018, 0.85, 2.57, and 5.64. If, on the other hand, we assume that 
the H,O of the solid entered into an exchange with the liquid and came 
to an equilibrium such that the D,O/H,0O ratio in the crystals is the 
same as in the solution, the mole percentages of D,O in the solutions 
would be 0.018, 0.76, 2.27, and 4.98. The hydrogen added by the 
H,SO, may be neglected since, in the most concentrated D,O sample, 
it reduced the concentration of D,O by only 0.002 mole percent. 

In table 1 the mole percentages of the original samples are compared 
with the percentages in the finished cells on the assumptions of no 
exchange and complete exchange in the crystals. The difference in 
results based on these assumptions is of little importance since the 
variations in AE per mole percent are within the limits of scattering 
of the values in the various groups. From the standpoint of the 
present work it is immaterial, therefore, whether the exchange is 
complete, except as the constancy of electromotive force over a long 
period of time may be involved. In this paper we have assumed that 
the second of the above alternatives is correct and the mole percentages 
' of D,O in tables 2 to 5 are expressed on this basis. In continuing the 
work, however, according to ges for constructing cells containing 
much higher percentages of D,O, this matter will become of increasing 
' importance and a more definite answer, based on experimental results, 
| is to be expected. i 

Another uncertainty in estimating the actual concentration of D,O 
in the finished cells arises from the possible partial dehydration of the 
cadmium sulphate crystals on exposure to the air prior to using them 
for preparing the saturated solutions. Partial dehydration of the 
crystals would result in less dilution of the D,O than is shown in 
table 1, but the uncertainty in the concentration would hardly be 
_ great enough to affect AE per mole percent by an amount large 
enough to be significant. 


III. EXPERIMENTAL RESULTS 


Repeated measurements of the electromotive force at 20° C for each 

cell are given in table 2. The constancy of the electromotive force of 
each of the 12 cells is shown by the small average deviation of single 
observations given in the last column. In no case is the average 
deviation more than 0.6 microvolt. The scattering of the mean 
values of individual cells in their respective subgroups is not more 
than 1.4 microvolts. The constancy and reproducibility of the cells, 
as indicated by the close agreement of the electromotive force values 
for each isotopic concentration, are important indications of the 
validity of the effect of the D,O. 
_ The mean electromotive force at 20° C of each subgroup is plotted 
in figure 1 with respect to the D,O mole percent of the water in the 
saturated solution. The relation of electromotive force to the D,O 
content is evidently linear within the range of the present work. The 
slope of the line indicates a decrease of electromotive force amounting 
to 2.9 microvolts per mole percent.’ 

The electromotive forces at 20 and 28° C are compared in table 3. 
No significant difference was found in the temperature coefficients for 

'Tn contrast with our results, Yoneda reported to the Comité Consultatif d’Electricité, 1937, a decrease 


of 27 microvolts for 1/50 percent of DyO. His results on three cells will probably be published in volume 18 
of the Procas Verbaux, Somite International des Poids et Mesures. . — 
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the various D,O concentrations. The mean change in electromotiy 
force going from 28 to 20° C, 374.0 microvolts, lies between the valy f 
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Fraure 1.—Relation of electromotive force to mole percent of D,O 


380.5 calculated from the international temperature formula! and > 
370.1 calculated from the formula of Vigoureux and Watts.’ 


‘ Report of London Conference on Electrical Units and Standards, schedule C (1908). 
5 Proc. Phys. Soc. (London) 45, 172 (1933). 
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TaBLE 3.—Change in electromotive force with temperature 


[Vol. 9 





Cell No, 


Mole per- 
cent of D:0 


Emf at 28.00° 
C (mean 
values) 


Emf at 20.00° 
C (mean 
values) 


Difference 


Mean differ. 





0. 018 


Vv 
1, 017899; 


v 
1. 018272 


BY 
373.7 


ence 


i, 
BV 
8744 


374.4 
374.9 


fA 
374s ~. 


. 018 975 72 
. 018 976 71s 


- 76 951 689 
. 76 95g 703 


- 76 965 700 373.5 


92s 663 373.5 
91s 665 374.7 
925 669 874.4 
373.7 


373.6 
373.9 


3743 


844 58; 
83s 574 
843 582 


S22 SSN 


3734 

















mH fo pope 





The results of measurements of the hysteresis of these cells ar 
shown in table 4. In the first test, the cells were moved from a bath 
at 20° C to another at 28° C. The electromotive forces were meas. 
ured after 2 hours in the 28° bath and compared with the mean of the 
equilibrium values subsequently observed at that temperature. In 
the second test, the cells were transferred back to the 20° bath and 
again measurements made after 2 hours at the new temperature were 
compared with the mean of the equilibrium values at 20° C. The 
hysteresis is somewhat less for cells containing added D,O than for 
the cells made with normal water. 

The internal resistances of the cells at 20 and 28° C and the ratio of 
the values at the two temperatures are shown in table 5. The internal 
resistances, R, were calculated from the relation, Pe a X10, 
where E is the electromotive force and E’ the terminal voltage ob- 
served when the cell was discharging through an external resistance 
of 10’ ohms. 


TABLE 4.—Hysteresis measurements 20 to 28 to 20° C 


[The cells were measured 2 hours after abruptly changing the temperature 8° C.] 





20 to 28° C 28 to 20° C 


Mole per- 
cent of 
D:0 





Emf 2 hr 
after 
change 


Equilib- 
rium emf 
at 20° 


Equilib- 
rium emf 
at 28° 





v Vv 
1. 017899; 1. 018268; 
897 


















































Emf of Cells Containing Deuterium Oxide 


TaBLE 5.—IJnternal resistance of the cells 





Resistance 
7 Mole per- Ratio Mean ratio 

Cell No. cent D30 Ru/R Ru/Ras 

20.00° C 28.00° C 








Ohms Ohms 
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IV. CONCLUSIONS 


Aside from the theoretical interest in the effect of deuterium oxide 
on the electromotive force of the cells, there are practical reasons for 
continuing the investigation. The results given in this paper indicate 
that the cells containing D,O are fully as reproducible and constant as 
those made with normal water. With respect to hysteresis they may 
be superior. Observations over a longer period of time are needed to 
determine whether the cells are less subject to the troublesome forma- 
tion of gas over the amalgam terminal. We believe the results so far 
obtained are sufficiently promising to justify extending the work to 
cover a wider range of concentrations of the deuterium oxide. 


The authors thank Dr. E. R. Smith, who furnished the deuterium 
oxide, for helpful suggestions and in particular for standardizing the 
samples. 


Washington, December 23, 1937. 


[AutHoRr’s Norr.—Because of the time which has elapsed since the manuscript 
of this paper was completed, it is desirable to state that subsequent measurements 
on these cells show them to maintain substantially the same relative differences. 
Slight increases in electromotive force have been observed in all cells, but the 
average increase for those made with ordinary water, 2.; microvolts, exceeds that 
of the other cells, 1., microvolts, which argues against the likelihood of lon 
continued exchange of D,O between the crystals and solution. March 25, 1938. 
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ANALYSIS OF PHOSPHATE ROCK 
By James I. Hoffman and G. E. F. Lundell 







ABSTRACT 


The methods used at the National Bureau of Standards in the analysis of the 
two Standard Samples of Phosphate Rock, 56a and 120, are described, and a 
summary of the results obtained by nine cooperating analysts is given. The 
' usual chemical analysis is “tg, ee mig bya —" spectrochemical examina- 
' tion, and one constituent, boric oxide, has mn quantitatively determined by 
» chemical and spectrochemical methods. 
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I. INTRODUCTION 


In analyzing the Standard Samples of Phosphate Rocks, 56a anj 
120, issued by the National Bureau of Standards, a number of jp. 
provements in old methods were made, and some new ones were de. 
veloped. Most of the methods that were used when the first Standan 
Sample of Phosphate Rock was issued are still in use, but many of 
them have been modified to such an extent that the publication of the 
new and revised methods seems desirable. This paper brings the 
previous methods! up-to-date, but the discussion given in the ab 
paper is not repeated. Certain modifications and suggestions by the 
analysts who cooperated in this work are incorporated in the diseys. 
sion, but only the methods used at this Bureau are given in section 
III, Methods of Analysis. 


II. GENERAL DISCUSSION 
1. MOISTURE 


The moisture content of air-dried samples of finely ground phos. 
phate rock may vary considerably with changes in atmospheric condi- 
tions. All of the present work was therefore done on samples dried 
for 1 hour at 105° C, so that the results reported by all the analysts 
might have a common basis. Tests at this Bureau indicated that 
variation of +5° C in the temperature is permissible in the case of 
these particular rocks, and that even if a sample is exposed for 40 
hours in an atmosphere of carbon dioxide saturated with water and 
then dried, its weight is practically the same as it would have been if 
it had been dried without this severe exposure.” 


2. PHOSPHORIC ANYHDRIDE 


Most of the cooperating analysts used either the alkalimetric or the 
gravimetric method recommended by the Association of Official Agn- 
cultural Chemists.* The alkalimetric method involves precipitating 
the phosphate as ammonium phosphomolybdate in acid solution, dis 
solving the precipitate in an excess of a standard solution of sodium 
hydroxide, and then titrating the excess with standard nitric acid. Ih 
the gravimetric method the ammonium phosphomolybdate is dis 
solved in dilute ammonium hydroxide, and the phosphate is then pre- 
cipitated as magnesium ammonium phosphate and ignited to the 
pyrophosphate, Me.P.0;. The volumetric method is used almost 
exclusively in control laboratories and yields fairly accurate results 
in a very short time if conditions are closely controlled, especially if 
the sodium hydroxide is standardized by means of a standard sample 
of phosphate rock. The gravimetric method, involving only one 
precipitation with magnesia mixture, yields results that are slightly 
high, usually about 0.2 percent. pe 

n the standardization of Samples 56a and 120, only gravimetne 

methods involving two precipitations with magnesia mixture were 

nao cr AB aL On AO, ae ses 
Tenn., states that more uniform results are obtained if their superphosphates and ‘‘metap! 


dried at room temperature for 24 to 48 hours over sulfuric acid in a partially evacuated desiccator, than i 


they are dried in an oven at 105° C. 18-2, 
3 Official and Tentative Methods of Ana of the Association of Official Agricultural Chemists, ps 


4th ed. (1935), published at Washington, D. 
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used at this Bureau, because previous work ‘ showed that such a 
rocedure is necessary if accurate results are desired. Two such 
methods were employed, One involves a precipitation of the phos- 
phate as ammonium phosphomolybdate, followed by two precipita- 
tions with magnesia mixture, and is described in detail elsewhere.’ 
The other is the recently published method ° in which two precipita- 
tions with strong magnesia mixture in the presence of ammonium 
citrate serve to separate the phosphate from all the other constituents 
of the phosphate rock without the preliminary precipitation with 
molybdate reagent. 

For comparison, the results obtained at this Bureau by the official 
_ gravimetric method of the Association of Official Agricultural Chemists 
are included in the Summary (tables 2 and 3). 


3. FERRIC OXIDE AND ALUMINA 


Total ferric oxide and alumina are seldom determined. The manu- 
' facturers of fertilizers are usually interested only in “soluble iron and 
alumina.” To keep the results on a common basis, the cooperating 
analysts have agreed that the percentages of soluble Fe,O; and Al,O; 
shall represent the amounts dissolved by transferring a 2.5-g sample 
to a 250-ml flask, adding 50 ml of diluted hydrochloric acid (1+1) 
containing 1 g of boric acid, covering the neck of the flask with a 
watch glass, and boiling for 30 minutes. Tests have shown that the 
boric acid has practically no effect on the results, but that attack on 
the glassware by the fluorine in the rock is somewhat lessened. 
| By analysis of synthetic solutions containing known amounts of 
alumina, the method for soluble alumina described in the succeeding 
pages has been found to yield accurate results. Some of the cooperat- 
ing analysts obtained satisfactory results by shorter methods, usually 
| of an empirical nature. Such methods are not objectionable, espe- 
cially if the general type of the rock to be nalveed is known. The 
preparation of the solution is probably the most important part of 
this determination, because it 1s essential that all the analysts start 
with solutions containing the same amount of dissolved material. 

Considerable quantities of alumina are retained by the precipitate 
obtained with sodium hydroxide if this reagent is used for separating 
iron and other elements from aluminum in materials containing appre- 
clable quantities of magnesium.’ The error caused by this retention 
may become serious in the analysis of some of the newer fertilizers of 
high magnesia content. 

In the determination of total ferric oxide and alumina a mixture of 
perchloric, hydrochloric, and nitric acids is recommended for decom- 
posing the material, because this mixture effects a rapid and fairly 
complete decomposition of phosphate rocks, and at the same time 
ates fluorine, which might cause errors in the determination of 
iron.* Nitric acid should always be used with perchloric acid in this 
attack because, without its use, serious explosions may occur as a 
result of the rapid interaction of perchloric acid and the organic matter 
contained in most phosphate rocks. 

‘BS J. Research 5, 279 (1980) R 


P200. 
(4 Aus, Oficial Agr. Chem. 8, 104 ipa). i co a a 
NBS IN 30 sgn of Phosphoric aniydride im phosphate rock, superphosphate @ phosphate.” J. Researc 


'G. E. F. Lundell, J. 1. — and H. A. Bright, Chemical Analysis of Iron and Steel, p. 545 (J. Wiley 
iG. E. F. Lundell, Applied Inorganic Analysis, p. 775-776 (J. Wiley & Sons, New 
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4. CALCIUM OXIDE 


Total calcium oxide was determined by precipitating the caleiyn 
as oxalate in a slightly acid (pH, 3 to 4) solution. Bromphenol bly, 
was used as the indicator in adjusting the acidity, and two precipi. 
tations were made to free the calcium oxalate from phosphate, iron, 
and aluminum. 

5. MAGNESIA 


Most of the analysts determined magnesia by the method ® published | 


at this Bureau in 1932 or by some slight modification of it. In this 
method calcium is removed by precipitating it as sulfate in alcoholic 
solution, and the magnesia is then determined in the filtrate (after 
evaporating the alcohol) by precipitating it with phosphate in the 
presence of ammonium citrate, and igniting the magnesium ammo- 
nium phosphate to magnesium pyrophosphate, Mg.P,O,. 


6. FLUORINE 


If a phosphate rock is fused with sodium carbonate and the cooled 
melt is leached with water, much fluorine remains in the insoluble 
residue.” For example, by such a treatment of 0.5-g samples of 
Phosphate Rocks 56a and 120, 1.4 and 1.6 percent of fluorine out of 
total of 3.5 and 3.8 percent, respectively, were retained by the insoluble 
residue. The work of Reynolds and Jacob shows that three fusions 
and extractions usually do not yield more than 90 percent of the 


fluorine. As it is obviously impracticable to separate it by sucha F 
procedure, analysts who determine fluorine in mes singe, materials | 


almost always resort to a method involving distillation. 


All the analysts who determined fluorine distilled it by the method 


of Willard and Winter ™ or by some modification of it. Perchloric 
acid was used by cooperating analysts 7 and 9 (tables 2 and 3), 
Egbert Janes and H. H. Edwards maintained temperatures between 


135 and 140° C in the distilling flask, by dropping water into the solu- | 


tion during the distillation. A flask fitted with a dropping funnd 
and thermometer similar to that shown in figure 1 was used. The 
analysts of the U. S. Department of Agriculture maintained temper- 
atures between 120 and 150° C by passing steam into the flasks during 
thedistillation.” It was observed in both of the laboratories mentioned 
that free sulfur or pyritic sulfur in the phosphate rock causes high 
results for fluorine if the distillate is titrated with thorium nitrate. 
In the laboratory of the Department of Agriculture this difficulty is 


overcome by redistillation. Janes and Edwards state that this inter | 


ference can be prevented by igniting the sample at a low temperature 

in a silica dish before transferring it to the distilling flask. 
Because of the dangers attending the use of perchloric acid for 

decomposing materials containing organic matter,” es y if 


nitric acid can not be used as a preliminary oxidizer during the heating, F 


it was deemed advisable to develop a method in which perchloric 
is not used. At the National Bureau of Standards the fluorine was 


® The determination of magnesia in oan, BS J. Research 9, 487 (1932) RP 484. 

1” See also, D. 8S. Reynolds and K. D. Jacob, Ind. Eng. Chem. Anal. Ed. 8, 366 (1931). 

1H. H. Willard and O. B. Winter, Ind. Eng. Chem. Anal. Ed. 5, 7 (1933). 

1? D. 8. Reynolds, J. B. Kershaw, and K. D. Jacob, J. . Official Agr. Chem. 19, 156 (1936). 

13 Practically ot poeatets rocks contain organic carbon. See K. D. Jacob, W. L. Hill H. L. a 
and D. 8. Reynolds, Tech. Bul. 364, p. 51 (1933) U. 8. Department of Agriculture, Wash D.C. 
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therefore distilled from a sulfuric acid “ solution and determined in the 
distillate by the lead chlorofluoride method.” At the tempera- 
ture used in this distillation (160 to 170° C), a little sulfuric acid 

asses into the distillate, but this does not interfere in the lead chloro- 

uoride method. bit 5 

Willard and Winter indicated that they obtained complete recovery 
of fluorine from fluorite (CaF,), phosphate rock, or cryolite by direct 
decomposition and distillation with perchloric acid. By the method 
described in the succeeding pages, using the apparatus shown in 
figure 1, complete recovery of fluorine from fluorite was not obtained 
by the use of either sulfuric or perchloric acid. By first fusing the 
sample with sodium carbonate and then distilling, complete recovery 
was obtained. The difference (0.1 to 0.2 percent) between the results 
obtained with fused and unfused samples of phosphate rock is probably 
too small to warrant the introduction of this extra step in the method. 
For routine work, fusion probably is never necessary, but since it is 
not definitely known how the fluorine is combined in all types of phos- 
phate rock, the analyst should make determinations with and without 
preliminary fusions to find out whether direct treatment with acid and 
distillation will yield satisfactory recovery in the type of material 
being analyzed. 

7. SILICA 


Owing to the difficulty of separating all the fluorine from the insol- 
uble residue that is left when the cooled melt of a sodium carbonate 
fusion of a phosphate rock is leached with water, the determination of 
silica is even more difficult in this case than in glasses, enamels, or 


other materials that contain tittle or no a ama In fact, no 


reasonably rapid and accurate method is available. The one here 
described, although somewhat of a makeshift, is the best expedient 
available. The determination of silica therefore remains the most 
difficult and least satisfactory in the entire analysis. 

In the method described, two fusions with sodium carbonate and 
leachings of the cooled melts with water serve to separate most of the 
—_ from the insoluble residue, which still contains some of the 
silica. 

The silica in the alkaline solutions is separated from the fluorine in 
the usual way,'* by first precipitating most of it with zinc nitrate and 
then the remainder with zinc oxide dissolved in ammonium hydrox- 
ide. The silica in the insoluble residue is obtained by adding boric 
acid to lessen the effect of the fluorine which is still retained, and then 
dissolving the residue in hydrochloric acid and dehydrating in the 
usual manner. 

8 MANGANESE OXIDE 


By far the most satisfactory method for the determination of man- 


ganese oxide in phosphate rock is that in which it is determined color- 


imetrically after oxidation to permanganic acid by sodium or potassium 
periodate in strong acid solution. 


4 The use of phosphoric acid has also been suggested for this distillation. 

wy,T, ey —S. A va oe eg of fluorine and of silica in glasses and enamels contain- 
a c % 

"BS J. Research 3, 581 (1929) RP110. 
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III. METHODS OF ANALYSIS 
1. DETERMINATION OF MOISTURE 


Place 5.00 g of sample in a glass-stoppered weighing bottle of su¢) 
diameter that the depth of the sample does not exceed 0.5 cm. 
for 1 hour at 105° C in a well-ventilated oven. At the end of 1 how, 
loosely stopper the bottle, and cool in a desiccator. Momentarily 
raise the stopper, again put it in place, and weigh. The loss in weigh; 
is regarded as moisture. 


2. PHOSPHORIC ANHYDRIDE 
(a) PREPARATION OF THE SOLUTION 


Method A.—Transfer 2.500 g of the dried material to a 300-ml Ei. 


enmeyer flask, add 30 ml of hydrochloric acid and 10 ml of nitric acid, f 


Place a small watch glass over the mouth of the flask, and boil gently 
for 30 minutes. Rinse the glass, add 50 ml of water, heat to boiling, 
and filter through a paper of close texture, catching the filtrate ing 
250-ml volumetric flask. Wash the paper” and the Erlenmeyer flask 
thoroughly with diluted hydrochloric acid * (1+-4), cool the contents 
of the volumetric flask to the proper temperature, and add water until 
the flask is filled to the mark. Mix thoroughly and transfer 50.00 ml 
of the solution (representing 0.5000 g of sample) to a 300-ml Enrlep- 
meyer flask. 

Method B.—Transfer 0.5000 g of the dried material to a 300-ml 
Erlenmeyer flask and add 15 ml of hydrochloric acid and 3 ml of 
nitric acid. Place a small watch glass over the mouth of the flask, 
and boil gently for 30 minutes. Rinse and remove the watch glass, 
and without filtering, proceed with the addition of ammonium citrate, 
hydrochloric acid, and magnesia mixture. 

This method obviates the use of volumetric glassware and requires 
less time than method A. The results obtained are equally accurate 
if the material contains no appreciable amount of compounds of 
phosphorus insoluble in this mixture of acids. 


(b) DETERMINATION OF PHOSPHORIC ANHYDRIDE ® 


To the solution in the 300-ml Erlenmeyer flask, add 30 g of am- 
monium citrate, 10 ml of hydrochloric acid, and 100 ml of strong 
magnesia mixture.” Neutralize the solution with ammonium hydrox- 
ide, using litmus paper as indicator, and add 3 ml of ammonium 
hydroxide in excess. Dilute the solution to 225 to 250 ml with water, 
add a few glass beads, stopper the flask tightly, shake on a shaking 


17 By applying the treatment prescribed, no ag pe is volatilized, and only 0.03 percent of P10; 
remained in the insoluble residue in the case of the National Bureau of Standards —_ 56, and less than 
0.005 percent in Samples 56a and 120. Unless it is known that no peercerms remains insoluble, 

paper and insoluble residue in a, fuse with a small quantity of sodium carbonate, and 

cooled melt with hot water. Filter, acidify the filtrate with hydrochloric acid, and add the 

contents of the volumetric flask. Or, if the filtrate is acidified with nitric acid, the phosp 
precipitated with molybdate reagent, separately determined, and the percentage found 

obtained in the main determination. 

18 This denotes 1 volume of concentrated hydrochloric acid (sp er 1.18) diluted with 4 volumes of water. 
Diluted ammonium hydroxide (1+19) denotes 1 volume of ammonium hydroxide (sp Fi 0.90) diluted with 
19 volumes of water. If no dilution is specified, the concentrated reagent is meant. ‘This system of sped 
fying dilution is used throughout this paper. 

1” This method is described in detail in J. Research NBS 19, 59 (1937) RP1010. 

2% Prepare the strong magnesia mixture as follows: Dissolve 400 g of MgCls.6H30 and 300 g of ammonium 
chloride in 1,500 ml of warm water. When solution is complete, add ammonium hydroxi solu- 
tion is alkaline to litmus. Allow to stand for about 1 hour, filter, and add hydroc 
= the solution is acid to litmus. If properly prepared, the volume of the solution 

iters, 
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machine for about 30 minutes and allow to stand overnight. Filter 
through a popes of close texture, and wash the flask and paper once 
with a small quantity of diluted ammonium hydroxide (1+19). 
Discard the filtrate. li aad 

Dissolve any magnesium ammonium phosphate remaining in the 
flask in about 50 ml of warm diluted hydrochloric acid (1+4) and 
pour the solution through the paper to dissolve the CR OTR that 
was transferred to the filter, and to separate any insoluble residue that 
may have remained if method B was used for the preparation of the 
solution. Wash the paper and flask with more of the same acid, and 
add 0.3 g of citric acid and 1 ml of magnesia mixture to the solution, 
which should have a volume of about 100 ml. While stirring, add 
ammonium hydroxide until the solution is alkaline to litmus, and then 
add 5 ml in excess. Stir occasionally during 30 minutes, and allow 
to stand for 4 hours or overnight. 

Collect the precipitate on a 9-cm paper of close texture, wash with 
cool diluted ammonium hydroxide (1+-19), and transfer the paper and 
precipitate to a weighed platinum or porcelain crucible. Char the 
paper without flaming, burn off the carbon below 900° C, and finally 
ignite to constant weight, preferably in a muffle furnace with pyro- 
metric control, at 1,050 to 1,100° C. 

The weight of the magnesium erent, Mg:P,0,, multiplied 
by 63.79 and divided by the weight of the sample, yields the percentage 
of P,Os. 

Tests have shown that in the analysis of phosphate rock the‘time 
of standing in the first precipitation can be reduced to 3 hours and 
in the second to 2 hours, if the solutions are chilled before shaking or 
stirring and then allowed to stand in ice water or in a refrigerator at 
5 to 10° C. With small quantities (less than 10 mg of P,O;), the 
longer time of standing is essential, and even then it is best to chill 
and shake alternately two or three times. 

An alternative umpire method involving a preliminary precipitation 
of the phosphate with molybdate, followed by two precipitations with 
magnesia mixture is described in a previous publication.” 


3. TOTAL FERRIC OXIDE 
(a) PREPARATION OF THE SOLUTION 


Transfer 2.000 g of the dried sample to a 300-ml Erlenmeyer flask, 
and add 30 ml of diluted nitric acid (1+-2), 5 ml of hydrochloric acid, 
and 10 ml of 60-percent perchloric acid. Place the uncovered flask 
on the hot plate, evaporate until fumes of perchloric acid appear, and 
heat just below its boiling point until nearly all the aci thna been 
volatilized. Do not evaporate quite to dryness. Cool somewhat, 
and drench the residue with 5 int of hydrochloric acid. Add 30 ml 
of warm water, and heat until salts are in solution. Remove the 
insoluble residue by filtering through a paper of loose texture, wash 
thoroughly with hot diluted hydrochloric acid (1+19) and reserve 
the solution. Ignite the paper and residue in platinum, treat with 
hydrofluoric and sulfuric acids to eliminate the silica, and fuse the 
nonvolatile residue with a little sodium pyrosulfate. Dissolve the 
cooled melt in 25 ml of diluted hydrochloric acid (1+19), and pre- 


cipitate the iron in the solution by adding a slight excess of ammonium 
*« The analysis ef phosphate rock. J. Assn. Official Agr. Chem. 8 194 (1924). 
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hydroxide. Collect the precipitate on a small filter, and wash ong 
or twice with water." Dissolve the precipitate in hot dilute 
hydrochloric acid (1+19), and add the solution to the reserve 


solution. 
(b) DETERMINATION OF TOTAL FERRIC OXIDE 


To the combined solutions add a few drops of a saturated solution 
of potassium permanganate, boil to expel chlorine, and then adj 
stannous chloride dropwise until the iron is reduced, as is indicated 
the disappearance of the yellow color. Avoid any large excess. W, 
the inside surface of the beaker with a jet of water, cool rapidly, and 
add at one stroke 10 ml of a saturated solution of mercuric chloride 
Agitate the solution for 1 minute, and then pour it into a 1-lite 
porcelain casserole containing 20 ml of manganese sulfate solution® 
and 450 ml of water. Titrate with a solution of potassium pe 
nate, not stronger than 0.1 N, and correct the result by subtracting 
the volume of permanganate required in a titration of the reagents 
alone. Calculate the percentage of Fe.Q;. 

In most cases it is not necessary to filter and to recover the iron in 
the insoluble residue. The iron can be reduced with stannous chloride gs 
soon as the salts have been dissolved in dilute hydrochloric acid after 
the evaporation with perchloric acid, provided that no brown or black 
particles remain undissolved. Mere visual observation will easily show Ff 
whether it is necessary to recover iron from the insoluble residue. 

The percentages of iron oxide in the insoluble residues from eight 
phosmpate rocks were as follows: NBS Standard Tennessee Phosphate F 

ock 56, 0.05; NBS Standard Tennessee Phosphate Rock 56a, 0.007; 
NBS Standard Florida Phosphate Rock 120, 0.007; Florida pebble 
rock, 0.005; Florida hard rock, 0.005; Tennessee blue rock, 0.01; 
Wyoming phosphate rock, 0.01; and Idaho phosphate rock, 0.01. It 
is evident that the recovery of iron is of consequence in only one of 
the eight insoluble residues. The amount of iron in the insoluble 
residue may depend on the fineness of the sample. 


4. TOTAL ALUMINA 
(a) PREPARATION OF THE SOLUTION 


Transfer 2.500 g of the dried sample to a platinum dish, and add 
30 ml of diluted nitric acid (1+1) and 10 ml of hydrofluoric acid. 
Evaporate to dryness on the steam bath, then add 30 ml of diluted 
nitric acid (1+1), 5 ml of hydrofluoric acid, and 15 ml of 60-percent 
perchloric acid. Evaporate nearly to dryness on the hot plate, cool 
wash the inside surface of the dish with a jet of water, and add 5 ml 
of 60-percent perchloric acid. js evaporate nearly to dryness, 
cool somewhat, and add 30 ml of diluted hydrochloric acid (1+2). 


Heat the solution on the steam bath, and if a residue remains, collect 
it on a small tight filter paper. Ignite the paper in a platinum cruc- 
ble, and treat the residue with a drop of sulfuric acid and about 0.) 
ml of hydrofluoric acid. Evaporate to dryness, ignite for 1 minute 
at 600 to 800° C, and fuse the residue with about 0.5 § of sodium 


carbonate. Dissolve the cooled melt in the original acid solution of 
the sample. 
" The queen psy ammonium by droaide is made to eliminate traces of platinum whieb would 


etermination o " al 
2 Dissolve 200 g of MnS04.4H30 fi Lo0oanl of a0, end add a cooled mixture of 400 ml of 2504, 1,200 
of HO, and 400 ml of approximately 85-percent phosphoric acid. 


o' 
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(b) DETERMINATION OF TOTAL ALUMINA IN THE PRESENCE OF LESS THAN 0.5 
PERCENT OF MAGNESIA 


Nearly neutralize the warm solution (60° C) having a volume of 
| about 200 ml, with a 25-percent solution of sodium hydroxide, and 
pour this solution slowly while stirring into 100 ml of a 25-percent 
solution of sodium hydroxide containing 0.5 g of sodium carbonate. 
Transfer the solution and precipitate to a 500-ml volumetric flask, 
allow to cool, fill the flask to the mark with cool water, and mix thor- 
oughly. Filter through a dry filter into a dry 200-ml volumetric 
flask. Two 200- a may be collected if duplicate determi- 
nations are to be made. 

| Determine the alumina in a 200-ml portion of the filtrate, repre- 
senting 1.000 g of sample, by the phosphate method as follows: 
Acidify the solution with hydrochloric acid, add 10 ml of the con- 
centrated acid in excess, and then 1.0 g of diammonium phosphate. 
Dilute to 350 to 400 ml, add methyl orange indicator, make just 
alkaline with ammonium hydroxide, and restore the pink color with 
several drops of diluted hydrochloric acid (14-5). Heat the solution 
to boiling, and add 15 ml of a 50-percent solution of ammonium ace- 
tate. Continue the boiling for 5 minutes, and then collect the pre- 
cipitate on an 1l-cm filter of close texture (No. 42 Whatman or 
No. 589 S&S blue ribbon), and allow the precipitate to drain. With- 
out washing, transfer the paper and precipitate to the beaker in which 
the precipitation took place, and stir the paper to a pulp. Add 
10 ml of hydrochloric acid and 1.0 g of diammonium phosphate, and 
dilute the solution to 300 ml. Add methyl orange indicator, and pre- 
cipitate the aluminum as before. Filter, and wash with a hot 5-per- 
cent solution of ammonium nitrate until 5 ml of the washings gives 
at most a barely distinguishable opalescence with an acidified solution 
of silver nitrate. Usually, 150 ml of wash solution is sufficient. 

Ignite the paper and precipitate in platinum or porcelain. Heat 
at a low temperature until the carbon is gone, then at about 1,000° C 
until constant weight is obtained. Weigh as AlPO, and multiply by 
0.418 to convert to the weight of Al,O;. A blank should be carried 
through all steps of the method and the weight of the AlPO, corrected 
accordingly. 

In routine analyses two steps of the method may be omitted. 
(1) In decomposing most phosphate rocks the residue that remains 
undissolved after the treatments with nitric, hydrofluoric, and per- 
chloric acids contains only a negligible amount of alumina and no 
calcium oxide. It therefore is usually permissible to ignore the in- 
soluble residue, and, without filtering, to proceed with the sodium 
hydroxide precipitation. (2) Results are only slightly higher if only 
a single phosphate precipitation is made. For routine work it is 
therefore permissible to omit the second precipitation with phosphate. 
Macerated paper must, however, be added, and the washing with a 
solution of ammonium nitrate must be thorough. 

The percentages of alumina in the insoluble residues were deter- 
mined in eight phosphate rocks as follows: NBS Standard Tennessee 
Phosphate Rock 56, 0.03; NBS Standard Tennessee Phosphate Rock 
56a, 0.02; NBS Standard Florida Phosphate Rock 120, 0.02; Florida 
pebble rock, 0.01; Florida hard rock, 0.005; Tennessee blue rock, 


cae Wyoming phosphate rock, 0.005; and Idaho phosphate rock, 
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(c) DETERMINATION OF TOTAL ALUMINA IN MATERIALS CONTAINING MOpp 
THAN 0.5 PERCENT OF MAGNESIA 


To avoid low results for alumina when 0.5 percent or more of mag. 
nesia is present, separate most of the magnesium from the aluminyn 
as follows: To the perchloric acid solution of the sample add 50 ml of 
hydrochloric acid, and dilute to 400 ml. Add methyl orange indicato; 
solution and then ammonium hydroxide until the indicator just tums 
renee Restore the pink color by geen hydrochloric acid dropwise, 

eat the solution nearly to boiling, and add 15 ml of a 50-percent 
solution of ammonium acetate. Boil for 5 minutes, filter, and wash 
the precipitate a few times with a 5-percent solution of ammoniyn 
nitrate. Dissolve the precipitate in 30 ml of hot diluted hydrochloric 
acid (1+2), and proceed with the neutralization and precipitation 
with sodium hydroxide. 

The magnitude of the error caused by magnesium is illustrated by 
the following determinations of magnesia in a ——_ standard pre- 

ared by mixing 2.500 g of National Bureau of Standards Argillaceous 

imestone * 1A with 1.5 g of P,O;. Three determinations by the 
method in which magnesia is not removed, (4 (b), p. 615) indicated 
3.84, 3.90, and 3.91 percent of alumina, respectively, in the limestone, 
whereas 4.14 and 4.17 percent, respectively, were obtained by the 
method in which magnesia is removed before the alumina is separated 
by means of sodium hydroxide (4 (c). 


5. SOLUBLE FERRIC OXIDE AND ALUMINA 
(a) PREPARATION OF THE SOLUTION 
Transfer 2.500 g of sample to a 250-ml flask, add 50 ml of diluted 


hydrochloric acid (1+1) containing 1 g of boric acid, cover the neck 
ob the flask with a watch glass, and boil moderately for 30 minutes. 
Remove the flask from the source of heat, add 20 ml of water, and 
remove the insoluble matter by filtration. Catch the filtrate in a 
250-ml volumetric flask, cool the solution, and dilute to the mark. 
Mix well, and remove suitable aliquot portions for the determinations 
of soluble ferric oxide and alumina. 


(b) DETERMINATION OF SOLUBLE FERRIC OXIDE AND ALUMINA 


Determine soluble ferric oxide in an aliquot portion of the prepared 
solution by the procedure described for total ferric oxide under 3 (b), 
p. 614. 

Determine soluble alumina in an aliquot portion of the prepared 
solution by the procedure described for total alumina under 4 (b) 
p. 615 or 4 (c), p. 616. 


6. TOTAL CALCIUM OXIDE 
(a) PREPARATION OF THE SOLUTION 
Prepare the solution as directed for total alumina under 4 (a) p. 614. 
(b) DETERMINATION OF TOTAL CALCIUM OXIDE 


Transfer an aliquot portion of the prepared solution, representing 
0.5000 g of sample, to a 400-ml beaker, add 50 ml of a saturated solu- 
tion of ammonium oxalate (about 4 g of (NH,),C.0,.H,0 in 100 ml of 


8 This limestone contains 14.11 percent of Sis, 1.63 percent of FeyOs, 4.16 percent of AlsOz, 0.16 percent 
of TiOs, 41.32 percent of CaO, and 2.19 percent of MgO. 





— 2 See 


ee Co SY or 


oa aS 3. 


= 


— oO 


Hofman Analysis of Phosphate’ Rock 617 
water) and dilute the solution to 200 ml. Add 2 drops of brom 
phenol blue indicator solution,™ and heat to 80 to 90° C, and then 
add ammonium hydroxide until the color in the hot solution changes 
from yellow to a light green but not to a distinct blue. This yields a 
solution having a pH value between 3 and 4. Digest on the steam 
bath for 1 to 2 hours, cool, filter through a paper of close texture, and 
wash the precipitate in the beaker and on the paper two or three times 
with ammonium oxalate-oxalic acid wash solution (2 ¢ of 
(NH,),C:04.H,O and 1 g of H,C,0,.2H,O dissolved in 1 liter of water). 
Actual determinations of the pH of the solutions by means of a Cole- 
man Model 3D pH Electrometer showed that satisfactory results can 
be obtained at pH values between 3 and 4.4, but that low results are 
obtained if the acidity is increased to a pH of 2.0 to 2.4. 

Transfer the paper and the precipitate to the beaker in which the 
calcium was precipitated, stir the paper to a pulp, and add 20 ml of 
diluted hydrochloric acid (1+1). Heat on the steam bath until the 
calcium oxalate is dissolved, add 25 ml of a saturated solution of 
ammonium oxalate, dilute to 200 ml, and repeat the precipitation as 
previously described. Collect the precipitate on a paper of close 
texture, and wash with ammonium oxalate-oxalic acid wash solution. 
Transfer the paper and precipitate to a platinum crucible with cover, 
and heat in the uncovered crucible until carbon has been destroyed. 
Cover the crucible, ignite at 1,200° C, and weigh as CaO+Mn,Q,. 
The crucible should be reheated and the weight checked because of 
possible absorption of moisture during the first ten orm) 

Part of the manganese in the original sample be found in the 


ignited calcium oxide. To correct for this, dissolve the ignited residue 


in 50 ml of diluted nitric acid (1+-4), add 5 ml of 85-percent phosphoric 
acid to the solution, and determine the manganese colorimetrically by 
the periodate method described under determination of manganese 
oxide (p. 622). Calculate the manganese to Mn,Q,, and subtract from 
the weight of the ignited residue. The amount of Mn,O, that accom- 
panied the calcium oxide in the case of Sample 56a (MnO=0.18 per- 
cent) was only 0.1 mg, and with Sample 120 (MnO=0.027) it was 
only 0.03 mg. Therefore, no serious error is introduced if the correc- 
tion for manganese is neglected in the analysis of materials containing 
less than 0.2 percent of manganese oxide. The averages of four closely 
checking results obtained by this method for total CaO in Samples 56, 
56a, and 120 were 44.9, 45.62, and 49.67 percent, respectively, as 
ee the certificate values, 44.8, 45.55, and 49.62 percent, 
respectively. 
he proper conditions for the precipitation of calcium in dilute acid 
oe can also be obtained by the method previously described by 
e authors. 


7. DETERMINATION OF MAGNESIA 


Transfer 2.00 g of the sample to a 250-ml beaker, cover, add 15 ml 
of diluted hydrochloric acid (2+1) and 5 ml of nitric acid, and boil 
gently for 10 to 15 minutes. Remove the beaker from the source of 
heat, add 6 ml of diluted sulfuric acid (14-1), remove the cover, 
and evaporate until fumes of sulfuric acid appear. Cool slightly, 


et 


* Prepared by dissolving 0.10 g of brom phenol blue reagent in 1.5 ml of 0.1 N sodium hydroxide and 


aluting to 25 ml. 
%« The Analysis of Phosphate Rock, J. Assn, Official Agr. Chem. 8, 203 (1924). 
57763—38—4 
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wash down the inside surface of the beaker with a jet of water, and 
again evaporate until fumes of sulfuric acid appear. Cool, add 10 
ml of water, stir thoroughly, and digest on the steam bath for 10 tp 
15 minutes. Remove from the steam bath, add 100 ml of 95-percen; 
alcohol, stir so that the calcium sulfate is well dispersed thro 

the liquid, and then allow to stand for 30 minutes or longer. Filte 
by means of suction through a tight plug of filter-paper pulp, using 
a Gooch crucible, carbon funnel, or Biichner funnel, ae wash fiye 
times with 5-ml portions of 95-percent alcohol containing 1 ml of 
sulfuric acid per 100 ml. 

Evaporate the alcoholic filtrate as far as possible on the steam bath, 
Transfer the solution to a 300-ml Erlenmeyer flask, dilute to 75 to 
100 ml, and add 2 g of citric acid and 15 ml of a 25-percent solution of 
diammonium phosphate, (NH,)2HPO,. Add ammonium hydroxide 
until the solution is alkaline to litmus, and then add 10 ml in excegs, 
Add 5 to 10 glass beads, tightly stopper the flask, and shake ona 
shaking machine for at least 1 hour. Allow to stand in a cool place 
for 4 hours or preferably overnight. Filter through a tight paper 
containing a little paper pulp, and wash thoroughly with diluted 
ammonium hydroxide (1+19), containing 50 g of diammonium phos- 
phate per liter. Pass 25 ml of hot diluted hydrochloric acid (1+19) 
through the paper into the flask in which the precipitation took place, 
transfer the solution to a 150-ml beaker, and wash the paper and 
flask thoroughly with more of the diluted acid. To the solution ina 
volume of 50 to 75 ml and containing no glass beads, add 0.5 ml ofa 
25-percent solution of diammonium phosphate, cool, and then add 
ammonium hydroxide, while stirring, until the solution is alkaline to 
litmus. Stir for a few minutes, then add 3 to 4 ml of ammoniun 
hydroxide and allow to stand for 4 hours or overnight. Transfer the 

recipitate to a small filter, and wash with diluted ammonium hydrox- 
ide (1+19). Transfer the paper and precipitate to a platinum cruc- 
ble, ignite slowly at a temperature below 900° C until the carbon is 
burned, preferably in a muffle furnace with ypromenne control, and 
then at about 1,100° C for 1 to 2 hours. Cool and weigh. 

The residue consists of Mg.P,O, and possibly Mn,P,0, and Ca, 
(PO,)2. If the alcoholic filtrate was clear, the tricalcium phosphate 
will not exceed 0.3 mg and can be neglected unless very accurate Ie- 
sults are desired.* e correction for mangeneen is made as follows: 
Dissolve the residue in 10 ml of diluted sulfuric acid (1+9), transfer 
the solution to a 300-ml Erlenmeyer flask, and add 50 ml of diluted 
nitric acid (1+3), 2 ml of 85-percent phosphoric acid, and 0. ae 
potassium periodate, KIO,. Boil for 15 to 20 minutes, cool, and dilute 
to a convenient volume. In another flask containing the same 
amounts of the reagents treated in a similar way, match the color by 
adding a standard solution of potassium permanganate. From the 
volume of the solution of permanganate required, calculate the weight 
of Mn,P,O, in the residue. Subtract this weight from the total F 
weight, and regard the difference as Mg,P,0;, which contains 36.2 
percent of MgO. 


% If a correction for calcium is to be made, dissolve the residue in 10 ml of diluted sulfuric acid (1+9, 

© te to about 5 ml, and then add 25 ml of 95-percent alcohol. Allow to stand for 1 to 2 hours, filter on 
paper, and wash with alcohol. Dry the in the funnel, and dissolve the precipitate of calcium 

sulfate in 20 ml of hot diluted hydrochloric acid G14-99).. Add a few crystals of ammonium oxalate, best, 

render the solution faintly ammoniacal, and allow to stand for about 1 hour. Filter, 

CaO. Calculate to tricalcium phosphate, Ca;(P0O,):, and subtract from the 

pgm FE ns oe hn coy ee me sane 

Before applying the periodate oxidation, filtrate must be evaporated heated 

acid appear, and then treated with nitric acid to make sure that reducing substances are absent. 
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8. DETERMINATION OF FLUORINE 


tnd | 


Fuse 0.50 to 1.00 g of sample with 4.0 g of sodium carbonate in a 
platinum crucible. ° a < the molten mass at a red heat for 5 minutes, 
then grasp the hot crucible with tongs, and immediately immerse the 
lower half in cold water. Tap the cold crucible to loosen the solidified 
mass, and transfer the cooled melt to a 125-ml Claissen flask (see 
fig. 1). Scrub the crucible and cover with a policeman and a little 
hot water, and transfer the washings to the flask, which rests in a 4-cm 
hole in an asbestos gauze EH, and is fitted with rubber stoppers, a de- 
livery bulb A for admit- 
ting acid and water, a 
thermometer D, and a 
condenser B, which de- 
livers the distillate into a 
400-ml beaker C. The 
end of the condenser is 
immersed in 20 ml of a 
20-percent solution of 
sodium hydroxide con- 
tained in the beaker C. 

Add 0.05 to 0.1 g of 20- 
to 40-mesh flake graph- 
ite or granulated coke to 
prevent bumping,” close 
the flask with thestopper 
containing the thermom- 
eter and the delivery ¢ 
tube A, admit 30 ml 
of diluted sulfuric acid 
(1+1), and cautiously 
heat the contents of 
the flask until frothin 
ceases. A good Tirrill 
or similar burner is a 
satisfactory source of 
heat. Increase the heat 
until the liquid in the 
flask boils vigorously, 
and continue the boiling a 
until the temperature of Figure 1.—Apparaius for distilling fluorine. 
the liquid has risen to 
160° C. Then without removing the source of heat, admit water 
through the tube A at such a rate that the temperature of the liquid 
in the flask remains between 160 and 170° C. Collect 300 ml of dis- 
tillate between these temperatures. 

Evaporate the alkaline distillate to a volume of 250 ml, and add 3 ml 
of a 10-percent solution of sodium chloride and 2 drops of bromphenol 
blue indicator.” Add dilute nitric acid until the color changes to 


wh and then sodium hydroxide until the color changes just to 
ue. 

















———— LL 
* This means of preventing bum by W. E. 
national Agri salen cs ing a ping: was gusgested y Dickinson, Chief Chemist of the Inter 


orporati Ik, 
ing to spared by dissolving 0.10 g of bromphenol blue reagent in 1.5 ml of 0.1 N sodium hydroxide and dilut- 
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Now add 2 ml of diluted hydrochloric acid (1+1) and 5.0 g of soliq 
lead nitrate, Pb(NO;)2, and heat on the steam bath. As soon as 
the lead nitrate is in solution, add 5.0 g of solid sodium acetate, 
NaC,H;0;.3H,0, stir vigorously, and digest on the steam bath for 
30 minutes with occasional stirring. Allow to stand overnight gj 
room temperature, and decant the solution through a paper of cloge 
texture (589 S & S blue ribbon or 42 Whatman). Wash the pm. 
cipitate, beaker, and paper once with cold water, then four to fiye 
times with a cool saturated solution of lead chlorofluoride, 8 and then 
once more with cold water. 

Transfer the precipitate and paper to the beaker in which the 
precipitation was made, stir the paper to a pulp, add 100 ml of diluted 
nitric acid (1+19), and heat on the steam bath until the precipitate 
is dissolved (5 minutes is ample). Then add a slight excess of a 0.2 V 
solution of silver nitrate,” digest on the steam bath for 30 minutes, 
cool to room temperature while protected from the light, filter, wash 
with cold water, and determine the silver nitrate in the filtrate by 
means of a standard solution of potassium thiocyanate,” using 5 m 
of a solution of ferric alum as the indicator.” Subtract the ameuih 
of silver nitrate found in the filtrate from that originally added, 
The difference will be that required to combine with the chlorine in 
the lead chlorofluoride, and from this the fluorine may be computed. 
One ml of 0.2 N nitrate is equivalent to 0.00380 g of fluorine. 

The results obtained for fluorine by the foregoing procedure in the 
analysis of synthetic mixtures and Standard Phosphate Rocks 56a 
and 120 are shown in table 1. 


TaBLE 1.—Resulis obtained for fluorine with and without preliminary fusion with 
sodium carbonate 








Fluorine found 


Fluorine pone 


























® NBS Standard Sample 79 containing 94.8 percent of CaF:. 
» Certificate value. 


28 To prepare lead chlorofluoride wash solution: (a) Dissolve 10 g of lead nitrate in 200 ml of vaier 
solve 1.0 g of sodium fluoride in 100 ml of water, and add 2 ml of hydrochloric acid; and, (c) mix solutions 
and (b). Allow the precipitate to settle and decant the su tant liquid. Wash four or five times 
200 m] of water by decantation, and then add about 1 liter o Gy water » the “ey voy anes alee nae 
for 1 hour or longer with occasional stirring. Pour through a filter and use a Sn 
more a pel to de id recipitate of lead chlorofluoride in the 1 beaker and pe ney more soe 
prepared as ni 

Rt. A amen cit by precipitating a measured portion with hydrochloric acid and weighing the silver 

c 

% Standardized by comparing with the standard — of silver nitrate under the same conditions 3 
obtain in the method. A convenient concentration is 0.1 N. 

31 A cold saturated water solution of ferric alum by ud chlorides) to which has been added sufficient 
colorless nitric acid to bleach the brown color of the water solution. 
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9. DETERMINATION OF SILICA 


Fuse 0.500 g of the dried sample with 5 g of sodium carbonate, 
transfer the cooled melt to a platinum dish, leach with hot water, and 
filter when disintegration is complete. Return the insoluble residue 
to the dish by means of a jet of hot water, add 50 ml of a 2-percent 
solution of sodium carbonate, boil for a few minutes, filter, and wash 
thoroughly with hot water. Reserve the combined filtrates (A). 
Transfer the insoluble residue to a platinum crucible, ignite, fuse with 
3g of sodium carbonate, leach the cooled melt with hot water, and 
filter when disintegration is complete. Reserve the filtrate (B). 

By means of a jet of water transfer the insoluble residue to the 
platinum dish, then add 1 g of boric acid, 5 ml of nitric acid, and 10 ml 
of 60-percent perchloric acid. Heat until nearly all the perchloric 
acid is volatilized, cool, drench the residue with 5 ml of hydrochloric 
acid, and then add 50 ml of hot water. Digest on the steam bath for 
5 minutes, and collect the precipitate on a small filter. Discard the 
filtrate and reserve the precipitate (C). 

To the combined filtrates, (A) and (B), which should have a volume 
of about 300 ml, add 1 g of zine oxide dissolved in 20 ml of diluted 
nitrie acid (1+4), boil for 1 minute, filter, and wash the precipitate 
thoroughly with hot water. Reserve the precipitate (D). Add a few 
drops of methyl red indicator solution to the filtrate, nearly neutralize 
with nitric acid, and evaporate to a volume of 200 ml. Finish the 
neutralization of the concentrated solution by adding diluted nitric 
acid (1+9) until the color is a very faint pink. Now add 1.0 g of 
zinc oxide dissolved in ammonium hydroxide and ammonium car- 
bonate,”” and boil in a covered platinum dish until the odor of am- 
monia is entirely gone. After all the ammonia is expelled, add about 
50 ml of warm water, stir, allow to stand for a few minutes, filter, and 
wash the precipitate with cold water. Discard the filtrate and reserve 
the precipitate (£). 

ith the aid of a jet of diluted hydrochloric acid (1+19), transfer 
the precipitates (C), (D), and (Z£) from the papers to the dish in which 
the last precipitation was made. Ignite these papers and those used 
in the previous filtrations, and add any residue so obtained to the 
contents of the dish. Add 25 ml of hydrochloric acid, and evaporate 
to dryness on the steam bath. Drench the residue with 10 ml of 
hydrochloric acid, and then add 100 ml of hot water. Digest on the 
steam bath for 15 minutes, filter, and wash with hot diluted hydro- 
chloric acid (1+19), and then with hot water. Return the filtrate 
and washings to the dish in which the evaporation was made, add 
10 ml of sulfuric acid, and evaporate until fumes of sulfuric acid are 
evolved. Allow to cool, add 100 ml of water, digest for 5 minutes on 
the steam bath, filter, and wash with hot water. 

Place the two papers containing the silica in a platinum crucible, 
char the paper without flaming, and finally ignite at 1,200° C. Cool 
in @ desiccator, weigh, and determine the silica by treatment with 
hydrofluoric and sulfuric acids in the usual manner. 

Synthetic mixtures consisting of 0.5000 g of National Bureau of 
Standards Limestone 1A+-0.050 g of fluorspar (CaF 2)+0.25 g of tri- 


that 
, mater 1.0 g of zine oxide and 2.0 g of ammonium carbonate to a small beaker, add 20 ml of water and 
ml of ammonium hydroxide (sp gr 0.90), and digest on the steam bath until a clear solution is obtained. 








calcium phosphate were prepared and analyzed for silica by the fore. 
going procedure. ‘Two determinations yielded 0.0713 and 0.0709 g of 
silica, respectively. The total content of silica was 0.0714 g. In the 
analysis of Standard Samples of Phosphate Rocks 56a and 120,77 
and 11.0 percent were obtained, as compared with the certificate values 
7.6 and 11.0 percent, respectively. 


10. DETERMINATION OF MANGANESE OXIDE (MnO) 


Transfer 0.500 g of sample to a 300-m! Erlenmeyer flask and add § 


5 ml of nitric acid and 15 ml of 60-percent perchloric acid. Heat to 
boiling, and continue the boiling until fumes of perchloric acid begin 
to be evolved. Remove from the source of heat, cool somewhat, and 
add 10 ml of diluted hydrochloric acid (1+-1). Now boil the solution 
until most of the perchloric acid has been volatilized. Reduction of 
the volume to about 5 ml is satisfactory. 

Cool the solution and add 50 ml of warm diluted nitric acid (1+4) 
in which 0.5 g of boric acid has been dissolved. Digest for a few 
minutes, cool, filter through a paper of loose texture, and wash the 
flask and paper a few times with warm water. To the filtrate and 
washings in a 300-ml Erlenmeyer flask add 5 ml of 85-percent phos- 
phoric acid and 0.5 to 1.0 g of potassium periodate. Boil for 10 min. 
utes, cool, dilute if necessary, and determine the amount of manganese 
in the solution by comparing the color with that obtained by dissolving 
and similarly treating a material of known manganese content. It 
is obvious that it is desirable, but not necessary, to select a standard 
whose manganese content is approximately that of the unknown. The 
National Bureau of Standards Standard Samples of Phosphate Rocks 
56, 56a, and 120 are satisfactory standards for materials containing 
up to 0.5 percent of MnO. 

No significant quantity of manganese remains in the insoluble 
residue if the National Bureau of Standards Standard Samples 56, 56a, 
and 120 are decomposed by this procedure. 


11. DETERMINATION OF THE ALKALIES 
Determine the alkalies by the J. Lawrence Smith method.* 
12. DETERMINATION OF SULFUR 


_Sulfur is seldom determined. Detailed directions are therefore not 
ven. 

Total sulfur was determined by (1) digesting a 2-g sample with con- 
centrated nitric acid containing AS» (2) evaporating with hydro- 
chloric acid to remove the nitric acid; (3) digesting the residue with 
hydrochloric acid (1+1) and filtering; (4) adjusting the acidity to | 
to 2 ml of hydrochloric acid in 100 ml of solution; (8) precipitating the 
sulfate by adding barium chloride; (6) filtering and washing; 
igniting the barium sulfate; and (8) correcting for sulfate obtained 
from reagents. No sulfur was found in the siliceous residue resulting 
from the evaporation with hydrochloric acid. 
~~" W F. Hilleb ork 

W F. Hillebrand and G. E. F. Lundell, Applied Inorganic Analysis, p. 786 (J. Wiley & Sons, New Y i 


N. Y.. 1929 ed.); G. E. F. Lundell and J. I. Hoffman, Outlines of Methods of C Pp. 
and 185 (J. Wiley & Sons, New York, N. Y., 1938 ed.). 
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The percentages of total sulfur (as SQ,) found in Phosphate Rocks 
56a and 120 were 0.73 and 0.32, respectively. 

Pyritic sulfur was determined by transferring a 5-g sample to an 
Erlenmeyer flask, acting 40 ml of diluted hydrochloric acid (1+1), 
boiling for 5 minutes, filtering, washing the residue thoroughly with 
diluted hydrochloric acid (1+-99), discarding the filtrate and washings, 
and then treating the residue as was done for total sulfur. 

No significant amount of pyritic sulfur (less than 0.005 percent of 
S) was found in Phosphate Rock 56a and only 0.01 percent in Phos- 
phate Rock 120. 

Soluble sulfide sulfur, that is, sulfide sulfur that can be evolved by 
treatment with hydrochloric acid, was determined by the method that 
is used for this determination in cements.** The sample is treated 
with hydrochloric acid and stannous chloride in a suitable flask, and 
the evolved gases are passed into an ammoniacal solution of zinc 
sulfate. The stannous chloride prevents oxidation of sulfides by the 
oxidized constituents of the rock during the evolution. The sulfide 
which is absorbed in the ammoniacal zinc sulfate is determined by 
acidifying the solution with hydrochloric acid, and titrating with a 
standard solution of iodine. 

Only 0.01 percent of sulfide sulfur (as S) was found in Phosphate 
Rock 56a and 0.05 percent in Phosphate Rock 120. The foregoing 
results indicate that most of the sulfur in these two rocks is present as 
soluble sulfates. 


13. SPECTROCHEMICAL ANALYSIS (BORON) 
(a) QUALITATIVE EXAMINATION 


Spectrochemical analyses of the two samples were made by B. F. 
Soribsler of the spectroscopy section of this Rareae: The arc spectra 
were photographed in the region 2470 to 3300 A, using carbon elec- 
trodes, and in the region 3200 to 4600 A, using copper electrodes. 
The spectra were examined for the sensitive lines of Ag, Al, As, Au, B, 
Ba, Be, Bi, Ca, Cb, Cd, Ce, Co, Cr, Cu, Fe, Ga, Ge, Hf, Hg, In, Ir, K, 
Li, Mg, Mn, Mo, Na, Ni, Os, Pd, Pb, Pt, Rh, Ru, Sb, Sec, Si, Sn, Sr, 
Ta, Th, Ti, Tl, U, V, W, Y, Zn, and Zr. No special effort was made 
to detect such elements as fluorine or sulfur, since other tests had 
already shown their presence. The elements found are as follows: 
Tennessee Phosphate Rock 56a. 

Major constituents: Ca, Al, Fe, Mg, Si. 

Minor constituents: Mn, Na, K, Cr, Sr, Cu, Ti. 

eae vie nannccceeescavcowa Pb, B, Sn, Ni. 


Florida Phosphate Rock 120. 


Major constituents: Ca, Al, Fe, Mg. Si. 
Minor constituents: Mn, Na, K, Cr, Sr, Cu, Ti, Mo, Ni, V. 
a i ea Pb, Sn, B. 


(b) QUANTITATIVE DETERMINATION OF BORON 


Boron is the only constituent that was quantitatively determined 
by the spectrochemical method. This was ine by comparison of the 
spectra of the two samples with the spectra of a set of standard mix- 
tures photographed under comparable conditions. The standards 


“H. A. Bright, J. Research NBS 18, 137 (1937) R P96s. 
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were prepared to give a graded series of known concentrations of 
boron, using as a base a sample of tricalcium phosphate practical} 
free from boron. The sensitivity of the test for boron is approx. 
mately 0.0001 percent. The precision of the method decreases 4s 
the concentration decreases, being within +10 percent of the amount 
present for concentrations above 0.005 percent of boron and within 
+20 percent below 0.005 percent. The boron reported as B,O, was 
as follows: 

Tennessee Phosphate Rock 56a: 0.005 percent of B,O3. 

Florida Phosphate Rock 120: Less than 0.001 percent of B,Q3. 

It is interesting to note that the spectrochemical tests indicated 
higher percentages of manganese, copper, and boron in the Tennessee 
Rock than in the Florida Rock, and higher percentages of nickel, 
molybdenum, vanadium, tin, and chromium in the Florida than in the 
Tennessee Rock. This is in agreement with the following percentages 
of these constituents found by chemical methods.*® 

Tennessee Rock 56a: MnO, 0.18; CuO, 0.002; B,O;, 0.005; NiO, 
0.0005; MoOs, 0.0002; V,0;, none found; 
SnO,, 0.001; Cr.Qs, 0.001. 

Florida Rock 120: MnO, 0.027; CuO, 0.001; B,O;, 0.002; NiO, 
0.003; MoO;, 0.002; V,0;, 0.015; SnO,, 
0.002; Cr2O3, 0.01. 


The ability to distinguish between such small percentages suggests 
the possibility of determining these and other minor constituents 
quickly and accurately by spectrochemical methods, provided that 
suitable standards are available. 


IV. SUMMARY OF THE RESULTS OBTAINED IN THE 
COOPERATIVE ANALYSIS OF PHOSPHATE ROCKS 


Tables 2 and 3 give a summary of the results obtained in the analysis 
of the National Bureau of Standards Phosphate Rocks 56a and 120, 
The analysts mentioned in this tabulation have rendered valuable 
assistance in this work, and the authors express their deep appreciation 
to them. 


% The boron was determined chemically by Lewis F. Rader, Jr., and W. L. Hill, Fertilizer Research 
Division, Bureau of msg | A and Soils, U. S. Department of Agriculture, by the use of a method involving 
distillation and titration. The method will be published in the Journal of Agricultural Research. (Now 


in press.) 
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Lundell 


—Results obtained in the rative analysis of the National Bureau of 
Standards Standard Sample of Phosphate Rock 56a 


[All results are based on a sample dried for 1 hour at 105° C4] 


TABLE 2. 





Analyst* 
(by number) 


as SO; 


Total FeO; 
Total Al,O; 
| Total sulfur 











©33. 05 



































Averages...} & 33.01 














s Seven analysts reported an av: of 0.53% of moisture, obtained by drying for 1 hour at 105° C. The 
results varied from 0.36 to 0.64%. ot much cance is to be attached to these values because the 
moisture content changes with atmospheric conditions. 

» Gravimetric.—Precipitation with molybdate followed by two precipitations with magnesia mixture 
(J, Assn. Official Agr. Chem 8, 195 (1924). 

¢ Gravimetric.—T wo Pag grog with magnesia mixture without preliminary precipitation with molyb- 
date (J. Research NBS 19, 59 (1937) RP1010). 

4 AOAC, volumetric. 

e AOAC, gravimetric. 

Gravimetric, — Single precipitation with magnesia mixture after molybdate precipitation (J. Assn. 


Official Agr. Chem, 8, 194 (1924)). 

«Special work at the Nationa) Bureau of Standards, and a careful consideration of all the methods em- 
ployed indicate that 32.90 should be regarded as the most probable percentage of P30; in this sample. 

b Special work at the National Bureau of Standards indicates that an uncertainty of at least +0.2% 
should be assigned to this value. 

Analyst 1 also found 0.002% of CuO, 0.0005% of NiO, 0.0002% of MoOs:, 0.001% of SnOs, 0.001% of 
Cr03, and no V20; in a 10-g sample. 

Analyst 7 also reported 0.005% of B:03. 3B. F. Scribner of the spectroscopy section of the National Bu- 
reau of Standards found 0.065% of Bz03, s emically. Analyst 9 also reported 1.56% of COs 
0.%% of organic matter, 0.0021% of As,O3, and 0.39% of combined water. 


*ANALYSTS, BY NUMBER, AS LISTED IN TABLE 2 


1. James I. Hoffman, National Bureau of Standards, Washington, D. C. 

2. E. W. Magruder, Chief Chemist, and W. A. Ryder, F. 8. Royster Guano Co., 
Norfolk, Va. 

3. W. Catesby Jones, Chief Chemist, John H. Elder, and James R. Lindsay, 
Commonwealth of Virginia, Department of Agriculture and Immigration, 
Richmond, Va. 

4, C. A. Butt, Chief Chemist, and C. M. Cartledge, International Agricultural 
Corporation, East Point, Ga. 

5. W. E. Dickinson, Chief Chemist, International Agricultural Corporation, 
Norfolk, Va. 

6. P. McG. Shuey, Shuey & Co., Savannah, Ga. 

7. W. L. Hill, L. F. Rader, Jr., T. H. Tremearne, H. L. Marshall, and D. S. 
Reynolds, phosphate section, Fertilizer Research Division, Bureau of 
Chemistry and Soils, U. 8. Department of Agriculture, Washington, D. C. 

8. John B. Smith, Chief Chemist, W. L. Adams, and E. J. Deszyck, Rhode 
Island State College Agricultural Experiment Station, Kingston, R. I. 

9, Egbert Janes, Chief Chemist, and H. H. Edwards, International Agricultural 

ration, Mulberry, Fla. 

10. L. E. Dupont, Chief Chemist, International Agricultural Corporation, Wales, 

enn. 
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TABLE 3.—Results obtained in the cooperative analysis of the National Bureay, of 
Standards Standard Sample of Phosphate Rock 120 
[All results are based on a sample dried for 1 hour at 105° Cs] ° 
{ 
é|é 
Analyst* oa = 2s ne 
(by number) S 3 3 a5 
& e&elela NV 
Ruy * | %| * 
1 {35:1 0.8 0.85} 0.87 
© 35. 38 
Fhe ea ng. Ek, aces aera 1 
eatery » 35.24} 90). 80}_._._- 
Spee. © 35. 46)_...._|_.__.. 92 
65... ss © 35.40|_...-.|_..-.- 92 
“Wewemeom 4 35, 43|_---.-|------ 85 
Visakataatibetbe £35.25] .89| .82|...... 
atone (pe) 
SES OR 435.43) 87} .89) .86 
WD. conceit ¢35.37) .92| .92)_._..- : ~ 
Averages...| © 35.33} .89| .87| .89| .80 40.62} .14| 3.76] » 7.40] .033| .14) .09] .a2 hs 
0 
fre 




















® Seven analysts reported an a of 0.53% of moisture, obtained by drying for 1 hour at 105°C. The 
results varied from 0.43 to 0.66%. Not much significance should be attached to these values, because the pr 


moisture content changes with atmospheric conditions, Wi 
» Gravimetric.— mpltation with molybdate followed by two precipitations with magnesia mixturs m: 
(J. Assn. Official Agr. em. 8, 195 (1924)). 


¢ Gravimetric.—T wo prasreriene th magnesia mixture without preliminary precipitation with ; 
molybdate (J. Research NBS 19, 59 (1937) RP1010). is 
4 AOAC,. volumetric. 
e AOAC, gravimetric. : 
f Gravimetric.—Single precipitation with magnesia mixture after molybdate precipitation (J. Assn. di 
Official Agr. Chem. 8, 194 (1924). 
« Special work at the National Bureau of Standards,°and a careful convideration of all the methods em- de 
ployed indicate that 35.20 should be regarded as the most b meeey percentage of P30; in this sam f 
» Special work at the National Bureau of Standards indicates that 7.6 represents more nearly the true iv 
— content of this material, and that an uncertainty of at least +0.2 percent should be assigned to this 
value. 
Analyst 1 also found 6.061% of CuO, 0.003% of NiO, 0.002% of MoOsz, 0.015% of V20s, 0.002% of 8n0; 
and 0.01% of CrsO3. Analyst 7 also reported 0.001% of B30. . F. Scribner of the epeeerenen ane 
the National Bureau of Standards found (0.001% of B203, spectrochemically. Analyst 9 reported 
1.96% of COs, 0.32% of organic matter, 0.0018% of As;03, and 0.12% of combined water. ) 


I 
ID 


*ANALYSTS, BY NUMBER, AS LISTED IN TABLE 3 7 


1. James I. Hoffman, National Bureau of Standards, Washington, D. C. 
2. . Chief Chemist, and W. A. Ryder, F. 8. Royster Guano Co, 
Norfolk, Va. 

3. W. Catesby Jones, Chief Chemist, John H. Elder, and James R. Lindsay, VI 
Commonwealth of Virginia, Department of Agriculture and Immigration, 
Richmond, Va. 

4. C. A. Butt, Chief Chemist, and C. M. Cartledge, International Agricultural 
Corporation, East Point, Ga. | 

5. W. E. Dickinson, Chief Chemist, International Agricultural Corporation, tr 
Norfolk, Va. ; 

6. P. McG. Shuey, Shuey & Co., Savannah, Ga. 101 

7. W. L. Hill, L. F. Rader, Jr., T. H. Tremearne, H. L. Marshall, and D.8. — at 
Reynolds, phosphate section, Fertilizer Research Division, Bureau f F fy, 
Chemistry and Soils, U. 8. Department of Agriculture, Washington, D. ©. 

8. John B. Smith, Chief Chemist, W. L. Adams, and E. J. Deszyck, gre 
Island State College Agricultural Experiment Station, Kingston, R. L no 

9. Egbert Janes, Chief Chemist, and H. H. Edwards, International Agricultunl 3 
go isinaer Mulberry, Fla. 1937 

10. L. E. Dupont, Chief Chemist, International Agricultural Corporation, Wales, a 
Tenn. int 
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MAXIMUM USABLE FREQUENCIES FOR RADIO SKY- 
WAVE TRANSMISSION, 1933 TO 1937! 


By Theodore R. Gilliland, Samuel S. Kirby, Newbern Smith, and 
Stephen E. Reymer 





ABSTRACT 


Graphs are presented to indicate the maximum usable frequencies for sky- 
wave transmission from June 1933 to December 1937. The graphs are given for 
March, June, and December to show equinoctial, summer, and winter conditions 
for each year. The factors which must be considered in deriving these graphs 
from vertical-incidence ionosphere measurements are outlined briefly. he 
principal factors are Snell’s law, variation of virtual height of the ionosphere 
with haan the curvature of the ionosphere, and the effect of the earth’s 
magnetic field. 

he method of applying these graphs to simple and complex transmission paths 
is indicated. 

The effects of sporadic EF reflections, absorption, and scattered reflections are 
discussed briefly. 

It is pointed out that the graphs and the ionosphere data from which they were 
derived may be used to estimate future diurnal, seasonal, and long-time variations 
of maximum usable frequencies. 
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I. INTRODUCTION 


The graphs of maximum usable frequencies for radio sky-wave 
transmission presented here were prepared from the vertical-incidence 
ionosphere measurements made by the National Bureau of Standards 
at Washington.? These graphs, for any particular month, are derived 
from the average condition of the ionosphere for that month. The 
graphs give results representative of the winter, summer, and equi- 
noctial periods from June 1933 to December 1937, inclusive. 


1A nary version of this paper was presented as r C6 of the meeting in London in November 
ae the special committee on radio-wave Propagation. . 
Theodore R. Gilliland, Samuel 8. Kirby, Newbern Smith, and Stephen E. Reymer, Characteristics of 
pd ere and their application to radio transmission. J. Research NBS 18, 645 (1937) RP1001; Proc, 
1174 ( ote 823 (1937), and a series of monthly reports beginning in Proc. Inst. Radio Engrs. 25, 
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The maximum usable frequency for any distance is defined as th 
highest frequency which can be used for radio sky-wave transmission 
over the given distance. Waves of frequencies higher than this 
penetrate through the ionosphere and are not returned to the earth: 
waves of lower frequencies are reflected and may be used for trang. 
mission over this distance. Graphs of maximum usable frequencies 
in terms of distance also give skip distances as a function of frequency, 
The maximum usable frequencies give more directly than skip dis. 
tances information usually desired in connection with radio traps. 
mission; i. e., determination of the maximum frequency which ¢gp 
be used over a given distance is of more general value to an engineer 
than the minimum distance over which communication can be had 
on a given frequency. 


II. RELATION OF VERTICAL TO OBLIQUE-INCIDENCE 
TRANSMISSION 


Radio waves incident obliquely upon the layers of the ionosphere 
are more easily reflected than those incident vertically, that is, they 
can be reflected from regions of smaller ionization density. In gen- 
eral, for a given frequency of the transmitted wave, the larger the 
angle of incidence, the less is the ionization density required for reflec- 
tion. Thus, if a given layer will reflect waves up to a certain fre- 
quency at vertical incidence, it will reflect waves of considerably 
higher frequencies at oblique incidence. 

he angle of incidence is the angle between the ray and the normal 
to the lower boundary of the ionosphere at the point where the ray 
enters the ionosphere. The angle of incidence should not be confused 
with the angle of departure, which is the angle above the horizontal at 
which the ray is projected from the transmitting station. Although, 
in general, a transmitting station may radiate energy at all angles 
above the horizontal, only certain ray-paths, corresponding to a few 
definite angles of departure, are useful for transmission over a given 
distance at a given time. At vertical incidence the angle of incidence 
is zero, and increases for oblique incidence. 

The angle of incidence depends on (1) the distance of transmission 
for one reflection from the ionosphere, and (2) the height of the layer. 
Given these two quantities, therefore, the angle of incidence may be 
calculated. The maximum usable frequency is given roughly, from 
this angle and the vertical-incidence critical frequency, by the follow- 
ing relation, known as the “secant law’’: 


f' =f. sec 0; 


j’=maximum usable frequency, } 

f.=critical frequency or highest frequency reflected at vertical 
incidence from a given layer, and 

6,=angle of incidence. 


The given distance will be within the “skip-zone” for all higher 
frequencies. 

The secant law does not take several important factors into account. 
The first and most important of these is the fact that the virtual 


where 
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height of a given layer of the ionosphere increases with frequency, 
especially near the critical frequency. As a result, it is difficult to 
choose the proper height to be used in calculating the angle of inci- 
dence for a given distance. A method whereby, for a flat ionosphere, 
this difficulty may be overcome has been described by Smith.® 

The second factor which must be considered is the curvature of the 
ionosphere. This introduces another modification in the calculation 
of the angle of incidence for a given distance and a given height of 
reflection. This results in a higher maximum usable frequency over a 
given distance than would be calculated on the basis of a flat iono- 
sphere. 
Stil another factor that must be taken into account is the effect of 
the earth’s magnetic field. In the presence of the earth’s field a radio 
wave in the ionosphere is split into two components known as the 


JUNE 1933 


20% 
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Ficure 1.—Mazimum usable frequencies, latitude of Washington, D. C., June 1933. 


ordinary ray and the extraordinary ray. These components are re- 
flected from different levels in the ionosphere, and consequently have 
different critical frequencies. At Washington, the vertical-incidence 
critical frequency for the extraordinary ray is about 800 ke/s higher 
than for the ordinary ray, for frequencies in the neighborhood of 
5,000 ke/s, and the separation is greater for lower frequencies. At 
oblique incidence the effect of the earth’s magnetic field complicates 
the calculation of maximum usable frequencies. The maximum 
usable frequency for the extraordinary ray is always greater than for 
the ordinary ray. Except for short distances, the difference is not 
great, and decreases with increasing frequency and distance. For 
practical communication problems it is sufficiently accurate to make 
the calculations for the ordinary ray. 

Methods of obtaining maximum usable frequencies, taking into 
account the earth’s magnetic field and the curvature of the ionosphere, 
have been presented by Smith.‘ These methods have been used in 
preparing figures 1 to 14. 

i 

niaion. J Research NBS 10, 20 C08), 
‘Newbern 8 ? 
transmission. 


oe ionosphere measurements to oblique-incidence radio trans- 
mith, The cre of vertical-incidence ionosphere measurements to oblique-incidence radio 
J. Research S 20, 683 (1938) RP1100; Proc. Inst. Radio Engrs. 
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Fiaure 4.—Marimum usable frequencies, latitude of Washington, D. C., June 1934. 
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2.—Mazimum usable frequencies, latitude of Washington, D. C., December 
1933. 
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3.—Mazimum usable frequencies, latitude of Washington, D. C., March 
1984. 
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III. TRANSMISSION PATH 





The geographical part of the ionosphere which controls long-distance 
high-frequency radio transmission is that part traversed by the wave 
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Figure 5.—Mazximum usable frequencies, latitude of Washington, D. C., December 
1934. 


in passing from the transmitter to the receiver. The neareast part of 
the ionosphere thus traversed may be as much as 1,000 to 2,000 km 
from the transmitter, and a similar distance from the receiver. Iono- 
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Figure 6.—Mazimum usable frequencies, latitude of Washington, D. C. March 


1935. 


sphere conditions at this point may be considered to be the same as 
conditions at the transmitter or receiver, at the same local time, pro- 
vided no great differences of latitude exist. Consequently, depending 


















IN MEGACYCLES 


FREQUENCY 


FIGURE 


IN MEGACYCLES 


FREQUENCY 


632 Journal of Research of the National Bureau of Standards vn» 


Figure 8.—Mazimum usable frequencies, latitude of Washington, D. C., December 
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7.—Mazimum usable frequencies, latitude of Washington, D. C., June 1935, 
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on the time of day, different frequencies might be necessary for trap. 
mission in different directions. 
Take, for example, stations at Washington, D. C., North Platt 
Nebr., and San Francisco, Calif., in the eastern, central, and westen, 
arts of the United States, respectively. The distances betwee 
ashington to North Platte and between North Platte and Sy 
Francisco are approximately 2,000 km. Also take the time as 07}; 
local time at North Platte in December 1937 (see fig. 14). Then th 
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Figure 11.—Mazimum usable frequencies, latitude of Washington, D. C., December 
1936. 


local time at the reflection point between North Platte and Washing- 
ton would be about 0800 and the local time at the reflection pomt 
between North Platte and San Francisco would be about 0630. By 
interpolating linearly between the hs for successive times in figure 
14 it is evident that, at this time, t a lietiont frequency (about ge 


kc/s) which can be used between North Platte and Washington 
much greater than the highest frequency (about 8,300 ke/s) which 
can be used over the different path between North Platte and San 
Francisco. This does not mean, however, that different frequencies 
should be used for two-way communication over the same pa he 
as North Platte to Washington and Washington to North Platte, 
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Figure 12.—Mazimum usable frequencies, latitude of Washington, D. C., March 
1937. 
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Ficurg 13.—Mazimum usable frequencies, latitude of Washington, D. C., June 1987. 
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since the maximum usable frequency is the same for transmission jp 
opposite directions over the same path. 

For single-reflection transmission the controlling portion of the 
ionosphere is halfway between the terminal points. kK or multireflec. 
tion transmission the ionosphere along the entire path, except for that 
portion within several hundred kilometers of the terminal points 
must be considered. Because of large differences in local time and 
latitude encountered in long transmission paths involving more than 
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Figure 14.—Mazimum usable frequencies, latitude of Washington, D. C., December 

19387. 


one reflection, widely different conditions sometimes prevail over 
different parts of these paths. In such cases the frequency will have 
to be lowered to satisfy conditions in the part of the path in which the 
critical frequency is lowest. Under such conditions absorption 1 
likely to occur in the part of the path with the higher critical fre- 
quencies. Normally the ionosphere does not vary greatly over & 
latitude range somewhat exceeding that of the United States. — 
The maximum possible distance for single-reflection transmission, 
corresponding to zero angle of departure, is about 2,400 km for 
E-layer transmission, and about 3,500 to 4,400 km for F;-layer trans 
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mission, depending on the virtual height of the layer. Practically, 
it is usually impossible to accomplish high-frequency transmission at 
this zero angle of departure over land because of absorption at the 
earth’s surface. If a practical lower limit of 3%° is assumed for 
the angle of departure over land, the maximum distance for single- 
reflection transmission is about 1,700 km for #-layer transmission and 
2,800 to 3,600 km for F;-layer transmission, depending on the height. 
Single-reflection transmission may often be possible at greater dis- 
tances, while at the same time multireflection transmission over the 
same path may be more efficient. 


IV. GRAPHS OF MAXIMUM USABLE FREQUENCIES 


The maximum usable frequencies given in figures 1 to 14 represent 
monthly average conditions. The graphs for December represent 
winter conditions which hold for several months centered on the winter 
solstice, those for June represent summer conditions which hold for 
several months centered on the summer solstice, and those for March 
represent the spring and fall equinoctial conditions or transition condi- 
tions between winter and summer. The seasonal variations in 
characteristics of the ionosphere are given in detail in the papers 
referred to in footnote 1. 

Day-to-day variations over short periods, such as a few weeks, are 
usually small (less than 10 percent from the average) and are given 
in detail in the papers referred to in footnote 1. Large variations, 
however, do occur during periods of ionosphere storms,®* which are 
usually associated with magnetic disturbances. At such times the 
virtual heights of the F; layer are increased and the critical frequencies 
decreased, resulting in lower maximum usable frequencies. The F, 
region becomes diffuse and unstable, resulting in high absorption and 
poor quality of transmission. 

The graphs are given for single-reflection transmission at the lati- 
tude of Washington. Linear interpolations may be made between 
graphs for successive times. The conditions given in the curves are 
for local time at the geographical part of the ionosphere where the 
waves are reflected. For multireflection transmission, the distance 
scale on the curves must be multiplied by the number of hops. For 
example, for two-hop transmission, 3,000 km on the graphs should be 
read as 6,000 km. For a path which does not involve latitudes differ- 
ing widely from that of Washington, the maximum usable frequency 
is then the lowest one of the several corresponding to the local times at 
the several points of reflection along the path. 

As an example of the use of the graphs, let us consider two-hop 
transmission from Washington to San Francisco, 4,000 km, at 0900 
EST in December 1937. The distance of each hop is 2,000 km, and 
the local times at the points of reflection are 0815 and 0645, at the 
eastern and western reflecting points, respectively. Interpolating 
between the graphs for successive times in e 14 at the 2,000-km 
distance, we obtain the maximum usable frequency to be 21,000 ke/s 
for the first, or eastern hop, and 10,000 ke/s for the second, or western 
hop. Thus 10,000 ke/s is the highest frequency that can be used at 

time for two-hop transmission over the given path. 
ant 8. Kirby, Theodore LS Giltend, Ebert B. Judson, and Newbern Smith, The ionosphere, sun- 


s 
samuel 8. Kirby, Theodore R. Gilliland, Newbern Smith, and Stephen E. Reymer, The ionosphere 
Wler eclipse, and magnetic storm. Phys. Rev. 60, 258 (1936). mn . ‘ ™ 
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The graphs for June show that both the # and F _ are effectiye 
in determining the maximum usable frequencies. For example, 9 
noon the F, layer will be effective out to a short distance, say 500 km 
because of its higher critical frequency; then the £ layer will be effes. 
tive out to 1,600 km because of its lower height and resulting | 
angles of incidence. Beyond about 1,600 km the intensity of ian 
hop E-layer transmission decreases because of the necessarily small 
angles of departure at the transmitting station. The maximum usable 
frequency is then that which will be transmitted by two-hop E., o 
single-hop F,-layer propagation. Either of these conditions gives q 
maximum usable frequency lower than does the single-hop H-layer 
transmission at about 1,600 km. Of the two conditions, the F; layer 
gives the higher maximum usable frequency for average conditions 
at Washington as indicated on the graphs. A similar decrease of 
maximum usable frequency occurs when the F-layer transmission 
must change from one-hop to two-hop propagation at about 3,500 
km. 


V. OTHER TRANSMISSION CHARACTERISTICS 
1. SPORADIC E 


Reflections from the E region are often found at frequencies con- 
siderably in excess of the maximum daytime critical frequency for the 
E layer. These differ in character from the regular refraction phe- 
nomena of the normal £ layer and occur at all seasons, but they are 
more prevalent during the summer months in Washington, sometimes 
persisting for several hours. Sporadic £ reflections seem to be of 
the nature of reflections from a sh boundary, as suggested by 
Kirby and Judson,’ rather than of reatiler refraction from a region of 
intense ionization. Because of the sporadic nature of these reflec- 


tions their effect has not been included in the graphs. Sporadic £, | 


however, accounts for good transmission at higher frequencies than 
those indicated by the graphs for a small percentage of the time and 
at irregular intervals. Reports of long-distance transmission at 


frequencies as high as 56 megacycles, which occur mostly during the [ 


summer, appear to be by way of the sporadic £ layer. 
2. ABSORPTION 


In contrast to the sharply defined maximum usable frequencies 
shown by the present graphs, there are less well defined minimum 
usable frequencies for practical high-frequency sky-wave transmis- 
sion for various distances. These lower frequency limits are deter- 
mined by absorption. The absorption varies in general with time of 
day, season, frequency, and length of path. Absorption occur 
mainly in the lower ionosphere, that is, in the Z layer and below. ‘The 
absorption is greater Sung the summer day than during the winter 
day. It is greater for the lower high frequencies, that 1s, those fre- 
quencies which are close to or below the maximum usable frequency 
for the E layer for a given distance. It is usually greater during the 
day than during the night, especially for these lower high frequencies. 
A high-frequency wave at large angles of incidence, corresponding to 
long distances, is absorbed like a much lower frequency wave at 
angles of incidence, corresponding to short distances. 


? Samuel 8. Kirby and Elbert B. Judson, Recent studies of the ionosphere, J. Research NBS 14, 460 (1989) 
Proc. Inst. Radio Engrs. 23, 733 (1935). ; 
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3. SCATTERED REFLECTIONS 


tive 
at Usually, complex scattered reflections are observed at frequencies 
cm, | considerably higher than the F,-critical frequencies.’ These reflec- 
eC tions provide weak, ag ag transmissions much above the maxi- 
bd mum usable frequencies for F-layer transmissions over distances of 
' several hundred kilometers. 
nal 
ible 


VI. CONCLUSIONS 


8 For practical applications the principal value of these graphs is to 
yer — estimate transmission conditions in the future either diurnally, season- 
yer — ally, or over longer periods. The ionosphere data of reference 1, 
ons — from which these graphs were derived, indicate that the diurnal and 
of — seasonal characteristics are regular and may in general be predicted. 
ion — In addition, there is a large variation of ionization densities with the 
500 — il-year sunspot cycle, as indicated by the increase of maximum usable 
frequencies from 1933 to 1937. At the end of the year 1937 the sun- 
spot cycle was near a maximum and is expected to return to a mini- 
' mum about 1944. In a rough way these graphs may be used for 
corresponding times on the descending part of the sunspot cycle. 
on. | Current ionosphere data are given by the National Bureau of Stand- 
the | ®tds in its weekly ionosphere bulletins broadcast from station WWV 
yhe- | cach Wednesday, and in reports on the characteristics of the iono- 
‘are | sphere at Washington, D. C., published each month in the Proceedings 
mes | of the Institute of Radio Engineers. 


of 
by WasuineTon, February 11, 1938. 
n of ' Samuel 8. Kirby and Elbert B. Judson, Recent studies of the ionosphere. J. Research NBS 14, 469 (1935); 


Proc. Inst. Radio Engrs. 28, 733 (1935). 
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BOILING POINTS OF BENZENE, ETHYLENE CHLORIDE, 
n-HEPTANE, AND 2,2,4-TRIMETHYLPENTANE OVER THE 
RANGE 660- TO 860-MM PRESSURE 


By Edgar Reynolds Smith and Harry Matheson 


ABSTRACT 


By the comparative ebulliometric method, using water as the reference standard 
and ebulliometers of the Swietostawski type, data were obtained from which were 
developed the following equations expressing the relation between temperature 
and vapor pressure from 660 to 860 mm: 

For benzene 


t= 80.094 + 0.042683 (p — 760) — 0.00002199(p — 760)? +-0.0000000250(p — 760)?, 
and 

p= 760+ 23.429(t — 80.094) + 0.2804 (¢ — 80.094)? 0.000383 (t — 80.094); 
for ethylene chloride 

t= 83.483 + 0.042037 (p — 760) — 0.00002201 (p — 760)? +-0.0000000180(p» — 760)!, 
and 

p= 760+ 23.789(t — 83.483) +-0.2946 (t— 83.483)? +-0.00156(¢ — 83.483); 
for n-heptane 
and t= 98.413+ 0.044849 (p — 760) —0.00002344(p — 760)? + 0.0000000145(p» — 760)5, 
and 

p= 760+ 22.297 (t— 98.413) +0.2599(¢— 98.413)? + 0.00246 (t — 98.413)8; 
for 2,2,4-trimethylpentane (isooctane) 

t=99.234 + 0.046511 (p— 760) —0.00002374(p — 760)? +-0.0000000182(p—760)°, 
and 

p=760+ 21.500(t — 99.234) + 0.2348 (¢— 99.234)? + 0.00126 (t— 99.234)*. 


In these equations, ¢ is the temperature in degrees centigrade at which the sub- 
stance exerts the vapor pressure p in standard millimeters of mercury. The use 
of these substances as reference standards for boiling-point measurements is dis- 
cussed. Equations for the relative volatilities of the two pairs, benzene-ethylene 
chloride and n-heptane-2,2,4-trimethylpentane, are given. The data show that 
Dihring’s rule and the corresponding reciprocal temperature relationship hold 
accurately over a limited range of pressures, even for substances of different types. 


CONTENTS 
I. Introduction 


II, Apparatus and materials 
III. Experimental results and conclusions 


I. INTRODUCTION 


The comparative method of measuring vapor pressures of liquid 
substances requires a primary standard reference Sra for which the 
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relation between temperature and pressure is known with satisfacto 
precision. Undoubtedly, water is the best available liquid for 4 
primary standard, for reasons which have been discussed by Swieto. 
stawski.'? Recent measurements of the temperature-pressure relg. 
tions of water between 660 and 860 standard mm pressure ® are of 
adequate precision for a reference standard, and cover a most cop. 
venient and useful range of pressures for ebulliometric measurements, 
According to these measurements the boiling point of water between 
660- and 860-mm pressure is given by the equation 


t= 100+ 0.0368578(p— 760) —0.000020159 (p—760)?+ (1) 
0.00000001621(p—760)°, 


in which t is expressed in degrees centigrade (Int.) and 7: in millimeters 
of mercury (standard). This equation is in excellent agreement with 
the data of the National Bureau of Standards ‘* as well as with those 
obtained recently by Zmaczynski® at the Physikalisch-Technische 
Reichsanstalt, and is apparently reliable to better than 0.001° on the 
International Temperature Scale for nearly 4° above and below the 
fixed reference point of 100°. 

For comparative measurements of pressure, it is desirable to have 
an equation in which pressure is an explicit function of temperature, 
To obtain such an equation, the data of Beattie and Blaisdell were 
treated by the method of least squares to obtain the constants in the 
following equation 


p=760+27.1313(t—100)+0.40083(¢— 100)?-+-0.003192(t—100)®, (2) 


which is valid over the same range of pressures. Equation 2 repro- 
duces the 38 experimental points from which it was obtained with an 
average deviation of 0.012 mm and a greatest deviation of 0.031 mm, 
and agrees with the NBS data with equal precision over the range in 
which the two sets of data overlap. It may be well to mention that 
the equation fails rapidly above and below the range for which it was 
obtained.® 

In the experimental procedure, successive measurements were 
made of the boiling point of one of the substances under investigation 
and the boiling point of water, in ebulliometers connected to the 
same manostat, by means of which the pressure was varied in steps 
from 660 to 860 mm. The pressures were calculated from the meas- 
ured boiling points of water by eq 2 and tabulated with the correspond- 
ing boiling points of the substance. From the data tabulated for 
each substance, the constants were evaluated, by the method of least 
squares, for equations of the type 


t=t,+a(p—760)+6(p—760)?+¢(p—760)® (3) 


p=760+ ¢(t—ts) +r(t—tz)?+8(t—ts)’, (4) 


pan. é Swietostawski, Comptes Rendus de la XII Conferénce de I’Union Internationale de Chimie, p. 8 

1 W. Swietostawski and E. R. Smith, J. Research NBS 20, 549 (1938) RP1088. 

3J. A. Beattie and B. E. Blaisdell, . Am. Acad. Arts Sci. 71, 361 (1937). 

4N.8. Osborne and C. H. Meyers, J. Research NBS 13, 1 (1934) RP691. 

‘ Zmaczynski’s data over this range are given in the reference under footnote 2. 

* In the reference of footnote 2, there is given an alternative equation involving the logarithm of the 
sure, which is valid over a greater range. The additional alternative of interpolation from the of 
Osborne and Meyers, footnote 4, is also pointed out. 


and 














ory 
r a 
ala 
of 
on- 


as- 
\d- 
for 
ast 


(3) 
(4) 


), 89 








oe al Boiling Points of Some Organic Liquids 643 
in which a, b, ¢, g, r, and s are constants, ¢ is the temperature in 
degrees centigrade at which the vapor pressure is p, and fy is the 
normal boiling point, which is determined as one of the computed 
constants. 


II. APPARATUS AND MATERIALS 


For the determination of boiling points two ebulliometers were used. 
One of these was a simple barometric ebulliometer containing re- 
distilled water for the reference standard. The other was a standard 
differential ebulliometer containing the substance under investigation. 
These ebulliometers were of the Swietostawski type 7 and have been 
described in other papers from this ogee 2 They were connected 
through drying tubes to a manostat, which was an 18-liter carboy 
placed in an insulated wooden case to insure uniformity of temperature. 
The pressure in the system was increased by the addition of nitrogen 
from a cylinder connected through a stopcock to the manostat, and 
decreased by evacuating through the same stopcock. The changes in 
pressure were adjusted roughly according to the levels of mercury in 
an open-arm manometer, a so connected to the manostat through a 
stopcock, final adjustments bean made as indicated by observations 
of the boiling point of water. The procedure involved the approxi- 
mate adjustment of the pressure to a desired value, say about 660 mm, 
followed by three temperature measurements. The temperature of 
the water boiling in one ebulliometer was first measured. The ther- 
mometer was then transferred, at equal short time intervals, from the 
ebulliometer containing water to the one containing the substance 
under examination and then back to the water, a temperature measure- 
ment being made in each position. Averaging the temperature 
readings of the water corrects for any possible slight drift in pressure 
in the manostat during the operation. Thus when water boils at this 
average temperature, t», the substance boils at the measured tempera- 
ture, f;, under the same pressure. This pressure is calculated from 
eq 2, and thus there are obtained a series of temperatures and cor- 
responding vapor pressures for the substance relative to water, as the 
primary standard. 

A calibrated platinum resistance thermometer (coiled-filament 
type)® and a calibrated Mueller thermometer bridge were used to 
measure the temperatures. 

Water—For the ebulliometric reference standard, distilled water 
from the laboratory supply was distilled from alkaline permanganate 
and again distilled Bae. 2 before use. 

Benzene.—The initial material was benzene which conformed to 
the specifications adopted by the American Chemical Society for the 
reagent grade of this substance.!° Two liters of this material was 
distilled fractionally in a 40-bulb vacuum-jacketed column of the 
Swietostawski type" equivalent to about 25 theoretical plates. 
The distillate was collected in eight fractions, each of about 250 ml. 
The normal boiling points of the third and fifth fractions were meas- 


Pa Swigtostawski, Ebulliometry, p. 7 and 16. (Chemical Publishing Company of New York, Inc., 
* M. Wojciechowski, J. Research NBS 17, 453 (1936) RP921. 
E. R. Smith and M. Wojciechowski, J. Research NBS 17, 841 (1936) RP947. 
no H. Meyers, BS J. Research 9, 807 (1932) R P508. 
ind. Eng. Chem., Anal. Ed. 4, 347 a? 
W. Swietostawski, Ebulliometry, p. 27 (Chemical Publishing Company of New York, Inc., 1937). 
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ured by the comparative method and found to be 80.093 and 80.094° 
in agreement with the value of 80.094° reported by Wojciechowski # 
for benzene of high purity. The difference between the boiling point 
and the temperature of condensation of these fractions, measured jn 
the differential ebulliometer, was 0.001°, which also indicated high 
purity. Since these fractions were of satisfactory purity, they were 
used for the measurements and the other fractions were not tested. 

Ethylene chloride—About 2 kg of this substance, obtained from 
Eastman Kodak Co., was distilled in the 40-bulb column and ¢ol. 
lected in a series of fractions, each of about 150 ml. Fractions 5, 7, 
and 10 had normal boiling points of 83.478, 83.481, and 83.482° (, 
respectively, according to single measurements, and each showed in 
the differential ebulliometer a difference of 0.002° between boiling 
point and temperature of condensation. Fraction 7 was used for the 
measurements. 

n-Heptane.—The initial material obtained from the Eastman 
Kodak Co. showed in the differential ebulliometer a difference between 
boiling-point and condensation temperature of 0.04°. After frac. 
tionation in the 40-bulb column, this difference was reduced to 0.001° 
for fraction 3, which was used for the measurements. A prelimin 
measurement of its boiling point gave 98.413° C. Since this result 
is slightly higher than the value of 98.365° obtained by Wojcie- 
chowski,“ and lower than the value of 98.52° found by Beattie and 
Kay,'* a sample of synthetic n-heptane made by the California 
Chemical Co., was obtained through the kindness of F. W. Rose, Jr, 
Research Associate of Project 6 of the American Petroleum Institute. 
This sample had been prepared by fractionation of the initial synthetic 
material through a 2.5-m locket-chain column, equivalent to about 
60 theoretical plates. The difference between its boiling point and 
condensation temperature was 0.001°, and its boiling point in the 
differential ebulliometer was 98.413°, in complete agreement with 
the sample used for this work. 

2,2,4-Trimethylpentane (isooctane).—The sample used for measure- 
ment was a portion of material made by Réhm and Haas, which was 
further purified by successive fractional distillation. The sample was 
kindly supplied by D. B. Brooks, of the automotive power plant sec- 
tion of this Bureau, who will describe its preparation and properties in 
detail in a subsequent paper. This sample had a freezing point of 
—107.31° C, and a difference of 0.001° between boiling point and con- 
densation temperature in the differential ebulliometer. 


III. EXPERIMENTAL RESULTS AND CONCLUSIONS 


Benzene.—The comparative boiling point data for benzene and 
water are given in table 1. The boiling point of benzene is shown In 
column 1, and the corresponding boiling point of water in column 2. 
The corresponding pressure, calculated from the boiling point of water 
by eq 2, is given in column 3. 

“ut Wolstechowskt, J. Meeearch N38 47, 408 (800) HVA. ‘The material used in Woicechowi 


between bo! point and temperature of condensation. 
“J, A. Beattie and W. C. Kay, J. Am. Chem. Soc. 69, 1586 (1937). 
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TaBLE 1.—Boiling points of benzene between 660- and 860-mm pressure 
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The expression obtained for the boiling point of benzene as a func- 
tion of pressure 1s 


t=80.094 +-0.042683 (p —760) —0.00002199 (p—760)?-++-0.0000000250 
(p—760)? (5) 


The average deviation of the 15 measurements from this equation is 
0.001° and the greatest deviation is 0.003°. The normal boiling point 
is 80.094°, in agreement with a previous determination in this lab- 
oratory by Wojciechowski.” It is interesting that the value obtained 
by Wojciechowski is an average for several samples of benzene purified 
by different methods and measured at pressures within a few milli- 
meters of normal, while in this work the boiling point is obtained as a 
constant to fit the equation according to measurements on one sample 
from 100 mm below to 100 mm above normal pressure. The expres- 
sion obtained for pressure as a function of temperature is 


p=760+23.429 (¢(—80.094) +0.2804 (¢—80.094)?—0.000383 (6) 
(¢—80.094)*, 


which reproduces the measurements with an average deviation of 
0.03 mm and a greatest deviation of 0.05 mm. General expressions 
— values of dt/dp and dp/dt may be obtained easily from eq 5 
and 6. 

Ethylene chloride——The data for ethylene chloride are given in 
table 2. The equations for the temperature-pressure relations of 
this substance are 


t=83.483 +0.042037 (p—760)—0.00002201 (p—760)?-+0.0000000180 
(p—760) (7 


say an average deviation of 0.001° and a greatest deviation of 0.002°, 
an 
p=760+-23.789 (¢— 83.483) +0.2946 (¢— 83.483)?+-0.00156 
(t—83.483)* (8) 


with an average deviation of 0.01 mm and a greatest deviation of 
0.04 mm from the experimental values. 


LS 
“ M. Wojciechowski, J. Research NBS 19, 347 (1937) RP 1029. 
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TaBLE 2.—Boiling points of ethylene chloride between 660- and 860-mm pressuy 
p 
Boiling point 
A Pressure 1 
thylene 
chloride Water 

°C °C mm 
79. 325 96. 129 660. 80 
79. 628 96. 613 672. 57 
80. 100 97.029 682. 84 
80. 644 97. 509 604. 85 
81. 209 98. 004 707.42 
81. 556 98. 309 715. 24 
82. 273 98. 939 731. 66 
82. 898 99. 487 746. 19 
83. 487 100. 004 760. 11 
83. 494 100. 011 760. 30 
84. 065 100. 511 773. 97 
84. 638 101. 012 787. 86 
85. 198 101. 504 801. 72 
85. 782 102. 014 816. 29 
86. 358 102. 518 830. 90 
86. 912 103. 000 845. 10 
87. 486 103. 502 860. 06 

















n-Heptane.—The equations for n-heptane, based on the data given 
in table 3, are 


TaBLeE 3.—Boiling points of n-heptane between 660- and 860-mm pressure 

















Boiling point 
Pressure 
n-Heptane Water 
°C A, *, mm 
93. 432 95. 893 655. 11 
94. 046 96. 400 667. 38 
94. 675 96. 920 680. 15 
95. 341 97. 468 693. 82 
95. 650 97.725 700. 31 
96. 005 98. 018 707. 83 
96. 634 98. 536 721. 13 
96. 946 98.795 727.88 
97. 251 99. 046 734. 48 
97.748 99. 456 745. 36 
99. 020 100. 500 773. 67 
99. 672 101. 084 788. 49 
100. 480 101. 696 807.18 
100. 870 102. 016 816. 35 
101. 504 102. 537 831. 46 
102. 121 108. 045 846. 42 
102. 740 108. 551 861. 54 











t=98.413-+0.044849 (p—760)—0.00002344 (p— 760)? ++0.0000000145 
(p—760)8 (9) 


which reproduces the 17 experimental values with an average deviation 
of 0.001° and a greatest deviation of 0.003°, and 


p=760+22.297 (¢(—98.413)+0.2599 (¢—98.413)?+0.00246 — 





with an average deviation of 0,02 mm and a greatest deviation of 
0.07 mm. 
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9,2,4-Trimethylpentane (isooctane).—The data for 2,2,4-trimethyl- 
pentane are given in table 4. 
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Taste 4.—Boiling points of 2,2,4-trimethylpentane between 660- and 860-mm 
pressure 














Boiling points 
Pressure 
2,2,4- Tri- Water 
metnyip ane; 
Al”. °C mm 

94, 318 96. 090 659. 86 
94. 949 96. 592 672. 08 
95. 471 97.007 682. 31 
96. 088 97. 501 694. 64 
96. 706 97.992 707.10 

97.111 98. 318 715. 47 - 
97. 860 98.912 730. 94 
98. 598 99. 496 746. 44 
99.310 100. 061 761.65 
99. 869 100. 503 773.74 
100. 501 101. 003 787.61 
101. 117 101. 491 801. 36 
; 101. 768 102. 003 815. 99 
: 102. 404 102. 506 830. 57 
ven 108. 005 102. 980 844. 50 
108. 642 108. 481 859. 45 

















The equations obtained from the data are 


t=99.234+-0.046511(p—760) —0.00002374 (p—760)? 
+-0.0000000182(p—760)° (11) 


with average and greatest deviations of 0.001 and 0.003°, respectively, 
and 
p=760+21.500(t—99.234) +-0.2348 (t— 99.234)’ 
+0.00126(t— 99.234) (12) 


with average and greatest deviations of 0.02 and 0.07 mm, respectively. 

In table 5 are listed the normal boiling points and the values of 
(dt/dp)r¢0 for the four substances studied. It is noteworthy that these 
values have not been obtained from single measurements in the 
vicinity of normal pressure, but have been calculated as constants in 
the equations which reproduce 15 to 17 experimental values over the 
range 660- to 860-mm pressure. Of these four substances, benzene is 
the only one which would be suitable for a secondary ebulliometric 
reference standard. On prolonged boiling ethylene chloride undergoes 
a slight hydrolysis if a trace of water is present, while n-heptane and 
2,2,4-trimethylpentane have boiling points within the range where a 

45 secondary standard is superfluous. 


TaBLE 5.—Normal boiling points and values of dt/dp 
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Relative volatility.—For calculating the number of theoretical plate 
required in a fractionating column to obtain a given separation of gp 
ideal two-component mixture, the equation 


%o(1—2,) 
In(1—2) (13) 
og & 


is often useful. In this equation, n is the number of plates, R is the 
ratio, assumed constant, of the vapor pressures of the two pure cop. 
stituents, and 2) and z, are the mole fractions corresponding to n=( 
and n=n, 2 being greater than z,. Beatty and Calingaert have 
shown that the mixture n-heptane-2,2,4-trimethylpentane is very 
nearly ideal from the standpoint of vapor pressure, and also that in 
any case where £ is nearly unity, a change in the value of RF creates 
marked difference in the palouldted value of n, but that in such g 
case the boiling range is small, so that the absolute change in R is algo 
small over the boiling range. The effect on n of a change in R is thus 
considered negligible. The values of R at the normal boiling points 
of n-heptane (98.413°) and 2,2,4-trimethylpentane (99.234°) can be 
calculated accurately by eq 10 and 12 to be 1.02356 and 1.02432, 
respectively, for this pair. The corresponding values of n, from 
eq 13, for a change in mole fraction from 0.05 to 0.95 are 252.9 and 
245.1, respectively, the difference over this range of temperature 
being negligible for practical purposes. Values of R for the systems 
n-heptane-2,2,4-trimethylpentane and benzene-ethylene chloride over 
a range of temperature are given in table 6. For these two pairs the 
relative volatilities are expressed by the following equations: 


log 


=> 


TABLE 6.—Values of R 





n-Heptane-2,2,4-trimethylpentane Benzene-ethylene chloride 





R R 
Tempera- p(n-heptane) . p(benzene) 
ture | 3@,2,4-trimethylpentane) p(ethylene chloride) 











1. 01995 
1, 02089 
1, 02181 

02272 


1. 02364 
1. 02456 
1. 02550 


1. 02647 
1. 02746 


BESS BSESES 
Aang coro orc 


=e 























For n-heptane and 2,2,4-trimethylpentane between 94.5 and 102.5° C 
R=1.02088—0.0008783t-+-0.0000091982?, 


where R=p(n-heptane)/p(2,2,4-trimethylpentane). For benzene and 
ethylene chloride between 79.0 and 84.0° C 


R=1.14596 +0.0004845t—0.000011152#, 


where R=p(benzene)/p(ethylene chloride). 


16H. A. Beatty and G. Calingaert, Ind. Eng. Chem. 26, 504 (1934). 
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Dihring’s rule-—An empirical relationship between the tempera- 
tures at which two liquids have equal vapor pressures, known as 
Dihring’s rule, has found considerable use for the interpolation of 
data. he relationship is expressed by the equation 


T,=kh,T;+(C,, (14) 


where k, and C, are constants, and 7’, and 7; are the temperatures at 
which the liquids A and B have the same vapor pressure. A similar 
relationship with a better claim to some theoretical basis involves the 
reciprocal temperatures 


1/T'4= (k2/T3) +C, (15) 


Although these two relations have been tested before,” the relative 
vapor pressure data of this investigation afford a more precise com- 
parison than has hitherto been made, over the short range of 200 mm, 
of the relative validities of eq 14 and 15. 

In preference to an algebraic elimination between equations of type 
3, the temperatures at which water, benzene, ethylene chloride, 
n-heptane, and isooctane exert the 11 vapor pressures of 660, 680, 700, 
_, . . 840, 860 mm were calculated from eq 1, 5,7,9,and11. From 
these 11 sets of values the constants of eq 14 and 15 were found for 
each of the 10 possible pairs of the 5 substances. The accuracy with 
which these linear equations of types 14 and 15 hold is illustrated by 
the example of benzene and water given in table 7. 


TaBLE 7.—Comparison of temperatures calculated from eq 12 and 13 for the pair 
benzene-water 





Temperature of— Tempera- 


At ture of At 
Benzene,| Water, Water, | ©21—ea 16 water, | eq 1—eq 17 
eq 5 eql eq 16 eq 17 








°C “> 
75. 581 96. 096 
76. 526 96. 914 
77. 448 97.712 
78. 350 98. 492 


79. 231 99. 255 
80. 004 100. 000 
80. 939 100. 729 
81. 767 101. 443 


82. 581 
83. 380 
84. 167 


°C 
96. 098 
96. 916 
97.714 
98. 494 


99, 255 
99. 999 
100. 728 


© ide 
se | 888 S222 S8ee° 
lige 
see 


| 
8 





ESS -OtFh: 33! 
sae P27. 
ess g8e8 


ane 











Average deviation 
Greatest deviation 


28 




















The equations for this pair are 
t»=0.86268t,+30 900, (16) 
in which ¢ is expressed in degrees centigrade, and 


ir = (0.772945/T,) +0.00049175, (17) 


cialis nig 
" For example, see E. F, Linhorst. J. Am. Chem. Soc. 61, 3262 (1929); A. M. White. Ind. Eng. Chem. 
%, 230 (1930); J. H. Perry and E. R. Smith. Ind. Eng. Chem. 25, 195 (1983). % 


57763 —38—_——6 
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in which T=t+273.160. Without tabulating the results of the tesi; 
with the other nine pairs, an idea of the relative accuracy of eq 14 and 
15 can be gained from the statement that for all pairs the average of 
the average and greatest deviations was found to be 0.002 and 0,004 
from eq 14, and 0.001 and 0.003° from eq 15, respectively. Thus, in 
regard to accuracy there is very little choice between the two, tho 
the reciprocal temperature relation seems to be slightly better, on 
average, for the 10 pairs that were tested. 


WASHINGTON, January 18, 1938. 
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SUITABILITY OF VARIOUS PLASTICS FOR USE IN 
AIRPLANE DOPES 


By Gordon M. Kline and Cyrus G. Malmberg 


ABSTRACT 














A number of different cellulose derivatives of potential interest as film-forming 
constituents for airplane dopes have recently become available commercially. 
An investigation has been undertaken at the National Bureau of Standards to 
determine the fundamental factors involved in the formulation of dopes con- 
taining these new derivatives to obtain optimum tautness and durability. Data 
are presented in this paper relative to the effect of varying the acyl or ethoxyl 
content and the viscosity of cellulose esters and ethers on the tautness of fabrics 
doped with them. It is concluded that the solvents and diluents govern to a 
large extent the tautening properties of the dope and the durability of the film 
deposited on the fabric. 
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I. INTRODUCTION 


The problem of producing a fire-resistant doped fabric to replace 
the hazardous cellulose nitrate product was investigated in 1933 at 
the National Bureau of Standards for the National Advisory Com- 
mittee for Aeronautics, and the results were reported in this journal.’ 
It was established that none of the synthetic resins then available 
would produce films of satisfactory tautness when used alone on the 
fabric, and that their addition to either cellulose nitrate or cellulose 
acetate dope resulted in a corresponding lowering of the fabric tension. 
Data were also obtained on the comparative rates of burning of 
fabrics doped with various cellulose nitrate and cellulose acetate com- 
positions, both with and without the addition of fire-retarding salts 
to the fabric. It was shown that cellulose acetate dope yields a 
relatively nonhazardous product when applied to untreated fabric, 


'G. M. Kline, Fire-resistant doped fabric for aircraft, J. Research NBS 14, 575-587 (1935) RP788. 
651 
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and that when it is applied to fabric treated with boric-acid-borg 
mixture, a product which is nonflammable under ordinary condition; 
is obtained. However, the investigation was not continued at tha 


time to include a study of the effect of varying the ingredients of, J 
cellulose acetate dope on the tautening characteristics and particy. f 


larly on the moisture sensitivity of the doped fabric. 
Cellulose acetate has now become available commercially in sever) 


grades, representing materials of different degrees of acetylation and 


viscosity. Likewise, other cellulose derivatives, possessing character. 


istics which indicate that they might make satisfactory dope base. § 
ese more recently developed prod. 

ucts are ethylcellulose, cellulose acetobutyrate, cellulose acetopm. f 

pionate, and cellulose nitroacetate. These products are also availabk § 


are being manufactured. Among t 


in various grades denoting differences in chemical composition and 
in viscosity. The proper formulation of dopes based on these ney 
materials to obtain the requisite balance of solvents, diluents, and 


—— both as to types and amounts, can be accomplished only F 


y a thorough laboratory study of the various factors involved. The 


Bureau of Aeronautics of the U. S. Navy Department has, therefore, F 
requested the National Bureau of Standards to undertake such ap | 
experimental study in order to develop a dope, based on these com- 


paratively nonflammable cellulose derivatives, which will compar 

favorably with or surpass cellulose nitrate dope with respect to the 

effect of high relative humidity on the tautness of the doped fabric. 
In the first phase of this work a variety of cellulose derivatives and 


synthetic resins were —_— to airplane fabric to study the relation 
plastic base used in the dope, the percentage fF 


of tautness to the type of plas 
of acyl or alkyl substitution in the various cellulose derivatives, and 


the viscosity * of these derivatives. These panels are also ps 
oped fab- F 


data on the relative resistance to weathering of the various 
rics, while tests are in progress to determine the effect of solvents, 


diluents, and plasticizers on the tautening property and to establish [ 


the moisture relations of the various doped fabrics. 


II. DETAILS OF TEST MATERIALS AND MEASUREMENTS | 


1. PROPERTIES OF THE PLASTICS 


The materials tested were furnished by various manufacturers 
whose interest and cooperation in this investigation have been very 


helpful. They also supplied data regarding the viscosity and chem- 


cal composition of the compounds, shown in table 1. In the manv- 
facturers’ reports the chemical composition was expressed as percent- 
ages. The number of equivalents of each substitution group present 
in the compounds and also of free hydroxyl groups were caleula 
from these percentage figures. These latter values are listed in table! 
and indicate immediately the degree of hydrolysis, which is closely 
related to the solubility’ and moisture-absorbing‘ properties of the 
cellulose derivative. 


1 The term “viscosity” is used in this paper to denote certain flow characteristics of solutions of the vari- 
on gineiien & meerenet. Oa meena Soe in the footnote to table 1. ses T-37eT 
3 W. Coltof, Solubility properties of certain highly polymeric substances. J. Soc. Chem. Ind. 56, 


1 \ 
: TOE. Sheppard and P. T. Newsome, The sorption of water vapor by cellulose and its derivatives. J. Phys. 
Chem. $3, 1817-35 (1929). 
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TABLE 1.—Viscosity and composition of test materials 





Sample des- 
ignation 


Viscosity * 


Composition 





Percentage basis 


Equivalents basis 





CA-al.... 


CA-a5.... 
CA-a6_... 
CA-al0..- 
CA-b1._.. 
CA-b2__-- 
CA-b3_... 
CA-b4_... 
CAP-al._. 
.| CAP-bl1.. 
CAP-b2.. 


CAP-b3.. 


CAP-cl_.- 
CAP-c2... 
-| CAP-c3... 


CAB-b1.. 
CAB-b2.. 
CAB-b3.. 


Sec... 

34.5 cp &__ 
470 
650 


ep... 
cp... 


35 


cp... 
cp... 


SS8 SS55S 
-a Ne OF DOF 





17.1 acetyl; 26.9 propionyl_ 
15.5 acetyl; 31.3 propionyl. 
15.5 acetyl]; 31.3 propionyl. 

15.5 acetyl; 31.3 propionyl- 


14.0 acetyl; 34.0 propionyl_ 
11.0 acetyl; 34.0 propionyl_ 
9.0 acetyl; 33.0 propionyl. 


.| 32.5 acetyl; 15.4 butyryl__. 


32.2 acetyl]; 15.2 butyryl__. 
32.0 acetyl; 16.1 butyryl_.. 


2.38 acetyl; 0.62 hydroxyl. 
2.48 acetyl; 0.52 hydroxyl. 
2.57 acetyl; 0.43 hydroxyl. 
2.51 acetyl; 0.49 hydroxyl. 
2.48 acetyl; 0.52 hydroxyl. 


2.32 acetyl; 0.68 hydroxyl. 
2.35 acetyl; 0.65 hydroxyl. 
2.32 acetyl; 0.68 hydroxyl. 
2.38 acetyl; 0.62 hydroxyl. 
ca.3 acetyl. 


2.37 acetyl; 0.63 hydroxyl. 
2.51 acetyl; 0.49 hydroxyl. 
2.39 acetyl]; 0.61 hydroxyl. 
2.29 acetyl; 0.71 hydroxyl. 


1.13 acetyl]; 1.34 propiony]; 
0.53 hydro xyl. 

1.08 acetyl; 1.64 propiony]; 
0.28 hydroxyl. 

1.08 acetyl]; 1.64 propiony]; 
0.28 hydroxyl. 

1.08 acetyl; 1.64 propiony]; 
0.28 bydroxyl. 


1.00 acetyl; 1.83 propiony]; 
0.17 hydroxy]. 

0.74 acetyl; 1. 72 propiony]; 
0.54 hydroxyl. 

0.58 acetyl; 1.60 propiony], 
0.82 hydroxyl. 


2.31 acetyl; 0.66 butyryl; 
0.03 hydrox yl. 

2.27 acetyl; 0. 08 butyryl; 
0.08 hydroxy 

2.28 acetyl; 0. butyryl; 
0.03 hydroxyl. 


Cellulose nitroacetate.--_. Medium... 


Cellulose nitrate 
Do 


Methylcellulose 


aon 47.5 ethoxy] 
47.7 ethoxy] 
48.2 ethoxy! 
48. 3 ethoxy] 


.| 52.0 ethoxyl 


.| 48.7 ethoxyl 
48.0 ethoxyl 
48.3 ethoxyl 
49.5 ethoxy! 


2.43 ethoxy]; 0.57 hydroxyl. 
2.44 ethoxy]; 0.56 hydroxyl. 
2.48 ethoxy]; 0.52 hydroxyl. 
2.49ethoxy]; 0.51 prs 
2.77 ethoxy]; 0.23 hydroxyl. 


2.52 ethoxy]; 0.48 hydroxyl. 
2.47 ethoxy]; 0.53 hydroxyl. 
2.49 ethoxy]; 0.51 hydroxyl. 
2.58 ethoxy]; 0.42 hydroxyl. 


EC-d2._.. 
EC-d3._.. 


EC-d6_... 
BC-b1.... 
CR-b1.-... 


Methyl methacrylate_...| MM-cl__- 

















*The viscosity values were obtained ~ the manufacturers’ laboratories by the following variety of methods 


(all parts by = unless otherwise noted): 

Cellulose te series CA-a: ASTM falling-ball method, using 20-percent solutions in acetone. (See 
ASTM Standards. 1936, part II, Nonmetallic Materials, published by the American Society for Testing 
Materials, yee ny ae ) 

Cellulose acetate series CA-b: ASTM falling-ball method, using 20-percent solutions in a mixture of 90 
parts of acetone and 10 parts of ott 1 alcohol. 

Cellulose acetopropionate series CAP-a: ASTM a method, using 20-percent solutions in acetone. 

sem acetopropionate series CAP-b: capillary meter, using 10-percent solution in acetone at 


ei series re! Same as poor, 


ders using a 20-percent solution in a mixture of 80 
viscometer method, using a ae solution in a mixture of 80 
ute ethyl alcohol by volume at 25' 
viscometer, using 8 hd yah in toluene at 25° O. 
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2. COMPOSITION OF THE DOPES 


Inasmuch as it was desired to compare the behavior of various gg. 
lulose derivatives with one another, the same solvent mixture anj 
plasticizer, triphenyl phosphate, were employed in a majority of th 
dopes. In most cases the proportions of the latter were one part of 

lasticizer to nine parts of cellulose derivative. Preliminary test, 
indicated that varying the concentration of the cellulose base in the 
solvent mixture between the limits of 5 and 10 percent had compan. 
tively little effect on the tautness produced. A concentration of 
about 6.4 percent was, therefore, used for most samples, although it 
was necessary to dilute some of the aluminum-pigmented dopes ty 
obtain a solution of suitable viscosity for brushing. Complete day 
concerning the formulas of the dopes are shown in table 2. Formula] 
was used in general for the cellulose esters and formula VI for th 
ethylcelluloses. The other formulas were introduced because of the 
special solubility characteristics of certain derivatives, such as celly. 
lose triacetate, or to compare with those above solvent compositions 
recommended in the literature or by the manufacturer. 


TABLE 2.—Formulas of experimental dopes 





Alu- 


Plas- i 
> minum} Sol- | Composition of solvent (parts by 
Formula tic Plasticizer pow- | vent weight) ' 


base der 





Per- Per- 
cent Percent cent ; 
6.4 | 0.6 triphenyl phosphate_| None Acetone 48, ethyl] acetate 20, methy! F 
ethyl ketone 20, diacetone alcohol 


5. 
6.0 | 1.2 triphenyl phosphate.| 4.8 ; Acetone 50, ethyl acetate 15, methy! 
ay ketone 15, diacetone alcohol 


IA, clear.....-------] 6.4 3. Acetone 48, ethy! acetate 20, methy! 
eee? ketone 20, diacetone alcohol 


IA, pigmented--.._- q ' \ Acetone 50, ethyl! acetate 15, methy! 

aad ketone 15, diacetone alcohol 
II, clear . 0.6 triphenyl phosphate_ 
II, «Tee ave me 1.2 tripheny! phosphate_ 


0.3 > ~ snd phosphate 
0. 


Acetone 50, ethy! alcohol 10, toluene 
40. 


IIB, clear.....----.- 
IIB, pigmented 
III, clear 
III, pigmented 

sh tS 
IV, pigment 


WERE ent n sss 


Acetone 16.5, ethy! alcohol 8.5, ben- 
is 36.6, toluene 18.0, xylene 9. 


Do. 
Do. 
Acetone 25, ethy! acetate 20, toluene 
oo alcohol 10. 
Do. 
Do. 


Ethyl acetate 25, moan thy! 
ketone 25, buty] acetate 25, Methy! 
Cellosolve 25. 


do 
0.6 Dow No. 6...-.------ 
1.2 Dow No. 6.......-..- 


Vik, BEINN wtpie non inne 4 | 0.6 triphenyl phosphate- 


8 BRS F SBR RF SRRSR S 


tripheny] phosphate. 


VI, pigmented 1.2 
None 


VILA, clear 
Vv P' 


Do. 
’ Do. 
IIA, pigmented _.-- 


Do. 

0.6 triphenyl phosphate- Methyl Cellosolve. 

1.2 triphenyl phosphate_ Do. 

0.4 triphenyl phosphate. Chloroform 70, tetrachlorosthane 
0. 
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TABLE 2.—Formulas of experimental dopes—Continued 





Plas- Sol- | Composition of solvent (parts by 
i. Plasticizer - | vent weight) 





ag 


Toluene. 
D 


0. 

Toluene 33.1, xylene 22.3, Cellosolve 
22.3, methyl amy] ketone 22.3. 

Toluene 26.8, xylene 24.4, Cellosolve 

was methyl amy] ketone 24.4. 
ater. 


Do. 


Fe 8 SSS 
onan #B Pww 


© 
— 























3. METHOD OF PREPARING THE DOPED PANELS 


Wooden frames (fig. 1), constructed in accordance with specifica- 


"tion SP-16 of the U. S. Naval Aircraft Factory, were used in these 


' tests. Cotton airplane fabric, weighing approximately 4 oz/yd’, 
_ was placed over the frame and a 50-pound weight was attached to each 
| side. In this way uniform tension in both directions of the weave of 

the cloth was produced. The fabric was tacked to the frame with 


FiacurE 1.—Construction of exposure panel. 


No. 4 upholsterer’s tacks, spacing them approximately 1 inch apart. 
Four coats of clear dope and, in the case of the pigmented panels, two 
coats of aluminum-pigmented dope were brushed on the panels. The 
nto weight of the clear doped fabrics after drying was 7.2 0z/yd* 
and of the pigmented doped fabrics was 8.7 0z/yd?. 


4. METHODS OF TESTING THE FABRICS 


The measurements of tautness of the fabric were made with a 
spring-loaded tautness meter, the construction of which is shown in 
gure 2. The instrument is essentially a Schiefer Compressometer,* 
modified to vy it to measuring the deflection of a doped fabric 


instead of the thickness and compressibility of textile fabrics and 
similar materials. 


‘H. F. Schiefer, BSJ. Research 10, 705-13, (1933) RPS561. The assistance of Dr. Schiefer in 
the tautness meter and in preparing figure 2 is gracefully’ ackuowledent, — 
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For measuring tautness, the instrument is er upon the fabric 
A, of a panel or wing section with the foot, C, midway between thy 
two ribs, B. The foot, C, is of spherical shape with a 0.5-inch radix 
of curvature. It is fastened to the bottom of the spindle, D, of th, 
lower dial micrometer, E. The rod, F, is fastened to the top of th 
spindle, D, at G and to the top of a helical spring, H, at J, Thy 
bottom of the spring is fastened to the tube, J, at K. The upper diy 
micrometer, L, is fastened to the top of the tube at M. The spindle 
N, of the upper dial micrometer is attached to the rod, F, at O bys 
ball and socket union. The tube may be moved up or down relatiy, 
to the frame, P, by turning the knob, Q, of the rack and pinion, R f 
The frame, P, consists of an aluminum beam having a three-poin, 
support which rests on the framework of the panel or on any ty 





FigurRE 2.—Diagram of spring-loaded tautness meter. 


adjacent ribs of a wing. Graduations on the beam indicate the 
width of the rib spacing so that the instrument may be centered 
quickly and conveniently. By turning the knob the foot may b | 
lowered upon the doped fabric. The load which is applied to the 
specimen by the foot may be ascertained from the upper dial reading, F 
which indicates the elongation of the spring, and a calibration curve f 
of the spring. The amount of deflection of the doped fabric whe fF 
loaded is indicated on the lower dial, which is graduated to read 
directly to 0.001 inch. : 

For the estimation of tautness of doped fabrics, the foot, is lowered 
upon the specimen by means of the rack and pinion until 
on the fabric is 3 ounces, as indicated by the upper dial. A reading 
of the lower dial is taken at this load. Then the load on the fabrics 
increased to 19 ounces and a second reading of the lower dial is t 
The difference between these two readings of the lower dial is the 
deflection in thousandths of an inch of the fabric under one pound load, 
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Malnbers| 


as tested. The less the deflection under a given load, the greater is 


the tautness, ; i 
Tautness measurements were made on the fabrics conditioned at 


the |) approximately 70° F. and 65-percent relative humidity. For the 
the cloth alone a 2-ounce loading is employed instead of 1 pound in order 
The | to avoid excessive stretching. The average deflection observed for 
‘di |) the undoped fabrics with the 2-ounce load was 45 mils with an average 
nde | variation of +2.5 and a maximum difference of +6 mils. A 1-pound 

by load deflects such fabrics approximately 260 mils. 
ative After drying for at least 1 week in the laboratory, the doped panels 
1, 2 | were placed in a room kept at 70° F and 65-percent relative humidity. 
int | Three measurements of tautness were made in this room at intervals 
‘tw | during a period of 1 week. The — were placed upon the roof 
of the Industrial Building of the National Bureau of Standards on 


| August 14, 1937, on racks inclined at an angle of 45° and facing south. 
' Measurements of tautness and brittleness were made at intervals 
under prevailing conditions of temperature and humidity. Brittle- 
“ess of the film, i. e., tendency to “ringworm”, was determined by 
pressing thumb firmly into the fabric covering at each corner of 
» the panel. 


5. RESULTS OF TESTS ON THE DOPED PANELS 


Data obtained on the doped panels in the conditioning room and 
on the roof are presented in table 3. 
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III. EFFECT OF VARIOUS FACTORS ON INITIAL TAUTNEgs 
OF DOPED FABRICS 


1. TYPE OF PLASTIC : 
The tautness values of the fabrics doped with the different plastic 


materials, measured in the conditioning room at 70° F and 65-percen; 


relative humidity, are summarized in table 4. It will be noted tha 


most of the cellulose plastics produced tautnesses ranging between §) [ 


TABLE 4.—Summary of tautening properties of various types of plastics used in dopes 
containing aluminum pigment 





Number oo 
of tion 
under 
1-lb. load 


Material samples 
tested 





Mils 
Cellulose acetate (hydrolyzed) 
Cellulose triaceta te. 

Cellulose a 
Cellulose acetobutyrate 


S888 I2SS 


Ethylcellulose. 


8 


Benzylcellulose 
EE SESSA Mrem ee He 
Methyl methacrylate resin__- 











tO tT 


33 








* One panel read 76. 


and 97. Cellulose triacetate and the cellulose acetobutyrate, whieh 
was also practically a triester, had average tautness values of 60 and fF 
77, respectively. Great difficulty was experienced in obtaining a taut | 
fabric with benzylcellulose. Plasticizer was omitted from most of the 
ethylcellulose dopes because of the comparatively poor tautening qual- | 
ities of this type of cellulose derivative and its reported characteristic 
of forming flexible films without it. However, it was necessary to 
add plasticizer to the high ethoxyl compound (panels 74 and 75) 
because the unplasticized films of this material ringwormed spor- 
taneously 1 day after doping. The resinous materials, such as chlor § 
inated rubber and methyl] methacrylate polymer, did not tighten the 
fabric as well as did most of the cellulosic materials. Solutions of 
butyl methacrylate and Methyl Cellosolve methacrylate polymes 
failed to tighten the fabric when applied in the same solvent mixture 

used for the methyl derivative. 


2. ACYL OR ALKYL SUBSTITUTION 


The data in table 3 indicate that varying the acety] content od 
cellulose acetate between the limits of 2.29 and 2.57 equivalents did 
not have an apparent effect on the ability of the compound to tighten 
the fabric. Likewise, the results obtained with ethylcellulose cot 
taining from 2.43 to 2.77 ethoxyl equivalents did not show a 00 
tion between ethoxy] content and tautening property. There ww 
not sufficient variation in the acyl composition of the cellulose aceto- 
propionates and cellulose acetobutyrates to indicate its effect 
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| tautness. ‘Tests to determine this relationship are now under way 
- on materials of this type having a greater variation in acyl content. 


3. VISCOSITY 


| The absence of any effect of the size of the cellulose molecules, as 
| indicated by the viscosity of solutions of these materials, on the 
_ tautening property is shown by the data in table 3. Cellulose acetates 
| varying in viscosity from 2 to 94 seconds, as measured by the falling- 

ball test, did not differ markedly in their tautening ability. The same 
| was true for ethylcelluloses whose 5-percent solutions in a mixture of 
- 80 parts of toluene and 20 parts of ethyl alcohol gave viscosity values 
' of 20 to 210 centipoises. 


4. SOLVENT COMPOSITION 


' The results obtained with the different solvent formulas which 
were tested in this investigation lead to the rather startling conclusion 
' that the most important single factor involved in the tautening prop- 
' erty of a dope is the solvent composition employed. In the term 
' “solvent composition” are included those organic liquids which form 
' colloidal solutions of the cellulose derivatives, those which have only 
cosolvent and swelling action, and those which are merely diluents in 


“T the blend. It has been the usual practice to formulate the solvent 


» composition on the basis of such requirements as drying time, anti- 
| blushing characteristics, solvent property, and economy, without an 
| adequate appreciation that the choice of the solvent composition 
governs not only the tautening property of the dope but also the 
| flexibility (conversely the brittleness) of the film deposited on the 
fabric. 

| Tassie 5.—Effect of various solvent compositions on the production of tautness by 
cellulose dopes containing aluminum pigment 





Cellulose | Cellulose | Cellulose | Genuiose | Ethylecel- 
acetate | 























| Table 5 presents the data obtained on the effect of various solvent 
| compositions on the tautness of the doped fabric. It will be noted 
that the difference in the tautness values of the two cellulose nitrate 
panels listed in this table is much greater than the tautness differences 
shown in table 4 for the various types of cellulose derivatives or for 
compounds varying only in viscosity and percentage substitution. 
Since these panels were prepared, the effect of various solvents and 
mixtures of solvents and nonsolvents on the shrinkage and flexibility 
of films of cellulose derivatives has been studied in detail. The solu- 
tions are peored into glass Petri dishes with and without a Cellophane 
lining. The film in contact with the glass provides information on the 

esive property of the composition and is also used to follow the 
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rate of drying. The film in contact with Cellophane is used to obtain 
an estimate of the amount of shrinking characteristic of the dope; in 
only a few instances has tenacious adherence of the film to Cellophane 
been noted. These tests are still in progress, but it may be noted here 
that they indicate that in order to an & maximum tautening effect 


it is necessary to formulate a dope so that a minimum of active solvent [ 


will be present during the final drying stage. . Some solvent action jg 
necessary in this stage to prevent precipitation (blushing) of the 


cellulose derivative. The selection of this solvent is an important 4 


factor in avoiding the formation of a brittle film, which Sheppard and 
Newsome * found to be associated with the presence of large crystallites, 
Thus, cellulose acetate dissolved in Methyl! Cellosolve gave the maxi. 


mum tautness for this cellulose derivative (table 5), but reference to | 
table 3 shows that the film is brittle and ringworms after a few days’ [ 
exposure on the roof. Likewise, cellulose acetobutyrate in formuk [7 
VII (panel 52), which contains 25 percent of Methyl Cellosolve, tautens | 
well but yields a brittle film. However, although our tests have not [ 
yet been completed, present indications are that it will be possible to [7 
formulate dopes with these cellulose derivatives which will tauten the | 


fabric satisfactorily without yielding brittle films. 


IV. EFFECT OF WEATHERING ON TAUTNESS OF DOPED | 


‘FABRICS 


In order to obtain comparative information on the effect of varia- 
tions in temperature and humidity on the tautness of fabrics doped 
with various cellulose derivatives, it would be desirable to have all the 
panels at about the same initial tautness. This was not achieved with 


the panel tests reported in table 3. However, the tautnesses obtained | 
with the cellulose compounds, using the solvent blend listed as formula | 
VII, were more nearly alike than with any other composition. The 7 


variation in the tautness values obtained for these panels under differ- 
ent weather conditions is shown in figure 3. 

The cellulose acetate panels were not markedly affected by the 
heat of the sun, but did slacken during periods of rain to a greater 
extent than any of the other cellulose derivatives. The cellulose 
acetopropionate panels also showed marked slackening in rainy 
weather after 3 months of exposure on the roof. The cellulos 
acetobutyrate and cellulose nitrate panels behaved about alike, 
becoming less taut when removed from the conditioning room into 
the sunlight and, in general, showing a slight additional decrease m 
tautness in rainy weather. The ethylcellulose and benzylcellulose 
panels were considerably poorer in initial tautness than those covered 
with the cellulose esters. Both of the cellulose ethers underwent 
decrease in tautness in sunlight. During periods of rain the ethyl- 
cellulose panels showed additional slackening, whereas the bengjl- 
cellulose panels had a tendency to tighten. 

It should be emphasized that the behavior of these cellulose films 
under various weather conditions with respect to both tautness 
moisture absorption is dependent to a large extent upon the type 
amount of plasticizer employed. For the above tests triphen 


*s. E. Snepeed and P. T. Newsome. Film formation with cellulose derivatives, J, 30c. Chem. Ind. 5, 
256T-261T (1937). 
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phosphate, present to the extent of 10 percent of the total nonvolatiles, 
was used throughout. In the course of the phase of our investigation 
dealing with the effect of plasticizers on tautness, both initially and 
upon exposure, compositions should become available which will 
have greater uniformity of tension under varying weather conditions. 


Kline 
Hie) 


CELLULOSE CELLULOSE CELLULOS 


NITRATE 


CELLULOSE 


ACETATE ETHYLCELLULOSE} GENZYLCELLULOSE 


WEATHER 
CLEAR 
RAIN 
CLEAR 
CLOUDY 
CLEAR 


CLOUDY 
RAIN 
CLEAR 


cLouox 
RAIN 


NUMBER OF DAYS EXPOSED 


CLOUDY 
CLEAR 


PANEL 13 PANEL 36 PANEL 52 PANEL 61 PANEL 66 


so 6100 so 86100 
OEFLECTION IN MILS 


Ficurn 3.—Variation in tautness of doped fabrics exposed on the roof. 


Formula VII was used for all of the dopes. ‘‘S’’ indicates slackness of the fabric. Weather conditions on 
the days when measurements were made were as follows: 





Tempera- 
ture 


Relative 
humidity 


Time 
on roof 


Tempera- 
ture 


Relative 
humidity 





°F. 


Percent 
74 35 
52 


NSRerees, 


























V. EFFECT OF WEATHERING ON FLEXIBILITY OF THE 
DOPED FILM 


The results of exposure tests on panels doped with the different 
cellulose derivatives, shown in table 3, must be considered as only 
tentative evidence of their relative stability. A formulation which 
= a flexible and durable film with one derivative may produce a 

rittle and unstable film with another cellulose compound because of 
differences in their solubility. However, certain trends can be 
observed in the data in table 3. The cellulose acetate films are, in 
general, very stable to sunlight. The only panels which gave evidence 
of embrittlement were three covered with clear dope made from low- 
viscosity cellulose acetate (panels 1, 5, and 25), one panel covered 
with clear cellulose triacetate dope (panel 23), and one panel covered 
with aluminum-pigmented dope Sets by dissolving cellulose acetate 


of medium viscosity in Methyl Cellosolve (panel 14), in which case 
the embrittlement is directly attributable to the solvent employed. 
The clear dope films made with cellulose acetopropionates of various 
viscosities failed on five (panels 36, 38, 42, 44, and 46) of the seven 





670 Journal of Research of the National Bureau of Standards [vu 


panels exposed. One panel covered with cellulose acetopropionate gf 
low viscosity and containing aluminum pigment (panel 43) failed jy 
17 days on the roof. However,'a sample of medium-viscosity celly. 


lose acetopropionate of more recent manufacture than the rest showe, [ 


good stability (panels 33, 34, and 35). The low-viscosity cellulos 


acetobutyrate dopes, both clear and pigmented (panels 48 and 49) | 
failed in 17 days on the roof, as did also a panel covered with medium. | 
viscosity cellulose acetobutyrate which was applied in a solvent mix. q 


ture containing 25 percent of Methy! Cellosolve (panel 52). Cracki 


of the films of clear cellulose nitroacetate (panel 55) and cellulog 
nitrate (panels 57, 59, and 62) was observed within 1 to 2 months of 


exposure. All 12 panels covered with clear ethylcellulose dope 


(panels 66, 68, 70, 72, 74, 76, 78, 80, 82, 84, 86, and 89) failed after re. | 
tively short exposure to the weather. Two panels covered with [ 
aluminum-pigmented dopes consisting of low-viscosity ethylcellulose; | 
dissolved in a mixture of acetone, ethyl acetate, toluene, and diace | 


tone alcohol (panels 73 and 75) failed in less than 2 months, wheres 
panels prepared with low-, medium-, and high-viscosity ethylcelluloses 
dissolved in this same formula (panels 67, 69, 71, 77, 79, 87, and 90) 
had not cracked after 3 months’ exposure. On the other hand, three 
panels of pigmented medium-viscosity ethylcellulose applied in a mix. 
ture of toluene, ethyl alcohol, and acetone (panels 81, 83, and 85) 
failed in 17 days, which indicated that this particular formula appar. 


ently tends to lay down a film which is inherently brittle. The clear § 


benzylcellulose film (panel 91) cracked and yellowed very soon after 
exposure. The pigmented benzylcellulose films (panels 92, 93, and 94) 
showed evidences of cracking after about 2 months of exposure, 
although the fabric was not tautened sufficiently for the typical 
“ringworm” type of failure to occur. The panels covered with methyl 
cellulose fondle 64 and 65) showed remarkable tautness during the 
first week of exposure, but the methylcellulose was completely removed 
by the first rain. This action was, of course, anticipated, but the 
rapidity with which the aluminum-pigmented methylcellulose film 
was dissolved by the rain was rather surprising in view of the compara- 
tively slow rate of solution of methylcellulose in the laboratory at 
ordinary temperature. The film of clear chlorinated rubber of low 
viscosity (panel 95) yellowed and disintegrated very rapidly upon 
exposure to sunlight. The pigmented film (panel 96) cracked ins 
manner similar to the benzylcellulose films after about 2 months’ 
exposure. The pigmented methyl methacrylate film (panel 97) did 
not show any evidence of cracking after 3 months on the roof. 


VI. SUMMARY AND CONCLUSIONS 


1. The tautnesses of airplane fabric doped with various plastics | 


dissolved in a variety of solvent mixtures were determined. It was 
observed that the most important single factor involved in the initial 
tautening property of a dope is the solvent composition. In order 
to obtain a maximum tautening effect, it is necessary to formulate 
a dope so that a minimum of active solvent will be present during 
the final drying stage. The selection of this solvent is also an m- 
poasene factor in avoiding the formation of a film which is initially 

rittle or which rapidly becomes brittle upon exposure out of door. 
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2. The highest initial-tautness values were obtained with cellulose 
triesters, such as cellulose triacefate and a practically completely 
acylated cellulose acetobutyrate. Varying the acyl or ethoxyl 
content of partially hydrolyzed cellulose derivatives did not have a 
ronounced effect on the ability of the compounds to tighten the 
abric. The initial tautening property is also apparently independent 
of the size of the cellulose molecule, as indicated by certain flow 
. characteristics of solutions of these materials. The tests for the 
majority of the cellulose esters were made with films containing 10 
percent of triphenyl phosphate. 

3. In exposure tests the cellulose acetobutyrate and_ cellulose 
nitrate panels behaved quite similarly, slackening somewhat when 
removed from the conditioning room into the sunlight, and, in general, 
showing a slight additional decrease in tautness in rainy weather. 
The cellulose acetate panels slackened during periods of rain to a 
greater extent than any of the other derivatives. The cellulose 
acetopropionate panels also showed marked slackening in rainy 
weather after 3 months of exposure on the roof. The ethylcellulose 
and benzylcellulose panels were considerably poorer in initial tautness 
and slackened still further upon exposure. However, the results of 
these exposure tests must be considered as only exploratory and 
preliminary to the testing of dopes formulated to develop optimum 
tautness, flexibility, and moisture resistance. Before dopes having 
these characteristics can be formulated, it will be necessary to obtain 
detailed information on the effect of various solvents, diluents and 
plasticizers on the properties of the film-forming plastics. 


This investigation was sponsored by the Bureau of Aeronautics, 
U.S. Navy Department, and the results are published by permission 
of the chief of that Bureau. The authors express their appreciation 
of the interest and suggestions during the course of this work of Lieut. 
Comdr. C. F. Cotton and J. E. Sullivan of the Bureau of Aeronautics, 
and of C. W. Woodrow of the Naval Aircraft Factory. 
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MINIMUM PERCEPTIBLE COLORIMETRIC PURITY AS A 
FUNCTION OF DOMINANT WAVE LENGTH 


By Irwin G. Priest and Ferdinand G, Brickwedde * 


[Eprror1aL N ote.—In February 1926, Priest and Brickwedde presented the 
results of their study of colorimetric purity before the Optical Society of America 
and an abstract, carrying the title given above, appeared in the J. Opt. Soc. 
Am. and Rev. Sci. Instr. 13, 306 (1926). The work had not been prepared for 
publication at the time of Mr. Priest’s death, and Dr. Brickwedde is now working 
in another field. Therefore, at the suggestion of the Director of the National 
Bureau of Standards, the writer has undertaken the preparation of this article.— 
Deane B. Judd.] 


ABSTRACT 


A 4° square, two-part photometric field, symmetrical about a vertical division 
and viewed through a pupil 3 mm in diameter, is illuminated in both parts by 
artificial sunlight at a constant brightness of about 3 or 4 mL (retinal illumina- 
tion, 70 to 90 photons) with a surrounding field of about 0.6 mL. Homogeneous 
light is added to one half, and a. simultaneously subtracted so that the 
field remains matched in brightness. oO adjustments of the mixture are made: 
(1) The least purity perceptible with certainty (pmax), and (2) the greatest im- 
perceptible purity (Pmin). The purity of these mixtures is then measured, in- 
creased accuracy being obtained by measuring a known large multiple of the 
homogeneous brightness. Values of Puss and Pmin have been obtained as a 
function of the wave length of the homogeneous component; these values are 
reported in detail, and some discussion of their interpretation is given. 


CONTENTS 
I. Introduction 
II. Procedure 
1. Adjustment of the mixture 
2. Photometric measurement of the components 
3. Precautions and auxiliary studies 
IIT, Results 


1. Comparison with other work 
2. Interpretation 


I. INTRODUCTION 


This work is undertaken as a part of a program of research on 
determination and correlation of the characteristics of vision out- 
lined by Priest 1 in 1921 when he listed as one of the functions to be 
determined, ‘Saturation sensibility—the least change in purity per- 
ceptible as a change in saturation.” This work deals with saturation 
sensibility near zero purity.? It is therefore of particular interest in 
the colorimetry of illuminants closely resembling sunlight in color 


*Research Associate at the National Bureau of Standards, representing the Munsell Research Laboratory, 
Jah 1925 to June 1926. 


G. Priest, The spectral distribution of energy required to evoke the gray sensat BS Sci. Pap. 17 

7H} (12) S417 (see Appendix, p. 260). 1 il ora emmection, p. 1%, 
‘For a definition of colorimetric ity, see Irwin G. Priest, Apparatus for the determination of colcr in 
terma of dominant wave length, purity, and brightness, J. Opt. Soc, Am, Rev. Sci. Instr. 8, 173 (1924); 
Irwin G. Priest, The computation of colorimetric purity, J. Opt. Soc. Am, and Rev. Sci. Instr. 9, 503 (1924). 
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and in the colorimetry of nearly achromatic surface colors views) 
under such illuminants; it bears importantly upon the question of 
tolerances in the specification of such colors and upon the determing. 
tion of nearest color temperatures * by means of which such illuminant 
are customarily specified. 


II. PROCEDURE 
1. ADJUSTMENT OF THE MIXTURE 
A 4° square, two-part photometric field, symmetrical about 4 


vertical division and viewed through a pupil 3 mm in diameter, was F 


illuminated in both parts by Abbot-Priest sunlight * at a constan 


brightness of about 3 or 4 mL (retinal illumination, 70 to 90 photons) [ 
with a surrounding field of about 0.5 mL.’ Homogeneous light wa f 
added to the right half of the field and sunlight simultaneously sub. 
tracted so that the right half of the field remained matched in bright. [ 


ness with the unchanged left half. Two adjustments of the mixtuy 
in the right half-field were made. First the proportion (p) of homo. 


geneous light in the mixture was increased until the observer per. | 
ceived a definite chromaticity difference between the sunlight in theleft | 
half and the mixture in the right; this is called a setting of the least | 


purity perceptible with certainty (Pmax). The observer then decreased 
the proportion of homogeneous light in the mixture until the chro. 
maticity difference just disappeared. This second adjustment is 


called the greatest imperceptible purity (Pum). The average be ff 


tween a setting of least perceptible difference and a setting of greatest 


imperceptible difference is thought to be a more reliable indication of 7 
sensibility than either one taken alone; it has sometimes been called 
the “doubtful difference.”* Accordingly, this average has also been | 


computed and is referred to as (Pmean)- 


2. PHOTOMETRIC MEASUREMENT OF THE COMPONENTS 


After each adjustment (either of least perceptible or greatest im- 


perceptible difference) the observer made a photometric measurement 
of the components of the mixture from which the purity (pmax OF Pun) 
may be computed. One of two photometric methods was used. 


In the equality-of-brightness method, the nonselective filter in the 


homogeneous beam was removed. Since this filter was usually either 
of 5- or 10-percent transmission, the homogeneous component by this 
— was increased either twentyfold or tenfold. Then the 
observer decreased the sunlight component of the right half of the 
field until the brightness was again judged equal to that of the left 
half which remained unaltered throughout. The reduction in su- 
light component was accomplished by increasing the distance of the 
lamp from a diffusing screen; it may therefore be expressed as a simp! 
function of this distance (essentially the inverse-square law). 
reduction is equal to a known multiple (usually about 19 or 9) of the 
homogeneous brightness, and thus permits the ratio of homogeneots 
3 R. Davis, A correlated color temperature for illuminants, BS J. Research 7, 659 (1931 Hoy: 


)R 
‘Irwin G. Priest, A precision method for producing artificial daylight, Phys. Rev. [2] 11, se 18); J, Opt. 


19 
Soc. Am. and Rev. Sci. Instr. 12, 479 (1926). D. B. Judd, Reduction of data on mirture of color stimu, F 


BS J. Research 4, 525 (1930) RP 163. on 

‘Irwin G. Priest, Apparatus for the determination of color in terms of dominant wavelength, purity, 
brightness, J. Opt. Soc. Am. and Rev. Sci. Instr. 8, 173 (1924). m1 
. enn). B. Judd, Sensibility to color-temperature change as a function of temperature, J. Opt. Soc. Am. 
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brightness component to brightness of the mixture to be computed. 
This ratio is the colorimetric purity of the mixture. 

The equality-of-brightness method is convenient but it is open to 
two criticisms: First, the field size (4° square) is larger than that to 
which the standard luminosity function recommended by the Inter- 
national Commission on IJumination applies,’ and second, the observ- 
er has to make the photometric match with a considerable chroma- 
ticity difference (sometimes amounting to 20 least perceptible differ- 
ences) between the two halves of the field. Accordingly, for about 
half of the work, the brightnesses of the components of the mixture 
were obtained by flicker photometry. 

In the flicker method a 2° diaphragm was substituted for the rec- 
tangular 2 by 4° opening constituting the mixture field. The non- 
selective filter in the homogeneous beam was removed, and the flicker 
disk interrupting the homogeneous beam with sunlight of brightness 
variable by an auxiliary photometer was set into rotation at a suitable 
speed. Adjustment of the auxiliary photometer for minimum flicker 
yields a measure of the brightness of this known multiple of the homo- 
geneous component. From this brightness and the known constant 
brightness of the photometric field the purity of the mixture (pmax or 
Puin) Was computed. 


3. PRECAUTIONS AND AUXILIARY STUDIES 


The wave-length scale was calibrated by reference to the lines of 
the helium emission spectrum. Stray light in the spectrum was 
measured and stray-light filters used in the extremes of the spectrum 
over the whole region within which stray light would introduce appre- 
ciable errors. This is particularly important for the short-wave 
region of the spectrum. 

A special check was made (observer IGP) of the slit-width errors. 
These errors can be shown by computation to be negligible in com- 
parison with the inherent uncertainty of liminal settings even for the 
rather large slit widths used (10 mu). The experimental check cor- 
roborated the computations. 

A special study was made (observer FGB) of the variation of p,,. 
and Pmin taken on different days over a period of 6 weeks. It was 
found that there was no regular variation from day to day, but the 
mean for a given day varied on the average from the mean on another 
day by more than would be expected from the theory of errors. This 
is probably to be ascribed to a change in criterion of what constitutes a 
“just perceptible” and a “just imperceptible” difference; such varia- 
tions are common, and indeed, characteristic of liminal settings. The 
conclusions from this study were used in a special determination of 
the wave length of the maximum value of Pmen (observer, FGB). 
This determination was carried out by making one setting only of 
Pusx 20d Pom at one wave length (565.0 mu) somewhat lower than 
that which preliminary determination indicated for the maximum, 
then a single determination approximately at the wave length of the 
maximum (567.5 my), and finally a single determination at a wave 
length (570.0 my) slightly above that for the maximum. After a rest 


'K. 8. Gibson and E. P. T. Tyndall, The visibility of radiant energy, BS Sci. Pap. 19, 131 (1923-24) 8475. 
marian Comm. on Ilium., 6th meeting, Geneva, 1924 (University Press, Cambridge, England) p. 67 and 


‘ * See, for example, D. B. Judd, Precision of color temperature measurements wnder various observing condi- 
ions; a new color comparator for incandescent lamps, BS J. Research 6, 1164 (1930) RP252. 
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pated the observer repeated the series starting with the higher way, 
ength. In this way, changes in criteria were avoided (long-time fly. 
tuations), and practice and fatigue effects (short-time fluctuations) 
reduced. From this work it was concluded that the wave length 
which mean is @ Maximum is for this observer (FGB, left eye) 568.04 
1.0 mu. 

III. RESULTS 


Figure 1 shows on a logarithmic scale as a function of wave le 
the result of every separate determination taken by IGP by the flicker 
method. Figure 2 similarly shows those taken by FGB. 

Table 1 gives the average values Of Pmean for each observer. These 
averages were obtained by plotting points representing the average 
of the groups of points shown on figures 1 and 2 on a logarithmic scaly 
together with similar points derived from data obtained with the [ 
equality-of-brightness method. Smooth curves were drawn among | 
these points, greater weight being given to the flicker determinations, 
The flicker method is not open to the two criticisms of the equality. 
of-brightness method discussed in the previous section, and it hes the F 
further advantage of greater precision in the measurement of the | 
smaller homogeneous Sebolimanaees. This greater precision is par- 
ticularly marked for the extremes of the spectrum (400 to 450 my, 
630 to 680 my); in these regions it sometimes happened that the 
equality-of-brightness method even gave negative values for pay in | 
spite of the nature of the experiment, which insures that some amount 
of homogeneous component greater than zero must be added. It wil 
be understood, however, that even when increased by ten- or twenty- 
fold these very small homogeneous brightnesses added to the sunlight 
component may fall below the differential limen for equality-o- 
brightness settings. It may therefore be concluded that the addition | 
to sunlight of small amounts of light from the extremes of the spectrum 
are detectable as chromaticity differences long before they become 
detectable as brightness differences. 

The only regular difference observed between the flicker and 
equality-of-brightness methods is that for the long-wave half of the 
spectrum the ficker method gave generally higher results. This is to 
be expected from the smaller field size, and since greater interest 
attaches to photometric results in accord with this smaller field size, 
the general course of the flicker determinations was followed closely 
in drawing the curve. The greater weight attached to the flicker 
determinations was therefore of two kinds: That attached to mor 
representative photometric observing conditions, and that to deter 
minations of higher precision. 
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TABLE 1.—Minimum perceptible colorimetric purity 





Pmean Pmean 
Average Wave Average 
Pmean length Pmean 
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The values of the first two columns in table 1 were computed from 
values Of Pmax 2Nd Pmin read from these adopted curves. . The values 
in the last column of table 1 are the arithmetical means of the values 
in the first two columns. 

Figure 3 shows ;(solid curve) a plot on a logarithmic scale of a 
smooth curve drawn through the values given in the last column of 

able 1. It also shows curves of pmsx (dotted curve, arithmetical 
mean Of Pmax for IGP and pmax for FGB) and pain (dashed curve, 
arithmetical mean of Pui for IGP and pain for FGB). 


IV. DISCUSSION 


1. COMPARISON WITH OTHER WORK 


When these data were obtained there was only one investigation 
with which to compare them; this was the pioneer work of Jones and 
Lowry’ on saturation scales extending from zero purity to the spectrum 
(unit purity) for eight different wave lengths. Although subsequent 
repetition of this work has in the main corroborated it, the variation 
with dominant wave length of the first saturation step from zero 
purity found by Jones and Lowry incidental to their main purpose is 
much less than that reported here. For example, Jones and Lowry 
give 0.035 as the size of the first purity step for 575 my and 0.025 for 


‘L. A. Jones and E. M. § 
Instr, 13, 25 (1920) Lowry, Retinal sensibility to saturation differences, J. Opt. Soc. Am, and Rev. Scly 
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440 my; this is a variation by a factor of about 1.4, while from table 
the corresponding variation is by a factor of about 18. 

Since 1926, however, Pmax has been redetermined by no less thay 
three independent investigations, those by Purdy.” by Martin 
Warburion, and Morgan," and by Wright and Pitt, Purdy inyegij. 
gated the change in Pass with field brightness, and therefore extended 
the scope of his wor 

particular; but he used stimuli of relatively few wave lengths (8) oye 

10 


LEAST PERCEPTIBLE 
COLORIMETRIC PURITY 


Average Results of 
Priest and Brickwedde 
Least purity perceptible 
with dertainty (.,.) 


——-— Greatest imperceptible 
purity (Pa) 


P; es Prax + P win 


450 5 600 700 
WAVE LENGTH mu 


Fieure 3.—Minimum perceptible colorimetric purity, average results of Priest and 
Brickwedde. 


the range 400 to 660 my.* Martin,“Warburton, and Morgan supply 
data for three observers, both for light adaptation and dark adapte- 
tion; their investigation refers to 21 wave lengths distributed over the 
range 440 to 680 my. Wright and Pitt obtained data at 19 wave 
lengths distributed uniformly over the range 470 to 650 my. In each 
of these investigations the results completely corroborate within 
experimental error those previously obtained by Priest and Brick 
wedde, with exceptions arising only from the fact that the Priest 
Brickwedde results cover a greater wave-length range more thoro 
than any other investigation. 

 D. MeL. Purdy, On the saturations and chromatic thresholds of the spectral colours, Brit. J. Psychdl 
(Gen. Sec.) 21 (pt. 3), 283 (1931). 

uL. C. M ’ L. Warburton, and W. J. Morgan, Determination of the Sensitiveness of the Eye! 
Differences in the Saturation of Colours, Medical Research Council, Reports of the Committee upon tbe 


Physio of Vision, XII1, Special Report Series, No. 188 (London, 1933). 
7 W. D. Wright and F.'H. G. Pitt, The saturation-discrimination of two trickromats, Proc. Phys. 80 


(London) 49, 329 (1937). 
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There can therefore be no doubt remaining as to the essential correct- 
ness of the data presented here completely for the first time.” 


2. INTERPRETATION 


Through the courtesy of Mr. Priest these data were, in advance of 
publication, placed at the disposal of Dr. Selig Hecht and the writer. 
By an extension of the Young three-component theory of vision 
Hecht has shown that these data are consistent with sensibility to 
wave-length change in the spectrum and roughly consistent with the 
mixture data of the OSA “excitation” curves, used for many years as 
a definition of normal vision.'* These data have formed an important 
basis for Hecht’s further theoretical papers, in which he has shown 
that they may be completely harmonized with the OSA curves.” 
Judd has shown by an empiric relation that these data are consistent 
with many other chromaticity-sensibility functions and are approxi- 
mately consistent both with the mixture data and coordinate system 
of the OSA “excitation” curves,'* and he has also shown by a special 
transformation of coordinates that these data are in precise agreement 
with the mixture data of the standard observer adopted in 1931 by 
the International Commission on Illumination.” 

From their appearance in summarized form in these publications, the 
Priest-Brickwedde data have continued to spread ® and in this way 
have already exerted a considerable influence on modern developments 
in visual theory. 

It is probably worth noting that Priest for some time held the view 
that Pmean Should duplicate very closely the luminosity function. This 
caused him to remark ™ “Least perceptible purity is largely dependent 
upon visibility; the relative amounts of energy required to introduce 
a hue difference show comparatively small variation with wave 
length,” and to view with suspicion the presence (see fig. 3) of the 
marked maximum quite definitely removed (570 mu) from the wave 
length of maximum spectral luminosity (555 my). He was inclined to 
the belief that if a truly achromatic stimulus were to be substituted 
for the rather greenish-yellow artificial sunlight then this pronounced 
maximum at 570 my would no longer be found. This belief led him to 
remark, ‘Further extensive investigation of least perceptible purity 
should be prefaced by a definite adoption of a neutral standard.” 

_ The present view of this maximum at 570 mz is that it indicates sun- 
light to be more confusible in chromaticity with 570 mz than with any 
. '’ Martin, Warburton, and Morgan (see footnote 11, p. gy | ae a discussion of possible sources of unre- 
— ad orn eyes data. See, also, D. B. Judd, Chromaticity sensibility to stimulus differences, 

u Selig Hecht, The development of Thomas Young’s theory of color vision, J. Opt. Soc. Am. 20, 231 (1930). 
The tentative theoretical curves in this paper do not completely harmonize mixture data for the normal eye 
with the Priest-Brickwedde data; consult D. B. Judd, The mixture data embodied in the tentative curves of 
Hecht’s theory of vision, J. Opt. Soc. Am. 20, 647 (1930); this deviation has been taken care of in Hecht’s subse- 
quent formulation of Young’s idea. 

Selig Hecht, The interrelation of various aspects of color vision, J. Opt. Soc. Am. 21, 615 (1931); A Quanti- 
tative Formulation of Colour-Vision, Report of a Joint Discussion of Vision, The Physical and Optical Socie- 


ties, June 1932; The Retinal Processes Concerned with Visual Acuity and Color Vision, Bulletin No. 4, 
Howe Laboratory of Ophthalmology, Harvard Medical School (Harvard University Press, Cambridge, 


ass., A it 1931. 

- D.B. Juda, Chromaticity sensibility to stimulus differences, J. Opt. Soc. Am., 22, 72 (1932). 

D. B. Judd, A Marwelil triangle yielding uniform chromaticity scales, J. Research NBS 14, 41 (1935); 
RP766; J. Opt. Soc. Am. 25, 24 (1985). 

‘See, for example, J. F. Schouten, wee einer ntitativen Vierfarbentheorie. I, Proc. Konink 
Akad. Wetenschappen Amsterdam 38, No. 6 (1935). abeuten, like Hecht, has shown an approximate 
ato ae mixture data for the normal eye and the Priest-Brickwedde data by resort to a 

isual theory. 

“TI. G. Priest and F. G. Brickwedde, The minimum perceptible colorimetric purity as a function of dominant 
wave length with sunlight as neutral standard, J. Opt. Soc. Am. and Rev. Sci. Instr. 18, 306 (1926). 
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other part of the spectrum. On the basis of this conclusion, a ney 
coordinate system was derived,” which has not only been successful in 
coordinating chromaticity sensibility data of various sorts but has 
also proved to be useful in practical specification of color tolerance 
and in the estimation for practical purposes of nearest chromaticity 
from a chromaticity series." The success of these colorimetric pro- 
cedures can largely be traced to the Priest-Brickwedde data. 


WASHINGTON, January 2, 1938. 


2» D. B. Judd, A Marwell triangle yielding uniform chromaticity scales, J. Research NBS 14, 41 (1935) RPM; 
J. Opt. Soe. Am. 25, 24 (1935). 

2D. B. Judd, A method for determining whiteness of paper, Paper ane J., Tech. Sec. 100, 266 (1935); 1m, 
154 4 (1990); Tech. Assn. Papers, series 18, 392 (1935); series 19, 359 (1936 

J. MeNicholas, Color and spectral transmittance of vegetable oils, J. Research NBS 15, 99 (1936) RP3u; 

oil AS Soap 12, 167 (1935). 

D. B. Judd, Estimation of geen differences and nearest color ae on the standard 1981 1¢] 
colorimetric coordinate system, J. Research NBS 17, 771 (1936) RP944; J. Opt. Soc. Am. 26, 421 (1984). 

H. J. MeNicholas, Selection o; [fom for signal lights, J. Research NBS i, 955 (1936) RP956. 

= B. Judd, Surface color, J ge Soc. Am. 25, 44 (1935). 
i ~ zt and W. R. Schaub, Rectangular uniform-chromaticity-scale coordinates, J. Opt. Soe, 

m. 27, 226 (1 

D. L. MacAdam, Projective transformations of ICI color specifications, J. Opt. Soc. Am. 27, 294 (1937), 





U. §. DEPARTMENT OF COMMERCE NATIONAL BuREAU OF STANDARDS 


RESEARCH PAPER RP1100 


Part of Journal of Research of the National Bureau of Standards, Volume 20, 
May 1938 


APPLICATION OF VERTICAL-INCIDENCE IONOSPHERE 
MEASUREMENTS TO OBLIQUE-INCIDENCE RADIO 
TRANSMISSION 





By Newbern Smith 


ABSTRACT 


An extension of the transmission curves described in a previous paper is given. 
The effects of the earth’s curvature and magnetic field are considered, as well as 
absorption or reflection from lower ionosphere layers. A chart is included for 
rapid calculation of transmission curves. 
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I. INTRODUCTION 


In a recent number of this Journal! the author outlined a graphic 
method of obtaining, from vertical-incidence ionosphere measure- 
ments, the limiting frequencies and virtual heights of reflection of 
radio waves incident obliquely upon the ionosphere. In that paper 
was developed a type of “transmission curve’ which, when super- 
imposed on a curve of frequency against virtual height, measured at 
vertical incidence, gave directly the virtual height of reflection for a 
given frequency at a given distance. A logarithmic sec ¢) trans- 
| mission curve was also described by which, for a given distance, one 
curve could be used to determine the heights, to a fair approxima- 
tion, at any frequency. 

It is the purpose of the present paper to evaluate in more detail 
the effect of the earth’s curvature upon the ray theory of ionosphere 
transmission, to indicate the modifications consequently necessitated 
in the application of the transmission curves, and to consider briefly 
— earth’s field may affect limiting frequencies and virtual 

eights. 


II. EFFECT OF THE EARTH’S CURVATURE 
The influence of the earth’s magnetic field will be neglected in this 
section. In the next section the effect of the magnetic field on radio 


'Newbern Smith. Extension of normal-incidence ionosphere measurements to oblique-incidence radio 
tranemission. J. Research NBS 19, 89 (1937) RP1013. 
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transmission will be considered for the case of the plane earth, and jt 
will be assumed that the results hold approximately for the curve 
earth, also. 

The elementary theory of propagation of electromagnetic waves jp 
an ionized medium whose surfaces of equal ionic density are planes 


‘ 
, N 


\ 
\IONOSPHERE 





Figure 1.—Variation, with height, of angle a straight line makes with the normal tp 
the earth’s surface. 


shows that, by Snell’s law, the waves will penetrate the medium until 
the refractive index yu’ is reduced to the value sin ¢;, where 4; is the 
angle of incidence of the waves upon the ionosphere, i. e., the angle 
the wave normal makes with the vertical on entering ,the ionosphere. 
To obtain the corresponding relation for a curved ionosphere it is 
necessary to consider the variation in the angle ¢ which a straight line 
makes with the vertical, or normal to the earth’s surface, at various 
altitudes (z). 
The geometry of figure 1 leads to the relation 


sin ¢’=sin ¢ (1+ 95575): () 


Assuming Snell’s law to be valid for a ray traversing an infinitely 
thin layer of the ionosphere, of thickness dz, we can use this relation 
to obtain the differential equation: 


d(u’ sing) _—_ dz 
yw sing  R+h+2 


Integrating this from the lower boundary of the ionosphere, where 
zp’ =1, ¢=¢, and z=0 (see fig. 2) up to the altitude z, we get 





sin 1 
l+Rar 
as the form of Snell’s law appropriate to the curved earth. 


z’ sin = 





smith} Ionosphere Measurements 685 


This means that the wave will penetrate the curved ionosphere 
until the refractive index »’ is reduced to the value po’, given by 


sg! = Si 
2 (3) 
1+RTy, 


where ¢; is the angle of incidence of the wave on the lower boundary 
of the ionosphere and 2 is the maximum height of penetration above 
this lower boundary. 

From the geometry of figure 1, we have also 


sin ¢,;=sin (1 +R5)" (1a) 


If we put ¢’ = and z=z2,, the values which these quantities have at 
the top of the equivalent triangular path in figure 2, this becomes 


sin ¢;=sin o(Itp ty ; (4) 


ts . 
¢ 
s 


R 
a ee ae | 
— eae : 
FicurE 2.—Transmission through curved ionosphere. 


2=height of equivalent triangular path; 2o=true height of reflection, at level where y’=sin ¢; ( 1 75) ’ 
R=radius of earth; and D=distance of transmission. 


and so the refractive index at the level of reflection must be 


(5) 


instead of uo’ =sin ¢o, as is the case for the plane earth. 
If, as is usually the case, the penetration ~4<<R+h, we may write 


this condition 
yo’ =sin (1+ 554) (5a) 
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The effect of the earth’s curvature is thus to cause the radio way, 
to be reflected from a region of smaller ionization density than would 
be necessary were the earth flat. 

In considering the problem of radio transmission over a distange 
it will often be convenient to speak of an “equivalent vertical-ing. 
dence frequency’’ for the given wave frequency and transmission path, 
This is the frequency which is reflected at vertical incidence at the 
same height and undergoes approximately the same absorption as does 
the actual wave at Solas incidence. It will be denoted by f, and 
the actual wave frequency will be denoted by f’. 

Since, for a frequency f’, 


is Ne* _.|/,—Jo 
V — may — 
where fo= xe 23 it follows that reflection of this frequency will take 


Ca 
po =a). 


where f,’ is the value of f, at the level of reflection, corresponding to 
an ionic density of Np. 





place at a level where 


Noe? 
mm 


f= 





’ 


and f is the frequency for which ».=0 at this same level of reflection, 
i. e., the equivalent vertical-incidence frequency. For a wave fre- 
quency f’ the equivalent vertical-incidence frequency f, in the absence 
of a magnetic field, is therefore 


f=f'V1—m”, (6) 


is the refractive index for the frequency f’ at the level of 


, 


where jo 
reflection. 
In the case of the plane earth, we had 


tgs uo’ =sin ¢=sin ¢o, 
which gives 





ie 
bye ty 


the well-known secant law. In the case of the curved earth, however, 
we must use the values of yo’ given by eq 3 or eq 5, and 








72 in? l+prh 
Vi—petme [3 em ;= | 1—sin? ¢&| ———] (7a) 

1+ Zo en. 

R+h R+h 


Under practical conditions 2 is almost always less than 400 km, 
and values of z, greater than 600 km contribute nothing to trans 


*? W. H. Eccles, Proc. Roy. Soc. Cantey) [A] 87, 79 (1912). 
4 J. Larmor, Jahrb. drab Telegraphie 


141 (1925). 
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ave — mission over at distances, so that we can always consider 
ud | zand2<&+h. Using this approximation, eq 7a becomes 





i | -Y1—%0?=c0s dry) 14 tan® ¢,=cos boy) 1-2 tan*¢,(7b) 
th. 
the Figure 3 gives values of earn 


bo” 


: y 


plotted against tan ¢, for various 
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1 
Fiaure 3.— Variation of Vi-n,? with z and tan 9». 
o 


ef, values of 2, assuming h=100 km (at the bottom of the EL layer). 
These curves may be used to determine the relation of f to f’ for 
various values of 2 and ¢;. By substituting ¢ for ¢, and a—z, 
for 2, similar curves may also be plotted to determine the relation 
of f to f’ for various values of 2, 2), and ¢po 

a) If, further, ~< (R+h) cot? ¢, or ae 2) <(R+h) cot?d) we 
may write eq th approximately, 


¥1—'2=cos | +p tan? 61 |=cos és] 1- Patan do | (7c) 





& form which it is convenient to use in some discussions. This ap- 
proximation leads to results good to 1 percent or better for E-layer 
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transmission, where 2 is less than 50 km, and for single-reflection 
F-layer transmission over distances less than 1,500 km, or multire. 
flection F-layer transmission where each reflection covers less 
1,500 km. For transmission over distances greater than 1,500 km 
for each reflection it is necessary to use the more exact expression, 

The correction term, or frequency Af which it is necessary to gyb. 
tract from f’/sec ¢) in order to obtain f, is the amount by which the 
equivalent vertical-incidence frequency for the curved ionospher 
differs from that for the flat ionosphere. 


This is 
af=f| cos go— V1 — po’? |? 
which is approximately 
Af=f’ cos of 1-12 FE tan? 60 | 


Ze— << (R+h) cot? go, 


. po Ze— & \ Sin? do 
Af=f R+hA / cos dp 


This correction term is zero at vertical incidence and for ~=2z,', i.e, 
reflection from a sharp boundary at the height 2. For 22,’ it in 
creases rapidly with distance. 

The vertical incidence frequency f derived above may be used to 
plot transmission curves of z,’ against f/, exactly as in the paper 
referred to in footnote 1. If the assumption be made that z, is the 
same as the z,’ measured at the frequency f, then these transmission 
curves could be superimposed on the curve of frequency against virtual 
height, measured at vertical incidence (hereafter to be called the 
(z,’, f) curve), to give directly the z, for the oblique-incidence case, 
The actual relation of z, to z,’ is discussed later. 

Actually, it is more convenient to use the sec ¢ transmission 
curves, wherein z, is plotted against f’/sec @o, and to refer the correc- 
tion to the (z,’, f) curve. On the assumption that z,=—z,’, A f can be 
calculated, for each point on the (z,’, f) curve, for each sec ¢p trans- 
mission curve. This transmission curve can then be superimposed 
on the (z,’, f) curve displaced toward the higher frequencies by the 
amount A f, and the oblique-incidence values of z, read off. 

The calculation of A f for each distance from z,’ and 2, without any 
assumption as to 2,, and for a ¢ calculated for an equivalent triangular 
path of height z,’, is discussed below. This value of A f can be used 
in the same manner as was described in the preceding paragraph, 
the sec ¢, curves applied, ¢) being calculated for an equivalent triangu- 
lar path of height z,’ instead of 2,. ‘ 

If the (z,’ f) curve and the sec ¢ curves are plotted logarithmically, 
soe the paper referred to,‘ we may calculate a factor 1+Af/f, a 
ollows: 





and, if 





Af cos do 
1l+— =———— 
“2 t a 1— po” 
ionosphere measurements to oblique-incidence radio transmission. J. Re 


4 Extension of normal-incidence 
search NBS 19, 89 (1937) RP1013. 
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The frequencies on the (z,’, f) curve can then be easily multiplied by 
this factor by adding the factor ay ager mete to the curve, and the 
log sec $ curves can be applied to the resulting corrected (2z,’, f) curve. 
It should be noted that in this expression the ¢, refers only to an 
equivalent triangular path of height z,’, and not z, The z, enters 
only into the expression for yo’, : 

A typical (2,’, f) curve corrected in this manner for @ given distance 
(2,000 km) is shown in figure 4. The virtual heights are lower on 
the corrected curve and the curve extends out to frequencies higher 
than the critical frequency for the ordinary ray. The correction to 


' the curve is different for different distances, increasing with the 


distance. This means that the virtual heights of reflection are lower, 
for greater distances, than if the curvature of the earth were not 
considered. The limiting frequency of transmission, or maximum 
usable frequency, over considerable distances is correspondingly in- 
creased, the difference becoming as much as 20 percent at the greater 
distances, departs of 

course on the form of the 
(z,’, f) curve. 

When the (2,’, f) curve 
has been corrected in the 
above manner it may be 
plotted logarithmically, as 
described in the previous pa- 

er mentioned above, and the o 
— = sec bo le ond 
be applied to give directly the ‘ ; 
maximum usable frequencies Figure 4.—(z,’, f) curve corrected for a given 
and virtual heights over trans- cor ga km a ee 
mission pathsof givenlengths. 

Since the first draft of this aper was written there has come to 
the author’s attention an excellent unpublished paper ‘“‘Skip-Distance 
Analysis,” by T. L. Eckersley and G. Millington, in the form of a 
contribution to the November 1937 London meeting of the Special 
Radio Wave Propagation Committee held in preparation for the Cairo 
Radio Conference. In this they undertake an analysis of radio trans- 
mission over a curved earth and obtain curves for determining the 
retardation of sky wave over ground wave at a distance, in terms of 
vertical-incidence measurements. They limit the analysis, however, 
to the case where 2 is very small compared with h. The analysis can 
thus apply only to E-layer transmission, since for F-layer transmis- 
sion, % must be measured from the lower boundary of the ionosphere 
in order to include the effect of retardation in the E layer. 

Their work, however, suggests a method of approximate analysis 
for larger values of 2 and a means for easily estimating the error 
involved in the assumption made above that the z, for a given fre- 
quency f’ was equal to the vertical-incidence 2,’ for the equivalent 
vertical-incidence frequency f. 

The development given by Eckersley and Millington begins with 
the following relation, for the ray path shown in figure 2: 


sa [Pap (Bassing 
, =f, a= f R+h+2’ 
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where ds is an element of the ray path making an angle with th 
vertical. Since, from eq 2 


this becomes 








They also used the relation 








L3-f[ dz - # 
An J, »w cosd J, y" sin? ¢; 
thal 2 2 
(1+p53) 


For the case where there is no magnetic field yp’ = i, where 


bs: 


f'dz 
[f*-10-9" 





sin’ ¢, 

2 2 
(1+953) 
Putting into this expression the wave frequency f’ in terms of the 

equivalent vertical-incidence frequency given above in eq 6. 

[°%-z43f'+3-- fdz re 

AK V¥1—p"*Jo J a yee sin’, 
ip" (i+ 2 y —fe (ii) 
R+h 

At this point Eckersley and Millington assumed that they were dealing 
with a thin layer, and that in consequence (1) 6’ was only slightly less 


(10) F 





sin ¢; (3ds ie: : 
than RthJa x” from eq 8, (2) aw wes only slightly greater than 


J **_fdz_ from eq 11. This simplification is justified for 


Vi—m'J> We 


E-layer transmission, but is unfortunately not justified for F,-layer 
transmission, since the lower boundary of the ionosphere must 

taken at the beginning of the E-layer in each case. We shall not, f 
therefore, confine our discussion to this limited case. It must also 
be noted that the simplification introduced by assuming << 
(R+h) cot? ¢; also is not valid except at comparatively short distances, 
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where ¢:=45° or less. For practical ea wr it is, however, justi- 
fable to assume << R+hA. Using this relation then, we obtain 


,__sin , B ds 22 

= Reha wR) 
| Eckersley and Millington combined the approximate forms of eq 8 
and 11 and introduced the value of /1—,’=cos ¢,, to obtain 
ya of"  fdz 
"Eth o WPI 


obtain the more exact relation 


We shal], however, combine eq 11 and 12 to 





sin ¢; 


sis, » (RN ee 
: aoa, WPa—A)—feN’ Rik)? (2) 








sin? ¢; ‘ ‘ 

2 Zz 
(+ ga) _ theta 4 
1— py”? 


1— 





1—A= 





2 
1 +p tan” dg; 


Gey 


If we let 
1+-3B a 





sin 1 





= R+h)y1— ma? fs J ; as (1p) (13) 


For values of z considerably less than (R-+-h) cot? ¢, B is usually a 
small number, and for z—>z (near the level of reflection) it is also 
usually small. It may attain, however, relatively large values. 

In general, the evaluation of B involves a knowledge of the dis- 
tribution of ionic density, and a precise evaluation of the quantit 
6’ involves graphic integration for each case considered. We shall 
for the present consider two principal cases, —0 and 2=z,’, the 
virtual height at vertical incidence for the frequency f=f’ J1—w”, 
ys We some cases of linear distributions of ionic density, where 
2/4=Jo . 

(a) a<(R+h) cot? ¢,. This case applies to E-layer transmission 
small distance or to F-layer transmission at short distances (¢, 
8 ‘ 


Combining eq 6 and 7 we get 


faf' cos o| 1+ py tan’ 6: | 
and from eq 12b 


1—A=142( fe) tan? 6. 
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From eq 12c 
2 


1 din aon 
ome We 20 tan? 
14+ B=1+( p5,)—73 tan’ 6, 

if 
if we assume that (1—2/2,)/(1—f,?/f”) is nowhere so large as to make R 
comparable with unity. This is equivalent to assuming that f/f 
does not approach 1 much more rapidly than does 2/2, an assumption [ 
that is valid except quite close to a critical frequency. Transmissions f 
involving equivalent vertical-incidence frequencies close to the critica] 
frequency are of interest, however, only when the distance of trans. f 
mission is short, and in this case ¢; is small, so that B is a small number, f 


anyhow. 
Thus 


tan¢df, —_% | a dz 
R+hL!_R+h Oo % J; ’ te 
Me 


6’ = 





Zz 


1+p77, tan’ ¢1 F\ R+ht 
(1-f) 





_=2,'» the virtual height measured at vertical inci. f 


dence for the equivalent vertical-incidence frequency f. Since the 
terms involving z/R-+-A are small, we may write 


z,’ tan ¢, 


=| et tan? 6,(1—0)—" |, (14e) F 


where 


Pe 
ee i nie ts , and 
Zo J0 Jo’ 1— 2 

fp 


fi 
O’= 2 f 2% 2dz 
Zy (R-+h) 0 te 
we 
For a linear distribution of ionic density 


z Ie’ , 
1 ma 2,,=2z%, and C=1, and 


of dz fe 22 "| 
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For % vanishingly small C’—>0 and eq 14 reduces simply to 


= tan 
~ R+h 


(b) 2-2’. This case applies to reflection from a fairly sharp 
boundary at the level 2, the refractive index being nearly unity up to 
nearly this level. An example of this would be reflection from the 
sporadic E region. The angular distance of the part of the ray path 
where » departs appreciably from unity is small and so the ionosphere 
can be considered as essentially flat. This is the condition assumed 
by the author in the paper already referred to (RP1013). Here 
| 2,=2, and 


(15a) 


tan ¢o=— in €. (16) 
Ranti—s 0 


+ 1—uo’?=cos o( 1-25 tan? #0) (16a) 


| For z,/—2 vanishingly small compared with (R+h) cot? ¢o, this 
reduces to the case for the plane ionosphere, where 





In this case 


/1—u9"=Ccos gy. 
(c) a= i. e., a linear distribution of ionic density. Here 
1—f,?/f?=1—2/z and eq 12a becomes 


= 22 
Q’ sin ¢; f de 1— pry 
0 


(R+-h)-¥1— po’ y1-4-2 
0 


Putting in the value of 1—A and integrating 


sin b4/ lps tan’ di, 9 2) 
(R+h) Ji—y” (1-3 ete 


Now since we may consider 2.<R-+A we can write, as in eq 7b 











, 
==2,’ 








V1—1"=cos br +7 tan? ¢, 
and thus eq 17 becomes 


oe tan ¢, _2 Be 
hee 22 [1 3R+h |’ (15) 


just as in the case where ~<(R-+-A) cot? ¢,. We must, however, 
write for the equivalent vertical-incidence frequency in this case 





t=f' cos or] 1 +R tan’ ¢, (17a) 
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instead of the value 


f=f' cos 6] 1 +p tan® | (14h) 


which we could use when ~<(R+h) cot? 9). 

The three examples just treated give an idea of what happens m 
transmission through the ionosphere, and of the ony distance 
6’ the wave travels in the ionosphere as a function of e of incidence 
¢, (or vertex angle of equivalent triangular path ¢), of true height of 
reflection z and of virtual height measured at vertical incidence z/ 
We must now consider the part of the path from the earth to the 
lower boundary of the ionosphere. If we consider the incident ray 
to traverse an angular distance (6—6’) in going from the earth to the 
— boundary of the ionosphere, the geometry of figure 2 tells ys 
that 

san doe (18) 
Rti—cos (@—0’) 





For a given ¢;, then, we may solve this equation for (@— 6’), and, 
by adding this angle to the angle 6’ already computed, we obtain the 
entire angular distance @ traversed by the wave from the ground to 
the point of reflection in the ionosphere. For h=100 km, as we ar 
assuming, it is sufficiently accurate to replace sin (@— 6’) by (6—@') EF 


and 1—cos (@— @’) by 5 (0— 6’)? in this equation. We can thus 
solve for (é— 6’), obtaining 


6— 0’ =cot 61-4) cot? — (18a) 


We can now write the expression for the total distance of trans- 
mission D in terms of ¢; (or ¢o) 2) and 2,’: 


D=2R(cot b—4 cot? bp) te (19) 


where 6’ is computed as above, from eq 12a. 
The time T required for the sky wave to travel from the transmitter 
to the receiver is, if c=velocity of the wave in vacuum, 


r=2(R sin ((—0’) #) 
c sin ¢; "ae 





If we put in the value of [aie obtained on the basis of the above 


analysis, and replace sin( @— 6’) by its value in terms of ¢,, we may 
express 7’ in terms of the quantities ¢, (or ¢o), 2 and z,’, as was done 
with D. 

The relation between the height z, of the equivalent triangular path 
and the virtual height z,’ measured at vertical incidence for the 
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equivalent vertical-incidence frequency f=f’./i—” can now be 
written. From the geometry of figure 2, 
: / 
tan ¢do=tan (¢,—6’)= ane , 


Futin cos 6’ 


»=(R+h)[sin 0 cot (¢,—6@’)—1+cos 6’] (20) 


The relation between z,, 2,’, 2, ¢1, do, 0’, D, f’, and f may be sum- 
marized for the cases discussed here. 
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me 
Bi 6’ cot (¢,—4’)—1+ cos 0’, 


ire 


f=f' cos di, 
g’ = 2 tan ¢; 
R+h ’ 
2, R(@—0’) 
cL sin ¢; 





+ 2 


a cos o;.} 


6—@’=cot 6:—a/cot? 6-3 


D=2Ré, and 
2,=(R+A)[sin 6’ cot (¢,—6’) —1+ cos @’], 








(b) linear gradient of ionic density («=32"’) 





f=f' cos oy lhe tan’? $1, 
= z,’ tan a ee Ze 
" Jee 3R+ hy 
R(e—6’) 
T== sin ¢; +r | 


6—@’=cot $:—/cot? od1— id 


D=2Ré, 
2,=(R+A)[sin 6’ cot (¢,—6’)—1+-cos 6’], 


f=f' cos do; 
tan 4 = sin 0 =; cure 0 
Rtgt!— cos 6’ —*__—+-1—cos 0 


ee (R+h) sin 6 


~ ¢€ gin do 
D=2Ré, and 


2,=2, 





(c) m=2,’, 








Case (a’) corresponds to the special case treated by Eckersley and 
Millington, except that, for simplicity of Porn Gen computation, they 
took R+h=R in the expression for 6, whic introduces an error of 
about 0.016h percent (A expressed in ’km)—not a serious error for 
h=100 km, as we have assumed. This error becomes appreciable, 
however, if it is attempted to extend this treatment, as they have 
done, to a height of 400 km or so. 

Martyn’s equivalence theorem,’ developed for a plane ionosphere, 
tells us that 

S.=8, , 
ot ee 
~ €COSd) CCOS 
D’=2z,’ tan ¢=2z,’ tan ¢,, and 
t=f' cos d=f' cos 41, 


where 


jens of incidence of waves upon the ionosphere, 
¢)=half the vertex angle of the equivalent triangular path, 
z,=height of equivalent triangular path, 
T’=time the wave spends in the ionos here, and 
D’=horizontal distance the wave travels in the ionosphere. 
We see that the relation ante, te valid in case (c) but not in any 


other case. The relation T=2 - ae is valid in cases (a’) and (6), but 


is not valid in general. The cahetion D’=2z,' tan ¢1 holds approxi- 
mately only in case (a’). And finally, f=’ cos ¢» in case (¢), f=f 


5D. F, Martyn, Proc. Phys. Soc. (London) 47, 332 (1935). 
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os ¢; in case (a’), and f is a complicated function in every other case. 
it is therefore concluded that the equivalence theorem, in the form 
given, cannot be applied to the curved-earth problem. 

Referring again to the summaries of cases (a’), (b), and (c), we may 
for a given value of 2,’ take a set of values of 2, for each of which can 
be calculated the variation of D with ¢,, in each of the three cases. 
We may also, for these values of 2,’ and z, calculate the variation of 
"f with ¢1. By eliminating ¢, graphically, we can determine the 
variation of f’/f with D for a given z,’ and 2, or, what is the same thing, 
the variation of f’/f with 2,’ for a given D and %. We can plot a 
family of transmission curves, with f as abscissas and z,’ as ordinates, 
for each value of D and several values of 2, corresponding to each of 

ese Cases. 

Oench curves may be superimposed on the (z,’, f) curve obtained at 
vertical incidence, and the (2, f) curve deduced therefrom, to give 
directly the wave frequency f’ corresponding to reflection at a given 
level 2, Which is characterized by a given z,’. The value of 2 appro- 
priate to reflection at the given height determines which curve of 
the family for a given D is to be used, and the point of reflection is 
determined as the intersection of this transmission curve with the 
2, J) curve. 

) ke value of z, depends only on D and 4,, and so lines of equal z, 
may be plotted on the curve sheet for each D, so that z,, as well as 
j', may be read off directly. 

Because of the approximations made in case (b) and the assumption 
of a linear variation of ionic density with height, this case is of only 
special significance. It will be assumed, until further investigation 
determines more precisely the variation of conditions with 2, that the 
curves vary smoothly between those calculated for 20 and those 
talculated for 2=2,’. 

A family of curves for each distance is rather cumbersome for 
rapid use. It is, as was said above, more convenient to use the log 
¢ ¢ transmission curves, and apply a correction to the (z,’, f) curve 
by multiplying each vertical-incidence frequency by the factor 
1+Af/f, where 


1 4 4f on do . 
tim 
This factor is obtained, for a given D, z,’ and 2, by determining cor- 
responding values of D and ~1—yo? for arbitrary values of ¢ or ¢. It 
8 unity for z,/=2) and is quite easily obtained for 2=0. For inter- 
mediate values of 2 it will be assumed that the factor 1+Af/f varies 


na manner similar to that determined from the relation eq 7b with 
4=2,,1. 6., 





008 $y __ _2(2s’—2) . 
Trt VS tat 
but drawn through the values for 2=0 and 2=z,’ determined in 


the more precise analysis, rather than those indicated on the assump- 


tion that z,=z,’. Cos ¢ is here calculated for an equivalent triangu- 
lar;path?of height z,’. ‘ie i 
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Figure 5 gives the approximate factors by which f must be mul. 
plied to give f’/sec $ for values of z,’ from 200 to 500 km and {y, 
distances up to 4,000 km. 


III. EFFECT OF THE EARTH’S MAGNETIC FIELD 


The presence of the earth’s magnetic field introduces some compli. 
cations in the use of these transmission curves. These complication; 
are often of minor importance compared with some of the unknow 
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FicurRE 5.—Approximate factors by which the equivalent vertical-incidence fr- 
quency f must be multiplied to give the frequency f'/sec 9 which is to be wsed in 
conjunction with the log sec do curves. 


factors (e. g., the geographic uniformity of the ionosphere over the 
transmission path), especially over long distances, but the effect of 
the earth’s field must at times be taken into account. The anisot- 
ropy of the ionosphere due to this field causes the effect of the field 
on radio transmission to vary with the length, direction, and ge0- 
graphic location of the transmission path. te 

One effect of the field is to cause the received signal to be split 
general into two main components, the one with the lower maximum 
usable frequency known as the ‘“o” component and the other s 
the ‘x’? component. The refractive index for a frequency /’, in the 
presence of a magnetic field H, whose components along and: trails 
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yerse to the direction of phase propagation are, respectively, H, 
and H, is given ® by 


sap alt i i alec 
Iarbas ty|| oes | +h fi? 
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> ce coat 
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e 
Su=xmc’ and 


N=ionization density. 
Upper sign refers to o-component; lower 
to z-component. 


The frequency of the wave whose z-component is returned from a 
given ionization density, at a given height, is different from the fre- 
quency of the wave whose o-component is returned from the same 
level. This frequency separation is in general a function both of fre- 
quency and distance and may often become negligible nee distances. 

For practical calculation, it may be assumed that only the field and 
direction of wave propagation in the region near the level of reflec- 
tion will appreciably affect the propagation of the wave. This assump- 
tion is probably better for the o- than for the z-component, and is, it 
must be emphasized, only a good approximation. 

With this limitation, therefore, the value of yu’? may be calculated 
for a given transmission frequency and transmission path. If we put 
u’=sin ¢) and deduce 2, from sin ¢o as was done in the previous paper, 
we may plot transmission curves, of virtual height against equivalent 
vertical-incidence frequency, for the o- and z-components. When 
these curves are applied to the corrected (z,’, f) curves, they may be 
expected to give reasonably good results. An example of this type 
of transmission curve is shown in figure 6. The one curve gives trans- 
mission conditions for the z-component, and the other for the o-com- 
ponent. The frequency used is well over the gyrofrequency f,= 
eH/2me so that the z-component is returned from a lower level than 
is the o-component and has a higher limiting frequency. 

It is not now justifiable to plot the transmission curves logarithmi- 
cally, since the form of the curves will vary with the transmission 
frequency f’. For practical purposes, however, a logarithmic curve 
may be used within a limited range of wave frequencies about the 
frequency f’ for which the curve is plotted; a practical limit might be, 
say, within +15 percent of this frequency. 

_ The logarithmic sec ¢» transmission curves may be used in estimat- 
ing the maximum usable frequency for each component over a given 
by adding or subtracting the separation between the limiting 
uencies for the two components, evaluated at that pogneney and 

i the trans- 


tance. This separation is in general a function only o 


—_—_—_— 
‘E.V. Appleton. J. Inst. Elec. Engrs. (London), 71, 642 (1932). 
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mission frequency and the quantity sec ¢), and may be estimated within 
the limits of experimental error in most cases. 

Figure 7 gives the frequency to be added or subtracted from the 
maximum usable frequency given by the logarithmic sec ¢, trans. 
mission curves for the cases #0 H+ 0,and H,=H,. The H,=0 
case is uninteresting save for single-hop transmission over the mag. 
netic equator, close to the magnetic meridian. For transmission in the 
continental United States H, is much less than H, and, indeed, is 
negligible over east-west paths, so that such transmission we may con- 


sider as essentially transverse transmission. In this case, 
ti 2 
0 

7272? 


tf "tap 


where 6=0 for the o-component and 1 for the z-component, and the 
logarithmic sec ¢) transmission curves give the correct maximum 
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FREQUENCY MC/S 
Fiaure 6.—Transmission curve, including effect of earth’s magnetic field at level of 
reflection. 
Dotted line =secant-law curve (f’ =f sec ¢s) for the given distance D and wave frequency’. Dashed line= 
Cer’, f) ve for the fied He § this being eet on which, at vertical incidence is a measure of 


the ionization density N. The o-component of the frequency f’ over the given distance is reflected at A, 
and the z-component at B. 


usable frequency for the o-component. For this reason the o-com- 
ponent lies along the sec ¢» axis in the H,=0 case in figure 7. 

For a maximum usable frequency for the o-component below the 
gyrofrequency fz= +/fr’+/,? the z-component always has a maximum 
usable frequency above fg. This is only important for H-layer trans- 
mission and ouiky in cases where the z-component is not too highly 
absorbed at frequencies near fz. This case must not be confused wit 
the case of a transmission frequency f’ less than fz, in which case the 
z-component is reflected from a level above the level where the o-com- 
ponent is reflected. 

Another effect of the anisotropy of the ionosphere due to the earth’s 
field is to cause a difference in the directions of phase and cna a 
pagation in the medium. This results in the wave’s ba ee ted, 
not at the level where the direction of phase propagation is horizontal, 
but where the direction of energy flow (group direction) is horizontal. 
This effect has not been considered in these curves, and is probably 
not > ape to the degree of accuracy to which these calculations 
are carried. 





Ionosphere Measurements 


‘ 
° 


f* f 


3 
z 
° 
= 
< 
< 
a 
W 
" 
> 
1°] 
Zz 
W 


fo 


FREQU 
fe 


1.00 1.01 1.05 1.10 1.20 1.50 3.00 @ 
SECANT $o 


Figure 7. f*— fo’ and f°— fy’ plotted against sec oo f*=frequency whose x-component 
is reflected ai the same level as is the o-component of f°. 
fe’ =f" 008 do 
The value of /” is is given, in Me/s, on each curve. Curves (a) Hi=0. Here /*—fe’=0 and the curves for the 
are indepen on the sec ¢ axis. Curves (0) H.=Hr0. Curves (c) Hr=.0. Here f*—fo' and f*—fy’ 
are in dent of sec ¢o, save near sec ¢o=1. 
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IV. BEHAVIOR OF THE WAVE BELOW THE POINT of 
REFLECTION 


The development thus far has been given for any distribution of 
ionization in the ionosphere below the level of reflection. The effect 
on the reflection of the passage of the waves through the lower iono. 
sphere, and for that matter through the lower atmosphere itself, hag 
been taken care of by using the vertical incidence virtual height in the 
calculation, since this effect enters into the measurement of this 
height. Nothing has been said thus far, however, about phenomena 
that happen below the level of reflection. 

An example of such a complication is given in figure 8. Here the 
transmission curve crosses both the / and the F;, virtual-height curves, 
Therefore, transmission may be expected by way of each layer. The 
angle of incidence is not the same for each layer, since the virtual 
heights of the layers are different. The wave reflected from the 
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FREQUENCY MC/S 
Ficure 8.—Transmission by way of both E and F3 layers. 


A=point of reflection on E layer; B=point of reflection on F; layer; solid line=transmission curve; and 
dashed line=(z,’, f) curve. 


higher layer has a smaller angle of incidence upon the lower layer than 
does the wave reflected from the lower layer. Under certain con- 
ditions, then, it can penetrate the lower layer, and the waves can reach 
the receiver by way of the higher as well as the lower layer. 

When the critical frequency of the lower layer is higher than the 
equivalent vertical-incidence frequency for the wave which would be 
returned from the upper layer, the wave will not penetrate the lower 
layer, but will reach the receiver by way of the lower layer only. 
Furthermore, if this equivalent vertical-incidence frequency is not 
below the critical frequency for the lower layer, but does lie in a region 
of appreciable absorption in this layer, the wave will be returned from 
the higher layer, but will be appreciably absorbed in so doing. Ex- 
perimental observations indicate the absorption of the actual wave to 
be roughly the same as the absorption of a wave of the equivalent 
vertical-incidence frequency measured at vertical incidence. 

It is possible to plot, on the sec ¢» transmission curves, lines cor- 
responding to the values of sec ¢, (see fig. 2) for different heights. 
Such a set of “absorption lines” is shown in figure 9. When the 
transmission curve is superimposed on the (z,’, f) curve the behavior 
of a wave below the reflection level may be estimated by the region 
of the (z,’, f) curve sheet through which the absorption line thro 
the reflection point passes. If this line passes through a region of 





es1OoTitaAt | O6UMCICM TO hCOkkm 


smth Ionosphere Measurements 703 


jbsorption or cuts a lower layer the wave will be absorbed or will not 
netrate through to the higher layer. 
These absorption lines are the lines sec ¢,;=const. for a plane earth 
snd for short distances on a curved earth. They curve toward larger 


VIRTUAL HEIGHT km 


re) 8 
FREQUENCY Mc/S 
Fiaure 9.—Absorption lines on logarithmic transmission curve. 


Dashed line=(z»’, f) curve; dotted lines=absorption lines. EF reflection takes place at point A and F’; 
reflection at point B. reflection at C is shielded by the £Z layer at D. 


' values of sec ¢, for lower heights in the case of greater distances over a 
curved earth, and approach the transmission curve itself for great 
distances. An example of the use of the transmission curves and 
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Figure 10.—Virtual heights at vertical and oblique incidence. 


(a) (z,’,f) curve showing E, F:, and F:layers. (6) oblique-incidence heights corresponding to (a). Dotted 
lines show parts of curves shielded by E layer. Note transmission of some frequencies by several layers. 
This shows how the Fi-layer reflections are relatively unimportant for transmission. (c) (z,’, f) curve 
showing So (X) above critical frequency for the Z layer. (d) oblique-incidence heights corre- 


sponding to (¢). Absorption here (X) completely blocks out some frequencies. All curves are for the 
o-component. 


absorption lines is shown in figure 10, which gives vertical-incidence 
curves and the corresponding frequency-height curves derived there- 
from for the transmitted wave over a distance. From these latter 
curves can be deduced, of course, the equivalent paths or group re- 
tardation of the waves over the transmission path. 
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V. ANGLE OF DEPARTURE AND ARRIVAL OF THE 
WAVES 


Disregarding any possible asymmetry of the wave trajectory due, 
the ots magnetic field, the angles of departure of the waves fry 
the transmitter and of arrival of the waves at the receiver are equi} 
The amount of energy which is radiated from the transmitter gt, 
given angle from the horizontal depends, especially for low angl, 
on the design of the transmitting antenna, and to a great extent », 
the nature of the terrain surrounding the transmitter. Similar facto, 
affect the energy absorbed from the wave at the receiving station, 

For a station located on or near the ground no energy will be radiate) 
below the horizontal, and but little until the angle y of departure ¢ 
the waves above the horizontal becomes appreciable. The minimy 
value of y at which sufficient energy is radiated (or received) tj 
produce a readable signal varies with the terrain, the power of th 
transmitter and the sensitivity of the receiver. Over sea water cm 
be very nearly zero; over land the minimum value may be seven] 
degrees. A fair average approximation may be that y must exceoj 
about 3.5°. 

A simple geometrical calculation gives y in terms of sec ¢p for varioy 
distances, and these may be noted on the transmission curves. Th 
point of reflection must then correspond to an angle y greater thu 
the minimum value assumed for the terrain at the receiver aif 
transmitter. ; 

The effect of this limitation is to limit the maximum distance fu f 
single-reflection transmission to about 1,750 km (over land) for th} 


E-layer and about 3,500 to 4,000 km (over land) for F2-layer tran 
mission. Where the transmission path exceeds these values in | 
calculations must be made on the basis of multireflection transmission, § 

The condition of the part of the ionosphere traversed by the ware 
determines the behavior of the wave. In calculating transmission 
conditions, therefore, the ionosphere around the middle of the path, for 
ee transmission, and around the middle of each jump, for 


multireflection transmission, must be considered. The latter require F 
ment means that usually, for practical calculations, one must consider 
the condition of ionosphere over all the transmission path except a few 
hundred kilometers near the transmitter and receiver. This question 
has been discussed in detail elsewhere.’ It must also be kept in mind 
that at times the o- and z-components are reflected from different 
geographical parts of the ionosphere. 


VI. DETERMINATION OF SEC ¢ 


Figure 11 gives an alignment chart for the rapid determination o 
the factor sec ¢o to be used in calculating the logarithmic transmissio 
curves. To use the chart, place a straightedge so that it pass 
through the desired virtual height and the desired distance laid of 
on the distance scale at the lower left-hand edge of the chart (i- 
creasing distances lie to the left). The ordinate of the intersectio 
of the straightedge with the vertical line corresponding to the same 


’'T, R. Gilliland, 8. 8. Kirby, N. Smith, and 8. E. Reymer. Characteristics of the ionosphere ond thet 
a to ws transmission. J. Research NBS 18, 645 (1937) RP1001; Proc. Inst. Radio Engrs. %, 


The weekly radio broadcasts of the National Bureau of Standards on the ionosphere and radio-transmissoa 
conditions. Letter Circular LC499. 
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desired distance laid off on the main distance scale (increasing 
distances lie to the right) gives the value of sec ¢. The 
relation of the point of intersection to the curved dashed lines of 
| equal ¥ gives the value of the angle of departure of the waves from 
| the horizontal. A point of intersection f above the y=0° line 
indicates an impossible case, where the ray would have to depart 
at an angle below the horizon. 

For example, a distance of 2,400 km and a virtual height of 300 
km corresponds to a sec ¢» of 3.07, and an angle of departure of 8.2°. 


VII. TRANSMISSION FACTORS 


When average transmission conditions over a period of time or 
when a variety of transmission paths are to be considered, or when 
> an estimate of the maximum usable frequencies is to be made without 
a precise knowledge of the ionosphere over the transmission path, it 
is convenient to have available a means by which the maximum 
usable frequencies may be quickly estimated from an approximate 
value of the vertical-incidence critical frequency. 
| The National Bureau of Standards is now beginning a compilation 
| of factors by which the critical frequency for the o-component, 
measured at vertical incidence, may be multiplied in order to obtain 
the maximum usable frequencies. These factors are based on average 
observations over a period of time, and may be applied either to 
average critical frequencies to give average transmission conditions, 
' or to a given observation of a critical frequency to obtain approximate 
' transmission conditions at a given time. 


VIII. CONCLUSIONS 


The type of transmission curves described above have been in 
use at the National Bureau of Standards for the past 2 years in 
studying the correlation of high-frequency radio-transmission con- 
ditions with regular ionosphere observations. The results of con- 
tinuous measurements of high-frequency broadcasting stations and 
observations on other high-frequency signals, as well as the results 
of some specific experiments have been compared with vertical- 
incidence data. Practically all the available data agree with what 
would be expected on the basis of the theory outlined above, and the 
exceptions may in most cases be accounted for. 


Wasuineton, November 22, 1937. 
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A PORTABLE APPARATUS FOR MEASURING VIBRATION 
IN FRESH CONCRETE 


By George L. Pigman, Floyd B. Hornibrook, and Jesse S. Rogers 


ABSTRACT 


A portable vibrograph was designed and constructed for measuring vibration 
in fresh concrete during the process of placing by vibration. An electromagnetic 
type of pickup unit, submerged with definite orientation within the concrete, 
picks up the vibrations and conveys an electric response through either of three 
electric circuits, and then through an amplifier to a cathode-ray oscillograph. 
Particle displacements, velocities, or accelerations can be measured and displace- 
ment, velocity, or acceleration wave forms obtained, depending upon the circuit 
chosen. Frequencies which form simple ratios with 60 cycles can be measured 
accurately, and other frequencies can be estimated. 

The integrating or differentiating circuit-amplifier-oscillograph system was 
found to cause a variation of less than 3 percent in the recording of a constant 
input for frequencies ranging from approximately 20 to 300 c/s. 

The pickup unit was calibrated by use of microscope observations of amplitudes 
of vibration of a steel plate, vibrating in a horizontal direction at approximately 
60 c/s, to which the pickup unit was attached. The amplitudes included in this 
calibration ranged from approximately 0.045 to 0.00035 in. Lower amplitudes 
were evaluated by extrapolation. 

Calibrations of the acceleration and velocity scales were made by replacing the 
pickup unit with an audio-frequency oscillator supplying a sine wave of known 
frequency and voltage, corresponding to calibrated amplitudes. Accelerations 
ranging from approximately 4%.» to 40 times gravity, and velocities ranging from 
0.01 to 2.5 f/s were included in the calibration. 

Typical curves are presented showing the variation of amplitude of vibration in 
concrete with distance from the vibrator. Photographs are given showing corre- 
sponding particle displacement, velocity, and acceleration wave forms produced 
by vibration of concrete with two different types of vibrators. 


CONTENTS 


I. Introduction 

II. Design requirements 

III. Theory and construction of electrical and mechanical systems 
1. General considerations 
2. Pickup unit 
3. Integrating circuit 
4. Differentiating circuit 
5. Indicating system 

IV. Calibration 


U 
I. INTRODUCTION 


Within the past decade, numerous studies have been made of the 
effects of the vibration of fresh concrete. Meager data, however, are 
available on the actual displacements, velocities, and accelerations,' 


LK 
' Throughout this paper the terms amplitude, displacement, velocity, and acceleration refer to particle 


motion unless otherwise specifically stated 
707 
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within the concrete. It was thought that an apparatus capable of 
measuring these vibration characteristics at any point within the cop. 
crete could be of value in numerous studies. Such studies might ip. 
clude the effectiveness of various types and sizes of vibrators, the effect 
of size and rigidity of forms, and the effect of different mixes and 
C/W ratios on the propagation of vibrations within the concrete, 
Accordingly, an apparatus was designed and built at the National 
Bureau of Standards for measuring the foregoing vibration character. 
istics, as well as the frequencies and the propagation velocities of the 
vibrations. 

This paper describes the apparatus in detail and also gives a brief 
description of various types of measurements which have been made, 


II. DESIGN REQUIREMENTS 


The following requirements were considered essential for an ap 
ratus to be used for measuring vibration in fresh concrete over a wide 
range of both field and laboratory conditions. 

1. The unit for picking up the vibrations should be small, light, 
portable, mechanically rugged, and adapted to being immersed to any 
desired depth in the concrete, with a definite orientation. 

2. The indicating device should give a continuous indication which 
can be read rapidly. 

3. The reading should be proportional to the actual displacement, 
velocity, or acceleration, as the case may be, at a given point and ing 
given direction in the concrete. 

4. With proper adjustment, the pickup unit should be capable of 


registering vibrational motion in either a horizontal or the vertical 
direction. 

5. The apparatus should respond accurately over the range of am- 
plitudes ? and frequencies encountered in the commercial vibration of 
concrete. 

6. A means of obtaining wave forms should be provided. 


III. THEORY AND CONSTRUCTION OF ELECTRICAL AND 
MECHANICAL SYSTEMS 


1. GENERAL CONSIDERATIONS 


An electrical method of measuring the vibration in the concrete 
appeared to be most promising in ap abt the requirements. There 
are a number of well-known methods through which vibrational motion 
may be detected by electrical means. The motion may be caused to 
affect an electrical circuit through change of resistance, capacitance, 
or inductance, or through electromagnetic or piezoelectric effects. 
While certain of these methods * ** have been utilized to measure one 
vibration characteristic (either displacement, velocity, or acceleration) 
it is desirable to measure all three characteristics with one detecting 
circult. 

It was found that all three of these vibration characteristics could 
be measured by using a pickup of the electromagnetic type, wht 
has a voltage output proportional to the velocity, and by combining 

* Amplitade refers to the maximum displacement from the position of rest. 
*U. 8. Bureau of Mines Tech. d 


Paper 518 and 556. 
‘CO. D. Greentree, Gen. Elec. Rev. 40, No. 9 (Sept. 1937). 
§M. Mary, Ann. Ponts et Chaussées 338 (March 1936). 
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this pickup with electrical circuits for integrating and for differentiat- 
ing the voltage from the pickup. In addition, the electromagnetic 
type of pickup unit is adapted to the desired simple, compact, and 
ed construction, and requires a minimum of adjustment and re- 
ibration. 
: A similar method of utilizing an electromagnetic detecting system 
in conjunction with electrical circuits to obtain displacements, ve- 
locities, and accelerations has been used by Meyer and Béhm.’ The 
principles of their circuits are woeceagmey ifferent, however, from 
those described in the present paper, a different frequency range was 
covered in their work, and less accurate results were considered satis- 
factory. 
Beale and Stansfield’ have described an electromagnetic pickup 
and a resistance-capacitance type of integrating circuit used in pressure 
studies of diesel-engine operation. Draper and Bentley * have de- 
scribed a similar pickup unit and integrating circuit used in studying 
aircraft vibration during flight. __ Biche ee 
The present apparatus uses a pickup and an integrating circuit 
similar to those used by the above-mentioned investigators. In 
addition, it contains a differentiating circuit, so that by a simple switch- 
ing arrangement either the displacement, velocity, or acceleration as a 
function of time may be shown on the oscillograph screen. 


2. PICKUP UNIT 


The mechanical design of the pickup unit follows the general prin- 
ciples of seismic vibration measuring devices. The unit consists 
essentially of an electromagnet (the seismic element) coupled through 
a spring of low stiffness to the vibrating member, the mass of the 
electromagnet and the spring constant being adjusted so that the 
electromagnet remains relatively motionless throughout the frequency 
range to be covered. The pickup coil is rigidly attached to the 
vibrating member (the case). Hence, the relative motion between 
the electromagnet and the pickup coilicauses an induced,electromotive 
force in the coil which is proportional to the instantaneous velocity 
of the vibrating member. 

The first instrument made was constructed so that with the proper 
adjustment, measurements could be made in eitherga horizontal or 
the vertical direction. During experimentation in the measurement 
of vibration in concrete, it became evident that for the most widely 
used types of vibrators, including the internal and form types which 
vibrate in a horizontal plane, it was sufficient for the instrument to 
respond only to horizontal vibrational motion. Hence a new instru- 
ment was constructed which, by making it suitable for horizontal 
measurement only, was made more compact and rugged. 

his second instrument is shown diagramatically in figure 1. The 
external dimensions are approximately 2% by 2% by 2% in. The 
weight of the pickup unit and supporting rod together is about 2% Ib. 
The bulk specific gravity of the pickup unit is approximately 3.2. 

As shown in figure 1, the suspension of the stationary magnet in 
this instrument consists of two phosphor-bronze strips, about 0.005 in. 


' Elek. Nachr. Tech. 12, 404 (1985). 
1E 160, 667 (1 2 
05, Leno. Belt he 116 Goes): & 8. 281 (1987). 
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in thickness. The width and thickness of the middle section of eagh 
strip was decreased to obtain a natural period of the magnet of about 
% second. A phosphor-bronze leaf spring was adjusted to exert g 
light pressure upon the surface of the magnet in order to produce 
sufficient damping to minimize the effect of transient vibration. 

A semiflexible joint between the pickup-unit case and the supporting 
rod (section A-—A’, fig. 1) allows free movement of the unit in the 
direction of measurement, but simultaneously provides high strength 
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FiaurE 1.—Schematic drawing of pickup unit. 





to allow the unit to be pushed, twisted, or pulled in any manner to 
permit rapid placement and alignment in the desired position in the 
concrete. The instrument is watertight, but its internal parts are 
readily accessible for any adjustment or repair. 


3. INTEGRATING CIRCUIT 


The integrating circuit * is shown in figure 2. The pickup coil, P 
moving in a uniform magnetic field, has at any instant an induced 
electromotive force which is proportional to the pickup coil velocity, 
that is, e=K(dz/dt), where z is the displacement of the coil from its 


* For discussions of the theory involved, see Beale and Stansfield, Engineer 160, 667 (1935); also Draper and 
Bently, J. Aero. Sci. 3, 116 (1936); 4, 281 (1937). 
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osition of rest and Kis aconstant. For a single sinusoidal vibration 
of frequency «/2x the current in the circuit will be 


. é 
’="s 


VR?+ (@L—1]a0)” 


where FR is the resistance, in ohms, in series with the condenser of ca- 
B pacitance C farads, and L is the inductance in henries of the pickup 


oT err T © 
Pick-up Cc 


Coil T %6 amplifier input 


(1) 








P 





. ve) 
FiaurE 2.—Diagram of integrating circutt. 


coil. If Rand Care increased and Z decreased sufficiently to make the 
term (wL—1/wC)? small compared with R?, then approximately, i= 
¢/R=K (dx/dt)/R. . es 7 : 

In a complex wave the circuit introduces a phase distortion, due to 
the difference in phase lead of the various harmonic components. 
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Figure 3.—Diagram of integrating and amplifying circuits. 


Temporarily neglecting this phase distortion, the voltage across the 
condenser terminals at any instant, ¢, is then 


pul (*... K (de, K 
B= | idt= Fp [Sp tt= Zero (2) 


that is, EZ, is proportional to the displacement at the time, ¢. If E, 
then is applied t ugh Be amplifier to the deflecting plates of an 
oscillograph having a linear time axis, the trace on the oscillograph 





712 Journal of Research of the National Bureau of Standards yy, 


screen will give an accurate picture of the displacement of the Pickup 
as i. beegren of time. i ” ‘ 

e integrating and amplifying circuits are shown in figure 
The value of C used in the integrating circuit is 1 pf; two values tf 
are included, 10,000 and 200,000 ohms, respectively. Using th, 
resistance of 10,000 ohms, making the time constant (CR) equal t, 
0.01, the sensitivity is approximately 20 times that when the resistang, 
R equals 200,000 ohms (CR=0.20). The amplifier has three regis, 
ae coupled stages using pentode tubes having high 
gain. The maximum usable voltage gain through the amplifier ; 
35,000, above which the noise level becomes excessive. 

The effect of frequency change on the output of the complet 
integrating-amplifying circuit is shown in figure 4. To obtain this 


10 100 
Frequericy, cycles per second 


Figure 4.—Effect of variation in ee on the amplitude reading for constan : 
amplitude input. i 


figure a beat-frequency oscillator was used as a voltage source ands | 
constant ratio of input voltage to frequency was maintained, becaux 
of the integrating action, to correspond to constant amplitude input. 
The variation in amplitude reading on the oscillograph was 
observed. It can be seen that, with a time constant of 0.20 for the 
integrating circuit, the frequencies can vary from about 20 to 300 
c/s without introducing an error greater than 3 percent in the oseill- 
graph reading. Similarly, with a time constant of 0.01 for the inte 
grating circuit, frequencies can vary from about 40 to 300 c/s without 
causing more than the same error. It should be noted that an ero 
of 3 percent is about the least error that can be realized in reading the 
oscillograph even under ideal conditions. 

Above 300 c/s the error increases, reaching 1 a agi te: at 1,000 cycles. 
In commercial vibration of concrete, practically all vibrators operate 
in the frequency range of 50 to 150 ods and it is probable that most 

ignificant harmonics occur below 1,000 c/s. 5 

The phase lead for a sinusoidal current of frequency w/2m is (= 
tan’ 1/wRC, and represents the maximum phase error for aly 
harmonic. Substituting the circuit constants in this expression, It! 
found that the phase distortion for the circuit with the time constant 
of 0.20 is negligible. In this case the angle 6 is about 2.3 degrees st 
20 cycles and decreases at higher frequencies. The circuit with the 
time constant of 0.01 (designed to be used only at frequencies above 
100 c/s) gives a value of 9 degrees for the angle @ at 100 cycles, Even 
this maximum distortion, representing a shift of 1/40 wave length, # 
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» wherein ¢; is an alternating voltage o 


| sinusoidal wave may be con- 
| sidered as representing a 
» single term of a Fourier’s 4 
’ series giving the actual wave. 
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normally insignificant, particularly as the higher harmonics are gen- 
erally of comparatively low magnitude in the original wave and con- 
tribute even tse proportionally to the shape of the integrated wave. 


4. DIFFERENTIATING CIRCUIT 


Several different electrical circuits may be used to perform what cor- 
responds to a mathematical differentiating process and thus to convert 
an alternating voltage representing a velocity-time wave into its cor- 
responding acceleration-time wave. The voltage measured across an 
inductor of low impedance, in series with a resistor of high resistance, 
is approximately the time derivative of the current in the circuit. 
Similarly, the voltage across a resistor of a low resistance in series 
with a capacitor of high impedance is approximately the time deriva- 
tive of the voltage across the circuit. However, the circuit adopted 
as most satisfactory for the present purpose is shown in figure 5, 

i feeined from the pickup unit, 
and é, is the voltage output after passing through the differentiating 
circuit. For convenience, it may be assumed that e,—£ sin wt, where 
Eis a constant, w=2zf, f is 
the frequency in cycles per RF 


F second, and ¢ is the time 
) in seconds. For a complex 


wave in a particular case this mm 


,.) 








The circuit shown in figure Fygure 5.—Diagram of differentiating circuit. 


' 5 is essentially an air-core 


transformer, with a high resistance in series with the primary. With 
the alternating voltage e,=E sin wt placed across the input, the cur- 
rent i in the primary circuit, with a small current in the secondary, is 
given by the equation 


E sin wt=iR+L9 (3) 


where F is the total resistance, and L is the inductance, in henries, of 

the primar circuit. 

‘ With R large with respect to Lw, the voltage across the secondary 
ecomes 


a=-p Ew cos wt, (4) 


or a-y Fi Kee (from eq 2) (5) 


where d*z/dé is the acceleration at the time, ¢, and M is the mutual 
inductance, in henries, of the primary and secondary coils. 
The phase distortion due to the difference in phase lag ¢=tan™! 

“ait is negligible with the given constants for all significant har- 
cs. 

_ The present circuit has an advantage over the other differentiatin 

circuits mentioned, in that the relation between the resistance an 
57763—88——10 
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the reactance of the coil in the input circuit is of small consequence 
a condition which becomes of greater importance at low frequencies 
Also, while a large voltage loss occurs in any of the three circuits, , 
portion of the loss may be regained in the present one by increagj 
the amount of winding on the secondary coil. The amount which 
the number of secondary turns may be increased is limited by the 
size of the coil and the amount of shielding required. 

A diagram of the differentiating circuit with the constants ani 


the switching arrangement used, is shown in figure 6. For the 
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Figure 6.—Diagram of the acceleration and velocity measuring circuits. 











constants chosen, as shown in the diagram, and with the input volt- 
age corresponding to constant acceleration, the relation between 
acceleration reading on the oscillograph and the frequency follows 
the curve shown for the amplifier alone in figure 4. 

Included in figure 6 is a voltage attenuating circuit with constants 
adjusted so that the desired range of velocity readings may be 
obtained to correspond to the displacement and acceleration scales. 


5. INDICATING SYSTEM 


A cathode-ray oscillograph with 5-in. screen serves as an indicating 
device, the screen being ruled in 0.10-in. divisions to permit rapid 
readings of the beam deflection. The oscillograph also contains a 
single-stage amplifier, which gives a voltage gain of about 25, making 
the total available voltage gain in the apparatus approximately 
900,000. 

For precise results, such high amplification might be undesirable. 
In the present case, variations of 5 percent are not considered ex- 
cessive, and, with a simple rapid method of checking the calibration 
at any time having been devised, described in a later section, 
magnitude of amplification used is not objectionable. 
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The oscillograph may be utilized to measure frequency of vibration, 
although it is often simpler to measure the eapeney of the vibrator 
directly with some form of frequency meter. If the frequency to be 
measured, however, forms a simple ratio with 60 cycles, it can be 
determined directly through the patterns it forms on the oscillograph 
screen when the sweep circuit has been synchronized with the 60-cycle 
power supply. The accurate determination of any other frequency 
by means of the oscillograph required the use of a source of alternating 
voltage of known variable frequency. ‘ 

Figure 7 shows the complete apparatus, assembled as it is used in 
making measurements. 


IV. CALIBRATION 


Calibration of the apparatus for displacement measurement was 
accomplished through the use of a steel plate (approximately 24 by 
5 by % in.) held in a vertical plane by being clamped at both ends 
to heavy concrete blocks. The pickup unit was bolted to the center 
of the plate; and the plate caused to vibrate horizontally by a rotating, 
motor-driven eccentric mass also attached to the center of the plate. 
The amplitude of vibration of the pickup unit was observed through 
a microscope, having a calibrated reticule in the eyepiece, sighted 
upon a fine crosshair on the top of the pickup unit. It was found 
that with the ends of the steel plate held in this manner, and the 
rotating eccentric driven at frequencies from 50 to 70 c/s, the vibration 
obtained was smooth and linear in character, and in general could be 
held to a very constant value. The amplitudes obtained for the 
calibration in this manner ranged from 0.045 in. to 0.00035 in. 

During the calibration process, the various resistors in the tapped 
gain-control of the amplifier were adjusted in value to give the 
desired beam deflections on the oscillograph for the chosen amplitude 
ranges. A 3-in. total deflection of the oscillograph beam was obtained 
for amplitudes of approximately 0.045, 0.015, 0.0045, and 0.0015 in., 
respectively, with 200,000 ohms in the integrating circuit, and for 
amplitudes of approximately 0.0020 and 0.0010 in., respectively, 
with 10,000 ohms in the integrating circuit. For each range a linear 
relation was obtained between the optically measured amplitude 
and the oscillograph beam deflection. 

Although in this scale calibration the frequency was maintained 
between 50 and 70 c/s, other tests were made at constant amplitude 
in which the frequency was varied from approximately 30 to 100 c/s 
(the frequency limits of the available vibration source) to check-the 
linearity of the pickup unit output. No deviation from linearity 
was found in that range. 

In order to utilize all of the available amplification, and to increase 
the range of the apparatus for the measurement of amplitudes too 
small to be conveniently observed by optical methods, an extrapola- 
tion was made. Observations were taken of the input voltages, at 
60 c/s, necessary to give a definite oscillograph-beam deflection as the 
controls were set, first for each of the calibrated amplitude ranges, and 
then for each of the uncalibrated ranges. Assuming that the straight 
line relation found for the ranges covered by direct calibration still 
holds, the calibration can be extended for the evaluation of lower 
amplitudes. In this manner four additional ranges were obtained, 
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the lowest of which gives a 3-in. total oscillograph deflection for ap 
amplitude of approximately 0.00002 in. 

he calibration of the acceleration and velocity measuring circuits 
was accomplished through the use of a beat-frequency audio oscillator 
with an output wave containing a harmonic content of approximately 
0.5 percent. With a constant voltage output from the oscillator, at 
a known frequency, the signal was applied to the calibrated integrating 
circuit and then to the uncalibrated acceleration and velocity measur. 
ing circuits, respectively, and corresponding readings were obtained 
under the three conditions. The amplitude corresponding to the 
given input signal was known from the first reading, and, with the 
signal being a sine wave of known frequency, the corresponding maxi- 
mum velocity and acceleration could be accurately computed. 

The acceleration scale was calibrated for six arbitrary ranges, 
covering accelerations from 2 to 1,250 f/s/s, or approximately ¥, 
to 40 times gravity. The velocity scale was calibrated for five ranges, 
covering velocities from 0.01 to 2.5 f/s. 

In order to detect any changes in calibration which might appear, 
due to variations in the integrating, amplifying, or oscillograph cir. 
cuits, a portable calibration checking device was constructed. This 
consists of an attenuator, supplied with 60-cycle alternating current, 
with which a known voltage may be supplied to the input of the 
integrating circuit. Immediately after calibration in each range, 
oscillograph beam deflection readings were taken with this attenuator 
in the circuit. A volume control, incorporated in the amplifier, per- 
mits the amplification to be adjusted at any time to obtain these 
standard readings. With this device, the calibration may be rapidly 
checked for any range before that range is used. 

In the calibration of the integrating circuit, observations were made 
of the oscillograph reading obtained with the pickup unit being 
vibrated with large implitilie at right angles to the direction in which 
the unit is designed to respond. It was found that even such | 
transverse vibration apparently had no effect on the reading. 
concrete, the pickup unit should thus respond only to vibrational 
motion in the direction of its axis. 


V. VIBRATION MEASUREMENTS 


A short laboratory study was made to test the adaptability of the 
vibrograph to the purposes for which it was designed. ‘Two commer 
cial internal vibrators,’ of rotating eccentric weight type, the outside 
diameters of the vibrating shafts being 1 and 2% in., repose 
were used for vibrating the concrete. The small vibrator is directly 
driven by a %-hp electric motor, and operates at frequencies from 
9,000 to 14,000 rpm. This vibrator is of a type comgaed to compas 

e 


laboratory specimens, and was not inten or field work. Th 
large vibrator is hydraulically operated, being driven by a 3-hp electric 
motor and can be operated at any frequency below 6,500 Foe _ Th 
amplitudes of vibration of the two a woe a A when suspended in air 
were approximately 0.012 and 0.050 in., respectively. There was no 
appreciable diminution in vibrator amplitude when either vibrator was 
immersed in concrete, according to measurements made just above 
the immersed section of the vibrator. 


© The vibrators were generously lent to the Bureau for these studies by the Viber Co., Ltd., and by the 
Klectric Tamper and Equipment Co. 
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For most of the study a clay concrete was used, proportioned 
1:2.2:3.8 parts by weight of clay, sand, and graded aggregate, 
respectively. The form used was 4 by 4 by 2 ft. high, inside dimen- 
sions, having rigid 4-in. concrete walls and base. During vibration of 
the concrete, ost no vibration was perceptible in the walls of the 
form. 

With the limited data obtained, no attempt was made to evaluate 
the effects upon amplitude distribution in the concrete of such variables 
as size and shape of vibrator, size, shape, and rigidity of form, the 
frequency of vibration and variations in mix. However, several typical 
curves are presented in which some of these factors are varied. 

Figure 8 shows the variation of amplitude of horizontal vibration 
at different depths in the concrete at various distances from the 


surface of 
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FiaurE 8.— Variation of displacement amplitudes of vibration with depth in concrete 
at various distances from the vibrator. 


vibrator. The latter;was held{in the same position throughout (near 
the midpoint of the concrete) while the pickup unit was moved along 
the diagonal of the form. This figure indicates the necessity for care- 
fully controlling the depth of immersion of the pickup unit during 
measurements, especially near the vibrator. 

Figure 9 shows the variation of horizontal amplitude with distance 
from the vibrator for two different speeds of the large vibrator and 
one speed of the small vibrator. In No fe these data, the vibrator 
was moved towards the pickup unit, the latter being maintained in 
one position at approximately the midpoint of the immersed section 
of the vibrator. 

, The accuracy of the points shown in figure 9 depends upon the degree 
in which the pickup unit will follow the actual motion of the concrete 
which surrounds it. This in turn depends upon (1) the ratio of the 
wave length of the vibratory wave in the concrete to the size of the 


Pickup unit, (2) the change in amplitude in the concrete over the 
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length of the pickup unit, (3) the density of the pickup unit, and (4) 
the ease with which the concrete will flow around the unit, or the 
“internal friction” of the concrete. It is probable that in this respect 
the pickup unit acts as a piece of aggregate of the same approximate 
dimensions. 

It will be shown later that with the frequencies ordinarily used in 
vibrating concrete the wave length is large with respect to the length 
of the pickup unit. The average specific gravity of the entire pickup 
unit is about 3.2, which is somewhat higher than the average value 


.006 


LEGEND 
© Large vibrator, 5700 rom 
+ ” «~~, $400 rprn 
x —? /2000 


%. 
v 
2 
Lv. 
2 
S 
y 
=, 


0 x 


0 / b vvddiccl &b 6 " ee 9, As: fi 4 
Distance from edge of vibrafor, inches 


Figure 9.—Variation of displacement amplitude with distance from vibrator for 
clay concrete contained in rigid concrete form. 


for concrete. However, as the seismic element remains nearly sta- 
tionary, the mass being accelerated is less than the total mass of the 
pickup unit, and unduly large inertia forces are not expected to be 
present. The “internal friction” of concrete is considerable, even with 
relatively large amounts of water present, and in the dry concretes 
generally best adapted to placement by vibration is probably partic- 
ularly great in magnitude. The amplitude-distance curves given In 
figure 9 for the large vibrator show that for distances from the source 
of vibration greater than 4 or 5 in. the amplitude gradient is not 
excessively large. Considering all of these factors it may be ass 

that beyond 4 or 5 in. the pickup unit registers some average value of 
the amplitudes distributed over its surface. For distances less than 
about 4 in. the absolute value of the amplitude reading is to be ques 
tioned. However, from consideration of the regularity of the poms 
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Figure 11.—Corresponding displacement, velocity, and acceleration wave forms, 
A, B, and C, respectively, with correct phase relations maintained, of vibration 
produced in concrete with a compressed-air vibrator. 
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FiguRE 12.—Displacement, velocity, and acceleration wave forms, A, B, and (, 


respectively, of vibration in concrete produced by a revolving eccentric-weight type 
vibrator. 
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in figure 9 and their relation to the amplitude of the freely suspended 
vibrator, it is probable that even in this region the readings give a fair 
approximation of the true amplitudes in the concrete. 

Figure 10 shows the variation of horizontal amplitude in portland 
cement concrete contained in a heavy wooden form, 2% by 4% ft. by 
2 ft. high inside dimensions. In this case, although the vibrator was 
moved along the middle line of the box, the form itself vibrated very 
perceptibly, and the concrete was affected much more by the vibrator 
than was the clay concrete contained in the rigid concrete form (see 
fig. 9). (The clay concrete in the rigid form appeared, in general, to 
be compacted only 4 to 7 in. from the vibrator; the cement concrete in 
the wooden form was nearly all compacted in a short time.) 

For the purpose of obtaining typical wave forms of several types of 
concrete vibrators, and to demonstrate qualitatively the actions of 
the integrating and differentiating circuits, records were made of the 
oscillograph traces obtained with the pickup unit immersed in clay 
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Figure 10.—Variation of displacement amplitude with distance in portland cement 
concrete contained in heavy wooden form. 


concrete to which vibration was being applied. Figure 11 shows 
wave forms obtained with a reciprocating compressed-air vibrator 
applied to the form, A, B, and C being respectively the displacement, 
velocity, and acceleration waves. The correct. phase relation was 
maintained between the photographs by synchronizing the oscillograph 
sweep circuit with the vibrator frequency by means of a second pickup 
unit. The fundamental wave of 4 has a frequency of approximately 
110 c/s. ° A close inspection of the three waves will reveal that, going 
from A to B to C in that order, for the same abscissa the ordinate of 
each succeeding wave is proportional to the slope of the preceding 
one. Thus the maximum peak of velocity wave B corresponds pre- 
cisely to the region of maximum slope of the displacement wave A. 

As a further check on the validity of the relation indicated in the 
— of the three curves, the areas under the acceleration and 
Velocity curves were measured with a planimeter and the corresponding 
Velocity and displacement wave forms plotted (with arbitrary adjust- 
ment of scale). The plotted wave forms agreed closely with the 
pictured forms, although lack of linearity of the time base (oscillograph 
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sweep circuit) prevented precise superposition of the plotted waye m 
the picture. 

The predominant part played in the velocity and acceleratiy, 
waves by the high-frequency components, which are relatiy, 
unimportant in the displacement wave, is to be particularly no 
and appears to be characteristic of the vibration produced by pistop. 
type vibrators. 

igure 12 gives the corresponding wave forms obtained with th 
large internal vibrator described in an earlier section, Speetae at a 
frequency of approximately 3,600 rpm. As before, A, B, and Cap 
respectively the displacement, velocity, and acceleration waves, by; 
in this case the correct phase relation was not maintained betweg 
the three traces. It may be noted that, although the displaceme; 
wave is superficially sinusoidal, high-frequency components, again, 
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Fiaure 13.—Variation of the velocity of propagation of vibration through concréle, 

proportioned 1:2.2:3.8 by weight, with change of water content of the concréle. 


are very evident in the acceleration wave form. The maximun 
acceleration associated with the harmonic waves alone, shown i 
figure 12 (A), is of the same order of magnitude as the maximum 
acceleration associated with the fundamental wave alone. It seems 
probable, therefore, that attempts to correlate efficacy of vibration 
with acceleration may yield misleading results when the values for 
acceleration are obtained by computations from amplitude and 
frequency data under the assumption that no harmonics are preset. 
Unless the absence of appreciable harmonics is demonstrated (i. ¢., 4 
truly sinusoidal displacement wave obtained), the possible effects of 
the harmonics should be considered. 


Another type of measurement that can be made with this instrument J 


is illustrated by figure 13, which shows the variation of veloci of 
propagation of vibration with change in the water content of portal 
cement concrete proportioned 1:2.2:3.8 by weight. When two pickup 
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units, one at the source of vibration and one several feet away in the 
concrete are used, and the oscillograph sweep is synchronized with the 
closest pickup unit, then as the second pickup unit is moved toward or 
from the vibrator, a phase shift occurs in the oscillograph pattern which 
is equal to the phase shift within the concrete over the distance moved. 
From this the wave length in the concrete can be determined. Know- 
ing the frequency f and the wave length }, the velocity can be calcu- 
lated according to the relation v=/x. 

The technique used in making these particular measurements was 
not refined and accordingly the values given are first approximations 


nl as 
: the trend indicated in figure 13 was found to hold for concretes 
having much greater water content than shown, and measurements of 
the velocity in water gave results within 15 percent of the accepted 
value of the velocity of sound in water. 

For the velocities indicated in figure 13 the wave lengths varied 
between 8 and 16 in., the frequency being 9,000 rpm. At a frequency 
of 3,600 rpm, the wave length for the same concrete would be from 
20 to 40 in. Thus, as mentioned in an earlier section, the length of 
the pickup unit is small compared with the wave length, for the driving 
frequencies encountered in most of the commercial vibrators now in 
; use. 

VI. CONCLUSIONS 


An electromagnetic pickup unit, and amplifier, and special electrical 
integrating and differentiating circuits used in conjunction with a 
cathode-ray oscillograph for measuring instantaneous particle veloci- 
ties, displacements, and accelerations, for measuring the frequencies 
and for showing the wave forms of mechanical vibediitins in fresh 
concrete are described. It has been shown that the apparatus will 
give a measure of these characteristics with sufficient accuracy for 
conditions found in present-day concrete-vibration practice, and, in 
addition, is adapted to a variety of more precise laboratory tests. 


Appreciation is expressed to members of the Bureau’s sound section, 
particularly to W. F’. Snyder, for their assistance in the development 
and calibration of some of the electrical circuits. 

Wasuineton, March 4, 1938. 


"| Wave velocities as low as 260 f/s in peat and 360 f/s in fine sand have been reported. R. Koehler and A. 
Ramspeck, Verdffent Inst. Deut. Forschungsges. Bodenmechanik. 4, 9 (1936) . 
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